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REFLECTOR  PATTERN  CONTROL  USING  VARIABLE  PHASE  PLATES 


Daniel  Jacavanco,  Sr. 

EM  Technology  Applications  Section 
EM  Techniques  Branch 
Electromagnetic  Sciences  Division 
Rome  Air  Development  Center 
Hanscom  AFB,  MA  01731 

ABSTRACT 

Results  are  presented  for  two  applications  of  a  novel  technique 
for  pattern  control  in  reflector  antenna  systems.  In  the  first, 
suppression  of  a  high  power  mobile  S-band  source  was  accomplished 
while  maintaining  a  continuous  link  with  a  desired  low  power  signal, 
fixed  in  position.  This  resulted  in  the  proper  operation  of  a  phase 
sensitive  receiver  to  lock  onto  the  desired  signal,  and  not  the 
stronger  interfering  signal  during  antenna  experiments. 

In  the  second  application,  a  direct  measure  of  the  angular 
variation  of  ground  clutter  at  the  RADC  antenna  test  facility  was 
accomplished  using  the  variable  phase  plate  technique.  This  new 
measurement  approach  offers  a  significant  improvement  over  the 
traditional  'field  probe'  method  for  characterizing  the  quiet  zone 
of  an  antenna  mount.  By  nulling  out  the  direct  ray  and  using  a  high 
gain  paraboloid,  local  scattering  centers  have  been  Identified  with 
greater  angular  resolution  than  the  field  probe  method.  New  scat- 
terers  have  been  identified  that  have  the  potential  to  seriously 
degrade  low  sidelobe  pattern  measurements  in  the  vertical  plane. 


1.  INTRODUCTION 


Preliminary  results  of  a  new  technique  for  pattern  control  in 
reflector  systems  have  been  reported  by  Jacavanco^ • Two  AFIT 
Masters  Theses  have  been  written  on  the  theoretical  aspects  of  the 
problem  -  Havens^  and  Rudisill^.  The  essence  of  the  technique  is 
depicted  in  Figure  1.  Two  small  aluminum  disks,  attached  to  rods 
are  held  in  place  on  the  surface  of  a  reflector  by  collet  chucks. 
Pattern  control  is  accomplished  by  adjusting  the  height  of  the  disks 
above  the  surface,  thereby  modifying  the  amplitude  and  phase  distri¬ 
bution  in  the  aperture  plane  of  the  reflector. 

Figure  2  shows  the  small  effect  the  placement  of  two  eight  inch 
diameter  disks  flush  with  the  surface  of  the  reflector  have  upon  the 
antenna  pattern.  Figure  3  shows  a  measured  antenna  pattern  after  a 
null  has  been  created  at  nine  degrees.  The  large  dot  is  the  orig¬ 
inal  level  of  the  antenna  response  at  that  angle  before  nulling. 
Cancellation  of  20  dB  is  obtained  with  little  degradation  to  the 
overall  antenna  performance. 

The  purpose  of  this  paper  is  to  describe  two  practical 
applications  of  this  novel  technique  for  pattern  control. 

2.  JAMMED  REFERENCE  ANTENNA 

In  order  to  measure  the  RF  performance  of  low  sldelobe 
antennas,  the  RADC  Electromagnetic  Test  Facility  uses  the 
Scientific  Atlanta  2020  Antenna  Measurement  System.  The  heart  of 
the  S/A  2020  system  is  a  phase-locked,  narrow  band,  high  dynamic 


range  receiver.  This  receiver  is  required  when  antennas  are 
measured  that  possess  very  low  sidelobes.  In  this  case,  the  dynamic 
range  of  the  receiver  has  to  be  much  greater  than  the  ratio  of  the 
main  beam  to  average  sidelobe  level  of  the  antenna  under  test  in 
order  to  measure  the  depth  of  the  pattern  nulls.  The  phase  refer¬ 
ence  for  the  receiver  is  obtained  through  a  separate  reference 
signal  antenna. 

The  following  problem  exists  when  adaptive  phased  array 
antennas  are  tested:  During  the  test  of  the  adaptive  performance, 
two  signals  are  presented  to  the  low  gain  reference  antenna  -  the 
desired  signal  and  the  jammer  signal.  The  signals  come  from  two 
different  sources  separated  in  angle,  but  at  the  same  frequency. 
However,  the  jammer  signal  power  level  is  usually  50  dB  greater  than 
the  desired  signal.  The  receiver  locks  on  to  the  stronger  signal 
(the  jammer),  loses  phase  lock  with  the  desired  signal  and  becomes 
inoperable  for  measuring  the  amount  of  desired  signal  power  received 
by  the  antenna  under  test. 

The  reflector-disk  pattern  control  method  provides  the 
isolation  required  between  the  desired  signal  and  the  jammer  at 
minimal  cost,  and  in  no  way  affects  the  operation  of  the  system  for 
the  measurement  of  non-adaptive  antennas.  The  horn  antenna  which 
receives  the  reference  signal  has  been  replaced  by  a  high  gain 
reflector  antenna  mounted  on  the  roof  of  the  measurement  facility 
and  pointed  towards  the  transmitter  building.  Two  small  moveable 


disks,  mounted  on  the  surface  of  tlie  dish  are  adjusted  In  position 
until  the  jammer  signal  is  cancelled.  What  follows  is  a  description 
of  how  this  was  achieved  manually  with  a  description  of  the  test 
results . 

A  six  foot  diameter  S-band  reflector,  attached  to  a  short  sec¬ 
tion  of  a  metal  tower  was  mounted  on  the  roof  of  the  main  building. 
The  antenna,  with  a  suitable  feed  horn,  was  pointed  towards  the 
transmitter  building.  An  S-band  link  was  established  at  3.4  GHz. 

See  Figure  4  for  the  plan  view  of  the  test  range.  A  jammer-equipped 
truck  was  outfitted  with  an  S-band  signal  source,  TWT ,  six  foot  re¬ 
flector  and  a/c  generator.  The  truck  was  deployed  on  the  road  to  the 
transmitter  building  at  approximately  the  same  depression  angle  used 
in  earlier  tests  on  an  80  element  linear  array  adaptive  antenna. 

The  first  step  was  to  measure  the  relative  strength  of  the  two 
signals  reaching  the  main  building.  For  this  measurement  they  were 
separated  slightly  in  frequency  so  that  they  could  be  viewed  simul¬ 
taneously  by  a  spectrum  analyzer.  A  small  S-band  horn  with  a  main- 
beam  wide  enough  to  include  the  desired  and  interfering  sources  was 
placed  at  the  same  position  as  the  S-band  dish.  Both  signals  were 
viewed  with  a  Tektronix  Spectrum  Analyzer.  A  precision  step  attenu¬ 
ator,  mounted  in  the  coaxial  line  between  the  small  horn  and  the 
spectrum  analyzer  was  used  to  compare  the  two  signals.  It  required 
50  dB  of  attenuation  to  reduce  the  jammer  signal  level  to  the  level 
the  desired  signal.  See  Figure  5. 


The  next  step  was  to  measure  the  relative  stengths  of  the  two 
signals  when  both  are  received  by  the  six  foot  reflector  antenna 
with  the  two  aluminum  disks  mounted  on  the  surface.  Again  the  step 
attenuator  was  Inserted  In  the  coaxial  line  between  the  feed  horn 
and  the  spectrum  analyzer.  The  paraboloid  was  boreslghted  in  order 
to  be  sure  that  the  malnbeam  and  not  a  sldelobe  Intersected  the 
transmitter  building.  Figure  6  shows  the  relative  strengths  of  the 
signal  and  jammer  when  received  by  the  dish.  As  expected,  the 
-20  dB  dB  sldelobes  of  the  dish  reduced  the  strength  of  the  jammer 
relative  to  the  desired  signal  to  -30  dB. 

Finally,  the  two  8-inch  diameter  aluminum  disks  were  adjusted 
until  the  jammer  signal  was  cancelled.  This  was  performed  manually 
while  viewing  the  output  displayed  on  the  spectrum  analyzer.  The 
first  30  dB  of  cancellation  was  not  difficult  to  obtain  -  the  proper 
position  of  the  disks  was  found  by  carefully  moving  each  one  indiv¬ 
idually  until  the  received  power  of  the  jammer  was  minimized.  The 
fine  tuning  of  the  disk  positions  to  obtain  the  last  10  dB  of  can¬ 
cellation  was  more  difficult.  There  was  no  effect  noted  upon  the 
desired  signal.  Using  the  step  attenuator,  the  amount  of  cancell¬ 
ation  achieved  by  the  final  placement  of  the  disks  was  measured  to 
be  ^0  -  43  dB.  See  Figure  7.  Thus  it  was  was  demonstrated  that 
an  unjammed  reference  signal  (desired)  was  made  available  for  the 
2020  Antenna  Measurement  System. 

Present  plans  include  the  Installation  of  computer  controlled 


stepping  motors  to  automatically  control  the  position  of  the  disks. 
When  driven  by  a  gradient  search  algorithm  that  minimizes  the  total 
power  output  of  the  reference  antenna,  antenna  measurements  over  a 
broad  range  of  frequencies  will  be  possible. 

For  the  purposes  of  display  only,  the  jammer  frequency  was 
slightly  offset  from  the  desired  signal.  In  practice,  the  two 
signals  would  be  at  the  same  frequency.  In  order  to  overcome  the 
limitation  of  a  minimum  total  power  output  algorithm  (l.e.  jammer 
power  equals  desired  power),  one  could  simply  shut  off  the  desired 
signal  while  the  nulling  of  the  jammer  signal  is  accomplished. 

3.  QUIET  ZONE  MEASUREMENT 

The  traditional  method  of  measuring  the  "quiet  zone"  of  an 
antenna  test  mount  is  to  plot  the  output  of  a  modest  sized  probe 
antenna  as  it  traverses  a  line  perpendicular  to  the  direction  to 
the  transmitting  antenna  in  a  far-fleld  range.  From  the  plot  of 
the  variation  in  amplitude  with  traverse  distance,  one  can  calculate 
the  strength  of  reflected  waves  from  the  maxima  and  minima  of  the 
interference  pattern  of  the  trace.  The  direction  of  arrival  of  the 
strongest  reflection  is  Inferred  from  the  periodicity  of  the  pattern. 
The  accuracy  of  this  measurement  of  free  space  standing  waves  is 
usually  limited  to  a  few  tenths  of  a  dB  variation,  corresponding  to 
a  quiet  zone  of  approximately  -35  dB. 

Figure  8  shows  why  reflections  from  the  ground  must  be  deter¬ 
mined  and  eliminated  (if  possible)  in  an  antenna  test  range.  It 


Is  a  plot  of  the  power  received  by  a  four  foot  diameter  paraboloid 
taken  In  the  vertical  elevation  plane  at  3.1  GHz.  The  "sky  side" 
of  the  plot  is  when  the  dish  is  pointed  above  the  transmitter, 
while  the  "ground  side"  of  the  plot  is  when  the  dish  looks  into  the 
ground.  The  disparity  between  the  two  halves  is  obviously  due  to 
ground  reflections.  For  the  purpose  of  this  discussion,  the  pattern 
has  been  smoothed  out  by  plotting  the  average  of  the  trace  over 
those  angles  where  the  ground  reflection  caused  excessive  pen 
fluctuation  (-12  degrees  to  -24  degrees). 

Figure  9  is  a  schematic  diagram  showing  the  relationship 
between  the  direct  ray  and  reflected  rays  during  an  antenna  pattern 
measurement.  Typically  the  transmitting  antenna  has  average  aide- 
lobes  of  -20  dB  relative  to  the  main  beam.  M.so ,  for  the  terrain  In 
question,  the  forward  scatter  of  the  ground  bounce  is  in  the  -20  dB 
range  also.  This  means  that  the  reflected  wave  In  Figure  9  Is  of 
the  order  of  -40  dB  relative  to  the  main  beam.  Franchl  and  Poles^, 
have  measured  the  quiet  zone  of  the  mount  in  question  to  be  -35  dB 
over  a  12  foot  dimension. 

For  the  following  experiment,  a  four  foot  diameter  paraboloid 
was  equipped  with  the  two-disk  nulling  device  described  by 
Jacavanco^  and  full  vertical  plane  measurements  of  the  pattern  were 
made  from  +45  degrees  above  Che  horizon  Co  -45  degrees  below  the 
local  horizon.  (Mount  limitations  prevented  further  movement). 

When  the  antenna  was  pointed  above  the  horizon,  two  jammers 


> 
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were  effectively  present.  However,  since  the  direct  ray  is  stronger 
than  the  reflected  by  40  dB,  it  is  considered  the  effective  janmer. 
The  cancellation  achieved  by  the  movement  of  the  disks  was  limited 
only  by  the  overall  system  noise  level  because  both  signals  enter 
through  the  sldelobes.  This  is  seen  in  Figure  10  as  the  series  of 
dots  at  the  bottom  of  the  plot  from  +36  to  +45  degrees.  The  nulling 
technique  works  well  while  the  main  beam  is  looking  up  because  it 
seeks  to  minimize  the  total  power  received  by  the  antenna,  and  only 
works  in  the  sidelobe  region.  It  may  be  called  a  manual  version  of 
the  gradient  search  algorithm  with  a  man  in  the  loop. 

The  situation  is  entirely  different  when  the  antenna  is  pointed 
into  the  ground  during  a  measurement.  The  direct  ray  becomes  a 
sidelobe  jammer,  because  it  enters  through  the  sidelobe  as  the  dish 
is  pointed  down.  Since  the  two  disks  are  sized  and  positioned  on 
the  reflector  surface  to  cancel  only  signals  entering  in  the  side- 
lobes,  the  direct  ray  is  the  only  ray  affected  by  the  nulling  tech¬ 
nique.  As  the  main  beam  scans  the  ground  during  the  vertical  move¬ 
ment,  the  energy  received  by  the  main  beam  is  not  affected  by  the 
nulling  action  of  the  disks.  In  fact,  the  cancellation  achieved  Is 
limited  by  the  amount  of  energy  entering  the  main  beam.  That  is  to 
say,  the  disks  are  adjusted  until  minimum  power  is  received  by  the 
dish.  At  this  point,  Che  antenna  is  receiving  equal  amounts  of 
of  direct  energy  through  a  null  in  the  sidelobe  and  reflected  energy 
by  way  of  the  main  beam. 
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The  highest  return  is  at  -36  degrees  and  corresponds  to  the 
grassy  area  in  front  of  the  antenna  building.  This  reflection  is 
-35  dB  relative  to  the  direct  ray  and  agrees  with  a  recent  measure¬ 
ment  by  Franchi  and  Poles^. 

The  second  highest  return  is  at  -24  degrees.  The  strength  is 
-41  dB  relative  to  the  direct  ray  and  is  the  reflection  from  the 
security  fence. 

The  third  highest  return  is  in  between  -40  and  -48  degrees. 

It  corresponds  to  the  flat  grassy  area  next  to  the  road  that  en¬ 
circles  the  building. 

The  following  argument  is  offered  for  the  scatter  in  the 
achieved  cancellation  on  the  sky  side,  close  to  the  main  beam  In 
Figure  10.  At  the  Ipswich  Facility,  the  town  has  located  a  large 
water  tower  across  the  street  from  the  transmitter  building.  It  is 
very  tall  and  extends  from  +9  to  +15  degrees  in  elevation.  It  is 
conceivable  that  the  cancellation  on  the  sky  side  was  limited  by 
reflections  off  this  tower  in  a  manner  identical  to  the  performance 
of  the  antenna  looking  into  the  ground.  Further  testing  is  needed 
to  verify  the  reflection  off  the  water  tower.  If  this  reflection 
is  at  -50dB  relative  to  the  direct  ray,  it  could  seriously  limit  the 
measurement  of  ultra-low  sidelobe  antennas  at  the  RADC  range. 

4.  CONCLUSION 

A  description  of  two  applications  of  a  novel  technique  for 
pattern  control  in  reflector  systems  has  been  described.  In  the 


first,  an  unjamined  coherent  reference  signal  has  been  provided  for 
the  proper  operation  of  the  Scientific/Atlanta  2020  Antenna  Measure¬ 
ment  System.  This  has  been  accoraplis’hed  without  the  need  for  an 
expensive  buried  cable  system,  and  works  better  (the  full  dynamic 
range  of  the  receiver  is  assured)  than  other  alternatives  that  were 
considered. 

In  the  second,  a  direct  measure  of  the  angular  variation  of  the 
ground  clutter  at  the  test  range  has  been  obtained.  For  the  first 
time,  the  locations  of  the  scattering  centers  that  contribute  to  the 
reflected  signals  at  the  measurement  site  were  identified. 
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A  NUMERICAL  METHOD  FOR  APPROXIMATING  ANTENNA  SURFACES 
DEFINED  BY  DISCRETE  SURFACE  POINTS 

Richard  Q.  Lee  and  Roberto  Acosta 
National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland.  Ohio  44135 

1.  INTRODUCTION 

Reflector  antennas  on  earth  orbiting  spacecrafts  generally 

cannot  be  described  analytical 1y .  The  reflector  surface  Is 

subjected  to  large  temperature  fluctuation  and  gradients,  and 

Is  thus  warped  from  Its  true  geometrical  shape.  Aside  from 

distortion  by  thermal  stresses,  reflector  surfaces  are  often 

purposely  shaped  to  minimize  phase  aberatlons  and  scanning 

losses.  To  analyze  distorted  reflector  antennas  defined  by 

discrete  surface  points,  a  numerical  technique  must  be  applied 

to  ccpute  an  Interpnlatory  surface  passing  through  a  grid  of 

discrete  points.  Although  numerical  techniques  for  analyzing 

reflector  surface  distortions  and  doubly  curved  reflector  sur- 
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faces  have  been  reported  In  open  literatures,  all  these 
techniques  are  rather  complicated  and  Involve  lengthy  computa¬ 
tions.  In  this  paper,  a  simple  numerical  technique  based  on 
Taylor's  expansion  and  finite  differences  approximation  Is 
presented.  By  applying  the  numerical  technique  to  approximate 
the  surface  normals  of  a  distorted  reflector  surface,  the  aper¬ 
ture  field  over  a  near-field  plane  Is  obtained  using  geometric 


optics  technique.  The  aperture  Integration  method  Is  then 

4 

applied  to  find  the  radiation  patterns. 

2.  FORMULATION 

The  proposed  numerical  technique  Is  a  quadratic  approxima 
tion  of  a  distorted  surface  defined  by  a  grid  of  nine  points. 
Let  (x,y,z)  be  the  coordinate  of  any  point,  P,  on  the  surface 
which  Is  described  by 


2  »  z(x,y) 


The  position  vector  of  P  Is  given  by 

A 

r  a  XX  ♦  yy  ♦  zz 

To  obtain  an  approximation  for  the  surface  z(x,y),  the  position 

vector  Is  replaced  by  Us  Taylor's  expansion  around  r^(x^, 

y  ,z  ),  the  closest  point  to  the  point  of  Incidence.  That  Is 
0  0 
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where  the  partial  derivatives  of  r  are  evaluated  at  "TpCxQ, 

terms  higher  than  the  quadratics  of  x  and  y 

3  3 

have  been  neglected  as  Indicated  by  0(x  ,y  ).  The  surface 
may  then  be  approximated  by 

z(«.»)  -  ?  •  -f  -  ^0  •  <«  -  "o’  S  *  <»  -  »o) 

2  2 
1  /  x2  a  z  ,,  ,  a  z 

'  ^  ^  ^ 

/  a^zl  3  3, 

♦  (y  -  ^0^  ^  r  o<*  -y  > 

All  the  partial  derivatives  of  z  are  evaluated  at 
^o’^o^'  normals  to  the  Interpolatory  surface  can  be 
found  by  taking  the  gradient  of  z(x,y)  as  follows: 


^  6 2  ^  dZ  ^  ^ 

"  “ i7*  *  ^ 

n  «f^»(x-x)^»(y-y)  f-f-  f  0(x^,y^) 
X  ax  '  0,2  0  axay  '  ' 

oX 

nz  = 


where  (n  ,n  ,n  )  are  the  x,  y,  and  z  components  of  n  with 
X  y  z 

the  normal  pointing  toward  the  source.  Finite  differences  can 
be  applied  to  approximate  the  partial  derivatives.  Since  the 


partial  derivatives  are  of  second  order,  a  grid  of  nine  points 
Is  required. 

The  partial  derivatives  may  be  expressed  as 
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where  (x^.y^.z^),  (x^.y^.Zj)  •••  <*8’^8^8^  coordi¬ 

nates  of  the  points  around  the  closest  point  (^q'^o'^o^’ 

3.  DISCUSSION  AND  RESULT 

In  applying  the  numerical  technique  to  compute  the  surface 
normals,  the  nine  closest  points  to  the  point  of  Incidence  are 
located  by  searching  an  array  of  discrete  surface  points  defin¬ 
ing  the  distorted  antenna  surface.  Ray-tracing  Is  then  used  to 


compute  the  reflected  field  from  the  approximated  normals.  The 
secondary  far  field  patterns  are  obtained  from  the  Fourier 
transformation  of  the  aperture  field  distribution. 

To  Illustrate  the  numerical  technique,  the  numerical  algo¬ 
rithm  has  been  programmed  and  applied  to  the  antenna  configura¬ 
tion  shown  In  figure  1.  The  surface  grid  points  for  an 
artificially  "enlarged"  parabolic  surface  are  generated.  Thus, 
the  boundary  points  along  the  rim  of  the  actual  reflector  sur¬ 
face  can  be  Included  In  the  computation.  The  computed  secondary 
far  field  patterns  for  the  symmetric  and  the  offset  cases  are 
shown  In  figures  2  and  3  respectively.  In  the  figures,  the 
E-plane  radiation  patterns  computed  from  an  analytic  expression 
of  a  parabolic  surface  are  drawn  with  solid  lines,  while  that 
computed  numerically  are  drawn  with  dots  and  broken  lines.  The 
E-plane  radiation  patterns  computed  numerically  from  a  grid  of 
50x50  surface  points  of  ~0.7  X  spacing  was  found  to  compare 
exactly  to  that  computed  from  an  analytic  expression.  For  the 
offset  case,  the  aperture  surface  Is  described  by  fewer  points; 
therefore,  a  higher  density  of  target  points  Is  required  to 
produce  the  same  accuracy.  With  a  scan  step  of  x/3  for  the 
symmetric  case  and  X/6  for  the  offset  case,  the  computer  time 
required  to  compute  the  secondary  radiation  patterns  numerically 
Is  25  sec  for  the  symmetric  case  and  17  sec  for  the  offset  case, 
respectively.  For  comparison  purpose,  the  radiation  pattern 


for  the  offset  case  has  been  computed  using  spline  function 
approximation.^  The  computed  E-plane  pattern  Is  shown  In 
figure  3. 
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Abstract 

To  verify  the  properties  of  transform  feed  antennas,  Rome  Air  Devel¬ 
opment  Center  has  built  and  tested  an  experimental  model  for  azimuth-only 
scan.  The  30-wavelength  line  source  array  Is  space  fed  by  a  cylindrical 
bootlace  lens  constructed  of  parallel  plate  waveguide  with  coaxial  mono¬ 
pole  elements.  The  transform  feed  Is  a  parallel  plate  waveguide  Rotman 
lens. 

This  antenna  Is  shown  to  have  22%  signal  bandwidth  at  5'°  scan 
when  time  delays  are  applied  at  the  Rotman  lens  inputs.  Low  sldelobe 
patterns  are  achieved  over  the  operating  band  and  scan  region.  Low 
sldelobe  difference  patterns  are  obtained  without  the  need  for  a  separate 
feed  network.  With  proper  design,  total  antenna  losses  could  be  limited 
to  less  than  3  dB  Including  element  mismatch,  spillover,  waveguide  and 


cable  losses. 


1.  INTRODUCTION 


Transform  feed  antennas  have  long  been  considered  a  solution  to 
the  Inherent  dispersion  of  phased  arrays  [1],  [2].  Also  known  as 
"completely  overlapped  subarray"  antennas,  they  combine  the  advantages 
of  wide  Instantaneous  bandwidth  and  pattern  control  [3].  The  Rome 
Air  Development  Center  (RADC)  model,  shown  In  Figs.  1-3,  Is  an  X-band, 
azimuth-only  scan  antenna,  designed  to  demonstrate  the  properties  of 
this  generic  architecture. 

The  key  feature  of  the  antenna  Is  Its  nearly  constant  gain  versus 
frequency  at  scan  angles  out  to  55°  from  broadside:  The  16  time  delay 
steered  subarrays  correct  for  the  natural  squint  of  the  phase-steered 
aperture  array.  In  addition,  amplitude  tapering  need  only  be  applied 
at  the  subarrays  to  produce  a  tapered  aperture  Illumination  for  low 
sldelobes.  Reasonably  good  difference  patterns  are  obtained  without 
the  need  for  a  second  feed  network.  Finally,  the  space-feed  structure 
can  result  In  a  very  efficient  antenna. 

This  paper  discusses  the  results  obtained  with  the  RADC  model: 
Signal  bandwidth  Is  greater  than  22%  at  the  55°  scan  angle:  Peak  and 
average  sldelobes  are  below  -20  dB  and  -30  dB,  respectively,  over  the 
8-10  GHz  band,  at  both  0°  and  55°  scan.  In  both  the  sum  and  difference 
patterns.  Losses  due  to  element  mismatch,  spillover,  waveguide  and 
cables  are  less  than  2.5  dB.  Recommendations  for  Improvements  are 
also  presented. 


2.  HARDWARE  DESCRIPTION 


Thero  are  three  major  components  in  the  completely  overlapped 
subarray  antenna  ( COSA) ,  shown  in  Fig.  3.  The  largest  is  the  cylin¬ 
drical  objective  lens,  whose  front  face  is  the  60-element  aperture 
array.  A  16-beam  Rotman  lens  feeds  that  objective  lens  through  a 
2-foot  radius  parallel-plate  waveguide  cavity.  The  Rotman  lens 
linear  output  face  is  centered  at  the  objective  lens  focal  point. 

The  amplitudes  applied  to  the  Rotman  lens  beam  ports  are  set  in  the 
third  major  component,  the  corporate  feed.  Amplitude  weighting  is 
set  witli  fixed,  coaxial  attenuators.  The  coaxial  line  stretchers 
shown  are  for  phase  trimming  only:  the  time  delays  are  in  the  coax 
cables  joining  them  to  the  beam  ports.  The  three  power  splitters 
shown  in  Fig.  3  are  commercial  strlpllne  units.  The  small  two-way 
splitter  was  later  replaced  with  a  "magic  tee"  hybrid  for  forming 
simultaneous  sum  and  difference  patterns. 

2.1  Objective  Lens 

The  purpose  of  the  cylindrical  lens  is  to  focus  an  incident 
wave  onto  the  Rotman  lens  feed.  Because  scanning  is  performed  by 
phase  shifters  in  the  transmission  lines  joining  the  two  faces,  it 
need  not  be  a  wide  angle  lens. 

A  cylindrical  lens  was  chosen  for  two  reasons  [4] :  the  line 
lengths  are  all  the  samej  and  it  allows  a  very  short  focal  length, 
resulting  in  a  compact  antenna.  Usually,  the  back  face  elements  are 


loc^iLi'd  directly  behind  their  corresponding  iront  face  elements 
[ 1  )  ['^'1.  which  was  our  initial  approach.  Unfortunately,  the  resulting 
element  spa.  Ings  along  the  cylindrical  face  varied  by  more  than  33%. 
Under  rhai  c  i  rcurac  tance ,  impedance  raatchinir  tlie  entire  back  face  with 
a  single  clement  design  (probe  U ngtii  and  ground  plane  spacing)  was 
not  possible. 

Tlte  final  design  of  the  lens,  as  Fig.  3  siiows,  uses  uniformly- 
spaced  back  face  elements.  As  described  in  [7],  tiiat  (carefully 
chosen)  spacing  does  increase  the  scanning  aberrations  of  the  lens, 
but  that  only  affects  the  antenna's  performance  at  the  extreme  limits 
of  scan  angle  and  frequency. 

To  test  the  cylindrical  lens  Independently  of  the  rest  of  the 
antenna,  we  connected  a  4-way  power  divider  to  the  center  4  focal 
array  elements.  With  the  outside  two  ports  weighted  -10  dB  relative 
to  the  center  two,  this  simulates  the  ideal  Rotman  lens  outputs  for 
broadside  scan.  Measured  E-plane  patterns  are  shown  in  Fig.  4. 


2.2  Rotman  Lens 

The  essential  structure  of  the  Rotman  lens  is  the  same  as  the 
objective  lens;  All  elements  are  coaxial  monopole  probes  with  their 
center  conductors  extending  Into  parallel-plate  guide.  Each  array 
of  probes  is  backed  by  an  aluminum  spacer  that  forms  an  image  plane 
(see  inset  of  Fig.  3).  In  contrast  to  the  cylindrical  lens,  the 
beamformer  in  the  COSA  has  to  be  a  wide-angle  lens.  Its  outermost 
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beams  Illuminate  the  edges  of  the  objective,  which  are  at  +  45° 
relative  to  the  focal  array's  broadside. 

Oar  design  used  g=1.137  [8],  with  perfect  foci  at  ^30°.  The 
coax  cables  are  phase-trimmed  to  the  proper  lengths  "w."  A  quadratic 
length  Is  subtracted  from  w  to  correct  for  the  shape  of  the  focal 
array:  Ideally  It  should  be  circular,  with  radius  equal  to  the  objec¬ 
tive's  focal  length. 

A  major  modification  to  the  Rotman  lens  was  the  Interleaving  of 
terminated  (dummy)  elements  with  the  16  active  beam  ports.  The  re¬ 
quired  angles  of  the  16  beams  dictated  a  spacing  of  about  •7Xq 
(fQ=9  GHz),  lomaslc  and  Hessel  [9]  have  shown  (Fig.  5)  that  the  element 
pattern  Is  severely  narrowed  for  spaclngs  greater  than  .6X.  The 
dummy  beam  ports  decreased  the  spacing  to  0.35Xo,  greatly  Improving 
the  Illumination  of  the  antenna  ports,  as  illustrated  in  Fig.  6. 

Unfortunately,  those  dummy  ports  absorb  half  the  power  incident 
on  the  beam  port  array,  with  a  consequent  3  dB  reduction  in  antenna 
gain.  On  the  other  hand  it  is  certainly  possible  to  sum  those  extra 
ports  in  a  second  corporate  feed.  Alternatively,  a  smaller  Rotman 
lens  might  reduce  the  length  of  the  beam  port  arc,  with  a  consequent 
reduction  In  element  spacing. 

Element  spacings  of  both  the  beam  port  and  the  antenna  port  arrays 
were  nonunlforra.  However,  unlike  the  cylindrical  array,  the  variation 
in  spacing  was  only  about  10%.  Impedance  matching  those  arrays  for  the 


average  spacing  has  apparently  worked  quite  well.  Patterns  of  a  few  of  / 

the  Rotman  lens  beams  are  Fig.  7.  /.'N'. 


2.3  Corporate  Feed 

The  corporate  feed,  or  variable  power  divider  Is  designed  to  apply 
the  proper  amplitudes  and  phases  to  the  Rotman  lens  inputs.  Fixed  co¬ 
axial  attenuators  set  the  amplitude  weights.  They  are  followed  by  co¬ 
axial  line  stretchers  for  minor  phase  adjustments.  These  were  required 
because  the  near-ln  sldelobe  levels  are  extremely  sensitive  to  errors 
in  the  corporate  feed. 

We  have  two  sets  of  cables  to  link  the  corporate  feed  outputs  to 
the  beam  ports:  one  for  broadside  scan;  the  second  for  35°  scan.  The 
first  set  is  phase  trimmed  to  correct  for  the  shape  of  the  beam  port 
arc,  whose  center  of  curvature  Is  not  at  the  center  of  the  antenna  port 
array.  These  cable  lengths  equalize  the  electrical  path  lengths  from 
the  power  divider  outputs  to  the  center  of  the  antenna  port  arc.  The 
second  set  Includes  this  correction,  but  also  implements  the  time  delay 
steering. 

The  pattern  of  the  Rotman  lens  feed,  when  driven  froo  the  16-way 
power  divider  is  Fig.  8.  The  dashed  line  of  this  figure  Is  the  ideal 
pattern.  With  the  Rotman  lens  feed  inserted  into  the  wavigulde  cavity 
behind  the  cylindrical  lens,  this  will  be  the  shape  of  the  aperture  dis¬ 
tribution.  The  variation  between  the  desired  and  actual  distribution  is 
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due  more  to  amplitude  errors  In  the  corporate  feed  than  to  any  fault  in 
the  Rotman  lens. 

3.  EXPERIMENTAL  RESULTS 

3.1  Broadside  Scan  Measurements 

Measured  sum  and  difference  patterns  are  shown  in  Figs.  9  and  10. 

Random  errors  are  -joraewhat  higher  at  8.25  GHz  than  at  other  frequencies, 
accounting  for  the  higher  average  sldelobe  level  as  well  as  filling  of 
the  central  null  in  the  difference  pattern. 

Fig.  11  is  an  analysis  of  the  antenna  gain.  The  measured  gain  is  the 
lowest  curve,  with  encircled  points.  The  other  of  the  two  lower  curves  is 
the  predicted  gain,  found  as  follows:  The  theoretical  maximum  gain  is 
that  of  the  flared  aperture,  and  is  given  by  Schelkunoff  &  Frils  [10]  and 
is  the  top  curve  on  the  plot.  The  two  crosshatched  areas  represent 
unavoidable  losses.  The  first  is  the  sum  of  mismatch  losses  of  all  five 
probe  arrays,  calculated  from  Toraaslc's  theory  [9].  The  second  is  the 
measured  loss  in  the  objective  lens  and  Rotman  lens  cables,  plus  the  cal¬ 
culated  loss  through  the  2  "  space-feed  region.  The  next,  unshaded  region 
is  the  3  dB  loit  to  the  dummy  Rotman  lens  beam  ports.  The  most  significant 
loss  in  the  system,  however,  is  the  6-7  dB  in  the  corporate  feed,  because 
we  have  used  fixed  attenuators  for  amplitude  weighting.  Use  of  directional 
couplers  instead  would  reduce  that  loss  substantially. 

A. 2  55°  Scan  Measurements 

For  55°  scan,  the  theoretical  maximum  gain  Is  reduced  by  2.4  dB  due 
to  cos  55°  scanning  loss.  In  addition,  the  aperture  array  has  a  much  higher 


active  impedance.  Except  near  10  GHz,  where  the  grating  lobe  begins  to 
appear  on  the  opposite  side  of  the  pattern,  the  predicted  and  actual  gains 
agree  to  with  1.5  dB,  indicating  that  all  significant  losses  are  accounted 
for.  The  reduction  in  antenna  gain  is  apparent  also  in  the  higher  average 
sldelobe  levels  in  the  patterns  of  Fig.  13. 

The  9  GHz  pattern  has  a  prominent  sldelobe  at  27°.  It  is  caused  by 
multiple  reflections  between  faces  of  the  objective  lens.  The  impedance 
mismatch  of  the  aperture  array,  on  transmit,  reflects  some  power  back 
toward  the  feed.  The  cylindrical  array  then  re-reflects  a  portion  of  that 
back  toward  the  aperture,  to  be  radiated  as  an  error  lobe.  If  was  the 
inter-element  phase  required  to  scan  the  array  to  55°,  the  angle  of  the 
error  lobe  will  be  6 g=sin”^ [ (3ip-2Tr )/kd] .  At  8.25  GHz,  the  array 
Is  matched  better  at  55°,  and  the  mismatch  error  lobe  Is  indistinguishable 
from  the  remaining  sidelobes.  At  9.75  GHz,  the  lobe  appears  at  40°. 

In  the  8.25  and  9.75  GHz  patterns,  there  are  pairs  of  sidelobes  at, 
respectively,  -5°,  -28°  and  +3°,  -18°.  These  are  caused  by  the  line- 
length  phase  shifts  at  the  aperture.  Those  cables  are  trimmed  modulo-Zn 
at  9  GHz.  At  all  other  frequencies,  the  inter-element  phase  shift  will  have 
periodic  discontinuities,  causing  error  sidelobes.  This  is  a  phenomenon 
associated  with  line-length  phase  shifters,  rather  than  any  peculiarity  of 
the  COSA  architecture. 

The  most  important  feature  of  these  patterns  is  that  the  main  beam  is 
at  exactly  55°  for  all  frequencies.  The  time  delay  steering  has  eliminated 
the  squint  that  would  have  shifted  the  beam  to  63.3°  and  49.1°  at  8.25 
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Figure  3.  Hardware  Outline  Drawing 


40  20  0  M 


ao  40  to 


Angl«  (deg) 
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Figure  5.  Array  Geonetry  and  Element  Patterns 
(after  Toinasic  &  Hessel) 
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HIGH  EFFICIENCY  NULTIBEAN  ANTENNA  DESIGN 


P.  INGERSON,  C.  A.  CHEN 
TRW  Antenna  Systems  Laboratory 
Redondo  Beach,  CA 

I .  INTRODOCTION 

In  many  millimeter  wavelength  satellite  communication  systems,  there 
is  a  need  to  provide  multiple  high  gain  independent  beams  covering  an 
angular  region  of  space.  It  is  usually  highly  desirable  to  provide  the 
multiple  beams  from  a  single  aperture  if  possible  to  reduce  size,  weight 
and  complexity  of  the  antenna  system.  In  addition  to  providing  maximum 
gains,  in  many  cases  it  is  also  required  to  have  low  sidelobe  radiation 
patterns  to  minimize  out  of  beam  interference. 

The  desire  to  provide  multiple  high-gain  beams  from  a  single  aperture 
(lens  or  reflector)  usually  leads  to  beam  peak  gains  which  are  typically 
-3dB  or  more  lower  than  one  can  achieve  using  a  single  optimum  feed  at  the 
lens  or  reflector  focus.  As  will  be  discussed,  this  is  due  to  a 
fundamental  limitation  of  conventional  focusing  lens  or  reflector  type 
antennas.  The  problem  is  to  provide  the  multiple  beams  from  the  antenna 
with  gains  comparable  to  the  optimized  single  beam  antenna. 

The  basic  limitations  of  focusing  antenna  systems  is  that  the 
secondary  radiation  pattern  is  directly  related  to  the  aperture  size.  In 
the  present  approach,  this  limitation  is  overcome  by  using  non-focusing 
(diverging)  reflector  or  lens  antennas.  These  antennas  can  be  designed  to 
provide  the  required  secondary  beamwidth  essentially  independent  from  the 
aperture  size. 
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Given  a  minimum  required  feed  horn  size  or  feed  separationr  a  single 
non-focusing  aperture  can  be  designed  using  the  method  presented  which 
will  allow  almost  any  beam  to  beam  separation  or  cross-over  level  desired 
while  providing  beam  gains  which  are  essentially  equal  to  the  optimum 
single  beam  of  a  focused  aperture. 

The  independence  from  feed  to  feed  separations  of  this  design 
technique  is  particularly  useful  in  the  cases  where  switches#  receivers  or' 
amplifiers  are  to  be  placed  as  close  to  the  feed  elements  as  possible# 
since  these  components  often  determine  the  minimum  allowable  separations. 

This  paper  describes  the  design  of  a  37  beam  EHF  (43.5  -  45.5  GHz) 
antenna  providing  high  ga^n,  contiguous  beam  coverage  of  the  earth  from  a 
synchronous  orbit  satellite.  The  angular  region  of  the  earth  from 
synchronous  orbit  is  approximately  ±  8.68  degrees#  requiring  about  a  2.9 
degree  beam  separation  and  a  3.3  degree  HPBW  to  provide  contiguous 
coverage  at  the  -3dB  level  or  better  as  shown  in  Figure  1. 

II.  THEORY 

A  fundamental  limitation  arises  when  attempting  to  produce  contiguous 
spot  beams  space  approximately  one  half-power  beamwidth  apart  from  a 
single  conventional  reflector  (or  lens  aperture) .  There  is  a  fundamental 
disparity  between  the  feed  horn  aperture  size  required  to  maximize  the 
gain  of  the  antenna  and  the  horn  size  necessary  to  permit  packing  beams 
approximately  one  beamwidth  apart.  It  can  be  shown  that  the  feed  horn 
diameter  required  to  maximize  gain  is  substantially  larger  then  the  horn 
diameter  that  will  permit  close  (-3dB  cross-over)  beam  packing  in  a  single 
conventional  focusing  reflector  or  lens.  Accordingly#  horn  aperture  size  * 
that  yield  maximum  gain  lead  to  beam  separations  much 
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larger  than  one  half-power  beanwidth  even  though  the  feed  horn  apertures 
are  contiguous.  This  fundamental  limitation  of  a  conventional  reflector 
and  lens  is  independent  of  the  focal  length  to  aperture  diameter  (F/D) 
ratio.  This  limitation  is  demonstrated  in  Figure  2  which  shows  the 
relative  gain  of  a  conventional  reflector  as  a  function  of  the  reflector 
included  angle  [  ^  ire).  Keeping  the  feed  pattern  constant,  corresponding 
approximately  to  a  1  inch  4  wavelength  (  x))  aperture  feed  horn,  d^nc 
was  varied  from  the  'alue  that  yielded  maxiamm  gain  for  a  3.3  degree  beam 
(  ^inc  2i  32  degrees)  to  other  values  of  einc  while  maintaining  the 
beamwidth  constant  at  3.3  degrees.  Figure  2  shows  that  e^^c  Bust  decrease 
to  approximately  15  degrees  before  -3dB  beam  cross-over  is  achieved.  In 
other  words,  if  one  uses  37  one  inch  diameter  feed  horns  in  a  single 
conventional  reflector  to  achieve  earth  coverage  with  37  3.3  degree  spot 
beams,  the  achievable  aperture  efficiency  of  each  beam  is  approximately 
3dB  less  than  that  obtainable  from  an  optimized  reflector  design  using  a 
single  1.0  inch  (4^)  diameter  feed  horn. 

II. b 

There  are  three  basic  approaches  that  can  be  applied  to  the 
fundamental  problem  associated  with  close  beam  spacing  (-3dB  beam 
cross-over)  to  increase  the  minimum  coverage  gain,  these  are: 
o  Increase  the  number  of  beams  used  to  provide  coverage 

o  Ose  more  than  one  reflector  aperture  to  produce  the  collection  of 

spot  beams 

o  Use  non-focusing  apertures  which  circumvent  the  fundamental 
limitation  problem. 


If  we  wish  to  fix  the  number  of  beams  from  system  compelxity 
considerations,  then  only  the  last  two  approaches  are  applicable 


As  pointed  out,  if  the  feed  born  diemnsions  are  selected  to  maximize 
the  reflector  gain,  the  beams  generated  by  two  adjacent  contiguous  feeds 
will  be  spaced  much  further  apart  than  one  half-*power  beamwidth.  One 
method  of  obtaining  a  pattern  of  contiguous  spot  beams  crossing  over  at 
the  -3dB  level  is  to  generate  adjacent  beams  from  separate  reflectors. 

This  technique  is  illustrated  in  Figure  3  which  shows  the  feed  born  sets 
assigned  to  three  separate  reflectors.  Reflector  Nc .  1  radiated  13  beams 
while  reflectors  No.  2  and  No.  3  radiate  12  beams  e< ch.  The  beam 
positions  produced  by  each  feed  set  are  also  shown  .n  the  Figure.  The 
closest  spacing  that  occurs  between  beams  produced  by  adjacent  feed  horns 
in  the  same  set  is  about  1.73  z  half-power  beamwidth.  This  larger 
separation  between  beams  produced  by  adjacent  feeds  allows  the  included 
angle  of  the  reflector  to  be  larger  and  much  nearer  to  the  optimum  feed 
horn  dimension. 

The  reflector  diameter  required  for  this  three-reflector  system  is 
about  5  inches  (20X )  for  each  reflector.  Because  the  feed  horn  dimension 
is  now  reasonably  close  to  that  of  an  optimum  single  feed  producing  a  3.3 
degree  beam,  the  gain  loss  from  excessive  illuminat:on  spillover  is 
reduced  to  about  0.6dB  as  shown  in  Figure  2.  The  edge  illumination  is, 
however,  still  high  and  hence  there  is  no  means  of  controlling  the 
sidelobe  levels. 

Sidelobe  level  control  with  high  edge  illumination  is  possible  only 
using  composite  multiple  feed  clusters  or  using  artificial  edge  tapering 
of  the  aperture  using  lossy  material.  The  composite  feec  technique  would 
require  additional  feed  elements  and  a  complex  ..nterconnecting  feed 
network.  The  artificial  edge  tapering  would  substantially  reduce  the 
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The  fundamental  limitation  of  focusing  aperture  can  be  circumvented  by 
utilizing  a  non-focusing  hyperbolic  reflector  (or  diverging  lens) .  The 
differences  between  focusing  and  non-focusing  apertures  are  depicted  in 
Figure  4.  The  shaped  reflector  result  can  also  be  approximated  quite 
closely  by  using  a  parabolic  surface  contour  combined  with  a  highly 
tapered  aperture  illumination  and  axial  defocusing  of  the  feed  horn  as 
shown  in  Figure  5.  With  axial  defocusing  the  secondary  beamwidth  can  be 
made  essentially  independent  of  the  reflector  diameter.  With  a  hyperbolic 
reflector,  any  feed  at  the  focus  point  located  at  the  right  side  of 
reflector  can  be  thought  of  as  producing  an  image  feed  at  its  image  focus 
point  whose  pattern  HPBW  is  reduced  by  the  reciprocal  of  the  magnification 
factor  M.  The  beamwidth  of  secondary  pattern  is  then  controlled  by  the 
magnification  factor  M  and  horn  aperture  and  is  essentially  independent  of 
the  diameter  of  reflector  unlike  the  parabolic  reflector.  The 
relationship  between  K,  horn  aperture  and  secondary  BPBW  is  shown  in 
Figure  6. 

Therefore,  by  using  the  hyperbolic  reflector,  we  can  optimize  the  gain 
and/or  sidelobe  performance  first  and  then  adjust  the  diameter  of  the 
reflector  to  permit  close  (-3dB  cross-over)  beam  packing.  All  37  beams 
with  -3dB  cross-over,  as  shown  in  Figure  1  can  be  produced  from  a  single 
20.7  inch  diameter  reflector.  The  extremely  low  edge  illumination  (-30dB) 
utilized  provides  negligible  feed  illumination  spillover,  higher  gain  and 
beam  efficiency,  and  very  low  sidelobes  for  the  beam. 

The  advantages  of  the  design  technique  aret 

Allows  high-gain,  high  efficiency  multiple  beams  from  a  single 
antenna  aperture  (lens  or  reflector) 


o 


o 


Independent  of  feed  size  and  feed  to  feed  separation  over  a  very 
wide  range 

o  Allows  very  wide  range  of  beam  to  beam  separations  or  cross-over 
levels 

o  Allows  efficient  and  simple  sidelobe  control  by  illumination 
taper 

III.  DETAIL  ANTENNA  DESIGN 

The  multi-beam  antenna  is  shown  in  Figure  7  and  the  antenna  system 
block  diagram  in  Figure  8.  The  antenna  system  comprises  a  single 
reflector  and  37  feed  modules.  The  single  reflector  is  an  offset 
hyperboloid  with  21.66  inch  projected  diameter,  a  5.75  magnification  and  a 
19.34  inch  focal  length.  The  feed  modules  are  arranged  in  seven  rows 
with  trianglular  planar  grid  patterns.  As  shown  in  Figure  9,  this 
provides  37  independent  beams  of  approximately  3.3  degrees  HPBW  covering  a 
10  degree  elevation  angle  POV  on  the  earth.  Bach  feed  module  comprises  a 
feed  horn,  polarizer  and,  through  RF  switches,  a  downconverter/amplif ier. 
The  feed  horn  is  a  1-inch  square  aperture  multimode  horn.  In  use,  the  IF 
outputs  of  the  downconverters  are  connected  to  an  IF  switch  which  selects 
a  number  of  possible  beams  for  further  signal  processing. 

To  meet  the  reliability  goals,  the  37-horn  feed  system  utilizes  44 
downconverter  amplifier  modules.  In  the  baseline  redundancy  scheme,  there 
is  one  redundant  downconverter  for  each  row  and  a  failure  in  any  mdoule  in 
the  row  can  be  switched  out  as  shown  in  Figure  10.  With  the  redundancy 
switches,  the  probability  of  surviving  over  10  years,  with  3dB  or  less 
degradation  for  all  37  beams  is  estimated  to  be  98  percent.  The  baseline 
redundancy  scheme  requires  81  ferrite  switches.  The 


dovnconverter/amplif ier  module  is  comprised  of  a  fineline  mixer,  bandpass 
filter,  and  PET  amplifiers. 

The  downconverter/amplif ier  module  is  1.0  x  1.0  x  2.5  inches  in  size 
and  as  shown  in  Figure  11,  with  feed  horn  and  switch  included,  weighs 
about  5  ounces.  An  IP  frequency  of  5.4  i  1.0  GHz  was  selected  so  that  the 
signal  can  be  readily  dehopped  after  the  IP  switch.  Power  from  a  local 
oscillator  (LO)  is  distributed  to  the  44  modules  by  a  waveguide  power 
divider.  Each  module  requires  about  1.1  watts  of  dc  power  and  10  mW  of  LO 
power  at  39.1  GHz.  The  LO  drive  source  and  IP  switch  are  not  considered 
part  of  the  antenna  study. 

Key  objectives  of  the  antenna  design  are  lightweight,  minimum  power 
consumption  and  high  reliability.  Consequently,  the  downconverter 
amplifier  modules  incorporate  high  efficiency  diodes,  require  minimum  LO 
drive,  and  are  configured  for  a  soft  failure  mode,  using  one  redundant 
downconverter  per  row. 

The  breadboard  unit  of  downconverter  is  shown  in  Figure  11.  A 
performance  summary  is  given  in  Table  1.  For  the  performance  estimates  of 
the  antenna,  the  measured  pattern  data  of  this  breadboard  feed  was  used  as 
input  to  a  computer  program  which  is  capable  of  calculating  the  secondary 
patterns  of  concaved  hyperbolic  surfaces  with  scanned  feeds.  The  final 
calculated  pattern  coverage  is  shown  in  Figure  12.  The  performance 
estimates  are  given  in  Table  2.  In  the  following  section,  the  extremely 
close  agreement  between  the  measured  and  calculated  results  show  that  the 
computer  generated  designs  are  realizable. 

IV.  EXPERIMENTAL  VERIFICATION  OF  BREADBOARD  DESIGN 

In  order  to  verify  the  predicted  performance  characteristics  of  the 
antenna,  a  breadboard  model  was  built.  The  breadboard  model  is  comprised 


of  an  offset  hyperbolic  reflector,  and  a  subsection  of  the  37  element 
array  of  12  feed  horns  and  polarizers. 


The  experimental  measurement  of  the  breadboard  model  was  divided  into 
two  parts.  The  first  part  was  the  RP  far  field  measurement  of  the  offset 
hyperbolic  reflector  with  a  single  feed  horn.  This  single  feed  horn  was 
moved  over  the  focal  plane  to  find  the  best  feed  locations.  Thirty-seven 
feed  locations  were  then  identified  to  obtain  optimum  coverage  contour 
patterns.  The  measured  coverage  contour  patterns  were  verified  by  the 
theoretical  computer  calculation  inputting  the  measured  37  feed  locations 
into  computer  software.  It  was  understood  that  performance  of  the 
reflector  with  a  single  feed  horn  may  be  different  from  that  with  the  horn 
within  cluster  due  to  mutual  coupling  within  cluster.  The  blockage 
effects  of  the  array  were  simulated  by  using  a  metal  plate  of  the  full 
array  size  about  the  feed  horn. 

The  second  part  of  the  breadboard  measurement  was  the  RP  far  field 
measurement  of  the  offset  hyperbolic  reflector  with  12  feed  horn  cluster. 
In  this  measurement,  the  effects  of  mutual  coupling  within  cluster  were 
shown  to  be  negligible. 

I V .  a  a£Jlflflfluj:fiaflnl.fl£.Qffflet,Hmei:bQllc^BfiflfictQr,Mlth-SinalB..Eead..HQ£n 

The  antenna  configuration  for  these  measurements  is  shown  in  Figure 
13.  The  dimension  of  the  offset  hyperbolic  reflector  is  21.66  inches  with 
a  focal  length  of  19.34  inches.  The  feed  horn  was  the  1-inch  square 
aperture  multi-mode  horn  which  has  excellent  pattern  characteristics  over 
a  wide  beamwidth.  The  polarizer  was  a  quarter  wave  plate  vane  made  of 
Duriod.  The  axial  ratio  of  this  horn  with  polarizer  over  frequency  band 
(43.5  to  45.5  GHz)  was  less  than  0.3dB.  The  VSWR  of  the  feed  horn 
polarizer  was  less  than  1.04  VSWR. 


To  demonstrate  the  coverage  of  the  antenna,  only  20  beams  were  needed 


due  to  symmetry  in  antenna  configuration.  In  measuring  these  20  beams,  it 
was  necessary  to  adjust  the  feed  horn  location  to  obtain  the  optimum 
coverage,  up  to  10  degrees  elevation  angle  field  of  view  on  the  earth. 
Figure  14  shows  the  typical  measured  contour  patterns.  The  boresight 
beam,  beam  number  one,  is  shown  in  Fgiure  14a.  The  sidelobe  is  well  under 
30dB  with  almost  perfect  circular  contours.  The  beamwidth  is  3.3 
degrees.  When  the  beam  is  scanned  to  position  number  4,  as  shown  in 
Figure  14b,  it  can  see  that  beam  contour  shape  is  distorted  slightly, 
although  sidelobe  is  still  well  below  -30dB.  When  the  beams  were  scanned 
to  the  upper  elevation  direction  (Figure  14c)  and  lower  elevation 
direction  (Figure  14d) ,  the  distortion  of  beam  shape  was  different  because 
the  non-symmetr ical  offset  hyperbolic  reflector.  The  sidleobes  were  still 
well  under  -30dB.  The  mutual  coupling  within  cluster  in  these 
measurements  was  not  included  beacuse  of  single  feed  horn  measurement. 

The  20  measured  beam  contours  were  then  combined  to  show  the  contour 
coverage  of  the  offset  hyperbolic  reflector.  Figure  15a  shows  the 
composite  contour  coverage  at  44.5  GHz.  These  figures  show  that  the 
offset  hyperbolic  reflector  can  have  the  coverage  up  to  10  degrees 
elevation  angle  view. 

The  20  experimental  optimum  feed  locations  found  experimentally  were 
then  fed  into  a  computer  program  to  calculate  the  theoretical  contour 
coverage.  These  theoretical  contour  coverages  were  then  compared  to  the 
measured  data  as  shown  in  Figure  15.  This  verified  that  the  computer 
program  developed  for  offset  hyperbolic  reflectors  is  very  accurate. 
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The  purpose  of  this  veasurenent  was  to  find  both  the  effects  of  mutual 
coupling  within  the  cluster  and  the  blockage  due  to  cluster. 

A  feed  cluster  with  12  feed  horns  as  shorn  in  Figure  16  was 
assembled.  This  feed  cluster  represents  one-third  of  the  37  feed  horns. 
From  the  measured  patterns r  it  is  shown  that  the  effects  of  mutual 
coupling  and  blockage  of  cluster  are  negligible.  The  24  measured  contours 
are  combined  as  shown  in  Figure  17  to  show  the  composite  contour  coverage 
of  the  offset  hyperbolic  reflector  with  the  feed  cluster.  Figure  17b 
shows  the  theoretical  calculated  contour  coverage.  A  comparison  of 
Figures  17a  and  17b  again  shows  the  contours  are  almost  identical. 
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Table  1.  Downconverter  Performance  Summary 


RELATIVE  DEAM  POSITIONS  FOR  CONTIGUOUS 
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FIGURE  4.  EXAMPLES  OF  FOCUSING  AND  NON-FOCUSING  ANTENNAS 


FRONT  VIEW 


MAGNIFICATION 


SECONDARY  HALF-POWER  BEAMWIDTH  OF  A  HYPERBOLA 
AS  A  FUNCTION  OF  HORN  APERTURE  AND  MAGNIFICATION. 


FIGURE  7.  CONFIGURATION  OF  OFFSET 
HYPERBOLIC  REFLECTOR 


ANTENNA  SYSTEM  BLOCK  DIAGRAM 


FIGURE  10.  REDUNDANCY  SCHEME  USING  FERRITE 
CIRCULATOR  SNITCHES 
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(c)  Beam  No.  12 


(d)  Beam  No.  20 


FIGURE  14.  TYPICAL  MEASURED  CONTOUR 
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1  IXTRODLC'nOX 


In  anienna  engineering  ii  is  often  necessary  lo  compute  certain  parameters  which  are  used  to 
assess  the  input  and  radiative  properties  of  antennas.  Among  these  parameters  are  input 
impedance,  directive  gain,  power  gain,  and  efficiency,  all  of  which  are  derived  for  a  specified 
antenna  geometry  from  a  known  feed  network  and  the  corresponding  induced  current  over  the 
anienna  surface.  The  equations  which  govern  the  behavior  of  antennas  are  continuous 
mathematical  operations  which  relate  continuous  currents  on  the  anienna  surface  lo  a  specified 
source  field  or  feed.  Since  these  equations  cannot  be  .solved  analytically,  they  must  either  be 
discretized  and  solved  numerically  using  a  computer,  or  solved  approximately  using  variational 
techniques.  In  either  case,  the  solution  obtained  is  an  approximation  to  the  original  problem,  and 
the  figures  of  merit  mentioned  above  will  al.so  deviate  f  rom  their  respective  values.  The  purpose  of 
this  paper  is  to  describe  w'ays  lo  minimize  the  difference  between  the  computed  and  actual  values 
of  the  aforementioned  figures  of  merit  for  thin-wire  antennas.  It  is  shown  herein  that  a  code  such 
as  the  Numerical  Electromagnetics  Code  (.NEC)  of  Lawrence  Livermore  National  Laboratory  [l]  is  a 
powerful  tool  in  the  analysis  of  thin-wire  antennas  provided  that  the  model  is  accurate  and  the 
results  are  interpreted  properly. 

The  directive  and  power  gains  are  directly  related  to  the  radiated  or  far  field  of  the  antenna 
and  thus  exhibit  little  variation  from  their  true  values  when  small  computational  errors  are 
superimposed  on  the  current  [2].  This  is  not  always  the  case  with  the  input  impedance  and 
radiative  efficiency  of  antennas,  as  they  depend  highly  on  the  terminal  properties  and  near  fields  of 
the  antenna  as  well.  The  subject  of  this  paper  is  an  assessment  of  the  computational  accuracy  of 
the  input  impedance  and  radiative  efficiency  of  thin-wire  antennas  using  the  Numerical 
Electromagnetics  Code. 
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2  NL  MFRICAI  SOI  I  HONS  OF  THIN  WIRI  ANTFNNAS 

When  solving  thin  wire  antenna  problem  .  or  a  digital  computer,  a  code  such  as  \1  C  is  used  to 
solve  a  discretized  form  ot  the  electric  field  integral  equation  which  relates  an  applied,  or  incident, 
a.xially  directeti  electric  field  to  the  axially  directed  induced,  or  scattered,  electric  field  over  the 
surface  of  the  antenna.  A  problem  of  frequent  interest  is  that  of  a  vertical  monopole  or  whip 
above  a  perfect  ground  as  shown  in  Figure  1.  i'he  surface  of  the  antenna  is  frequently  modeled  as 
a  perfect  electric  conductor  except  at  di.screte  points  where  lumped  loads  are  to  be  placed.  It  is 
assumed  that  the  antenna  to  be  modeled  is  sufficiently  thin  so  that  the  fields  produced  by  all 
ncnaxialK  directed  currents  are  several  orders  of  magnitude  smaller  than  those  produced  by 
axiallc  directed  currents,  and  hence  non-axial  currents  can  be  neglected.  Further,  it  is  assumed 
mat  the  antenna  is  sufficiently  thin  that  the  induced  current  on  the  antenna  can  be  modeled  as  a 
filamentarc  current  along  the  axis  of  the  antenna,  even  though  the  boundary  condition  that  the 
tangential  electric  field  must  vanish  is  still  enforced  on  the  surface  of  the  antenna.  The  latter  is  an 
application  of  the  technique  of  extended  boundary  condition  which  is  used  to  formulate  scattering 
problems  for  closed  surfaces  [3].  When  applied  to  thin  wires,  this  technique,  together  with  the 
foregoing  assumption  that  the  radial  currents  at  wire  ends  can  be  neglected,  is  known  as  the  thin- 
wire  approximation  to  .scattering  problems. 

With  no  loss  of  generality,  this  discussion  is  confined  to  an  antenna  in  which  all  currents  flow 
in  the  /-direction.  Similar  expressions  to  those  which  follow  could  also  be  derived  if  the  structure 
to  be  modeled  is  more  complex  and  contains  skew  wires  as  well. 

The  thin-wire  formulation  for  the  whip  antenna  of  Figure  1  yields  the  following  integral 
equation: 


/;.  = — —  f  nz  K;(z.z-)dz-  +Hz)z{z) 

yw€„  dz^  i 


where 


I(z')»  the  unknown  z-directed  current  along  the  axis  of  the  antenna 
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Figure  1.  Vertical  thin-wire  monopole  above  infinite  perfect  ground  plane. 


(i(z.z')  =  — - — =^;r=::r^=:r  the  ihin-wire  approximalion  to  the  free-space 
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(ireen  s  f unction  along  the  wire  axis,  often  referred  to  as  the  "thin-wire  kernel’ 
and  /.(z)  =  the  distribution  of  impedance  loading  along  the  antenna,  which  for 

.V 

lumped  loads  reduces  to  /  (r  )  =  ^  7,  8(7  —7,  ) 

and  k  =  =  free  space  \xave  number. 

The  left-hand  side  of  ( 1 J  is  the  source  term  which,  for  most  thin-wire  antennas,  is  modeled  as 
a  constant,  highly  localized,  axially  directed  electric  field  incident  on  a  small  region  of  the  surface 
and  zero  elsewhere.  This  is  a  field-theory  model  of  a  voltage  source  which  is  connected  to  the 
antenna  across  two  closely  spaced  terminals.  In  NEC.  two  voltage  source  models  are  available. 
Thev  are  the  "current-slope-disconiinuity"  (CSl)J  source  and  the  "applied-electric-field"  (AEF) 
source  [l]. 

The  CSD  source  model  consists  of  a  highly  localized  voltage  source  across  a  biconical  junction 
of  two  wires  as  show  n  in  f  igure  2  The  axial  electric  field  generated  by  this  model  is  approximated 
by  that  of  a  thin  biconical  transmission  line  with  a  point  voltage  source  across  the  junction  of  the 
two  cones  [2].  This  approximation  al.so  requires  the  charge,  or  "current-slope."  to  be  discontinuous 
across  the  junction.  The  field  produced  by  this  model  becomes  more  localized  as  the  length-to- 
radius  ratio  (h/a)  of  the  wires  is  increased  and  works  best  for  wires  with  h/a  >  1000.  This  is 
especially  important  in  the  computation  of  input  impedance  and  is  discussed  in  detail  subsequently. 
I’lots  ol  the  axial  electric  field  and  current  over  the  feed  region  of  the  monopole  of  Figure  1  were 
generated  using  \F('  and  are  shown  for  h/a  =  100  and  1000  in  Figure  3.  The  height  of  the  antenna 
is  one-tenth  wavelength  and  the  source  is  modeled  using  the  CSD  option. 


The  .AFf  .source  is  simply  a  constant,  sinusoidally  varying  electric  field  incident  on  a  small, 
specified  section  ot  the  antenna  and  zero  elsew'here  as  shown  in  Figure  4.  Physically,  this  source 


Figure  2.  Vertical  ihin-wire  monopole  with  Currenl-Slope-Disconlinuity  (CSD)  source  model. 


represents  a  constant  voltage  across  a  \\  ire  gap  at  the  antenna  feed.  This  gap  is  usually  electrically 
small;  thus  care  must  he  taken  to  ensure  that  the  section  over  which  the  AHF  source  is  applied  is 
very  much  less  than  a  wavelength  and  significantly  less  than  the  length  of  the  wire  itself.  This 
will  result  in  a  computed  current  which  varies  smoothly  and  is  nearly  constant  over  the  feed 
region.  This  is  shown  to  facilitate  accurate  computation  of  the  input  impedance,  even  for  wires 
with  h/a  <  1000.  Plots  of  the  axial  electric  field  and  current  in  the  vicinity  of  the  AEF  source  for 
the  vertical  monopole  described  above  are  presented  in  Figure  5  for  h  a  =  100  and  1000.  It  is 
interesting  to  note  that  the  axial  field  generated  in  response  to  the  .AFl'  source  looks  more  like  the 
desired  pul.se  as  h/a  increases.  The  improvement  in  appearance  of  the  axial  near  fields  in  response 
to  both  sources  as  h  a  increases  is  a  direct  consequence  of  the  increasing  validity  of  the  thin-wire 
approximation  described  previously  as  h  a  becomes  larger,  thus  making  the  wire  appear  thinner. 

The  right-hand  side  of  ( 1 )  consists  of  two  terms;  the  first  of  these  represents  the  scattered 
field  on  the  surface  of  the  antenna  due  to  the  unknown  induced  current  l(z),  while  the  second 
represents  the  electric  field  across  loads  distributed  along  the  axis  of  the  antenna.  At  this  point  it 
should  be  emphasized  that  I'quation  4  is  only  accurate  for  wires  w  ith  radii  that  are  several  orders 
of  magnitude  less  than  the  operating  wavelength  and  for  wire  segments  with  lenglh-to-radius 
ratios  in  excess  ol  S  when  using  the  NFC  code  [l].  If  the  latter  condition  is  not  satisfied,  a  multi- 
term  expansion  for  the  free-space  Creen's  function  can  be  used  to  give  more  accurate  results  for 
short,  thin  wires.  One  such  expansion,  known  as  the  "extended  thin-wire  kernel."  is  available  in 
NFC  and  is  discussed  in  detail  in  [l]. 

liquation  1  cannot  be  .solved  analytically,  and  thus  it  must  be  discretized  and  solved 
numerically  with  the  aid  of  a  digital  computer.  The  most  commonly  u.sed  computational  technique 
to  discretize  and  numerically  .solve  electromagnetic  scattering  problems  is  known  as  the  method  of 
moments.  The  moment  method  is  used  in  Nl'C  and  is  discussed  briefly  herein.  For  further 
information  on  the  method  of  moments  the  reader  is  referred  to  [l]  and  [4]. 


al  e-field  vs.  height  axial  e-field  vs.  height 


ire  3.  Scattered  axial  electric  near-field  over  the  source  region  for  the  monopole  of  Figure  2. 
H  =  10  meters.  M/X  ■=  0.1.  CSD  source  with  V,„  =  1.0  volt.  Twenty  uniform 
segments  with  thin-wire  kernel. 


Figure  4.  Vertical  thin-wire  monopole  with  Applied  Electric  Field  (AEF)  source  model 
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Discrelizalion  of  ( 1 )  begins  by  expanding  the  current  I(z)  in  terms  of  a  weighted  set  of  basis 
functions.  The  bases  are  chosen  to  most  accurately  represent  the  current  at  all  points  along  the 
wire.  Within  NEC  the  wire  is  divided  into  N  segments,  of  length  A„  .  The  current  on  each  segment 
is  represented  by  a  weighted  three-term  basis  function  such  as  that  given  by  Equation  (6),  for  the 
n  segment 

y„  (r  )  =  A,,  +  sin  k  (z  — r,,  )  +  C„  cos  k  (z  —z„  )  /z  —z„  I  <  Ln  H  (2) 


where  z„  is  the  center  of  the  n.'*'  segment  and  k  =  is  the  free-space  wave  number.  Each 

basis  function  contains  three  coefficients,  two  of  which  can  be  eliminated  by  enforcing  continuity  of 
the  current  (Kirchoffs  Law)  at  wire  junctions  and  a  separate  condition  derived  by  Wu  and  King  [5] 
which  approximates  the  distribution  of  charge  at  wire  junctions.  For  thin  wires  of  constant  radius, 
the  charge  is  continuous  across  segment  junctions.  Since  the  axially  directed  current  cannot  change 
abruptly  to  zero  at  wire  ends,  a  separate  condition  given  by  Equation  (3)  is  enforced  which  relates 
the  current  to  its  derivative  at  free  ends  [l]. 
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where  J„  and  J  j  are  Bessel  functions  of  order  zero  and  one.  Once  the  junction  and  end  conditions 
have  been  applied,  it  remains  to  .solve  for  the  remaining  N  coefficients.  This  is  done  first  computing 
the  inner  product  of  both  sides  of  the  integral  equation  for  the  expanded  current  with  a  set  of  N 
testing  functions.  In  NEC  the  testing  functions  are  impulses  located  at  the  center  of  each  segment. 

We  now  have  a  set  of  N  equations  for  the  N  unknown  coefficients  used  in  the  current 
expansion. 
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where  are  ihe  unknown  coefficients  and  /„  the  unknown  basis  functions  on  each  segment.  For 
lumped  loading,  the  impedances  Z„  are  the  loading  impedances  which  are  assumed  to  be  located  at 
the  center  of  each  segment,  and  therefore  contribute  only  to  the  self  impedance  or  diagonal 
elements  on  the  right-hand  side  of  (4).  The  coefficients  a„  are  computed  by  inverting  Equation  (4) 
and  using  a  procedure  outlined  in  [l]. 

The  use  of  impulses  as  testing  functions  is  known  as  point  matching.  The  advantage  of  using 
point  matching  in  conjunction  with  a  non-singular  approximation  to  the  free-space  Green’s 
function,  such  as  the  thin-wire  kernel,  is  that  the  inner  product  can  be  performed  simply  and 
analytically,  thus  making  for  more  efficient  use  of  the  computer.  The  disadvantage  of  using  point 
matching  is  that  the  discretized  operator  equation  (4)  is  only  solved  at  the  match  points  and  for 
irregular  distributions  of  current,  significant  sampling  error  can  result.  It  is  therefore  highly 
recommended  that  as  a  check  for  accuracy,  the  axial  near  field  be  computed  at  all  points  on  the 
antenna,  from  the  computed  current  distribution,  and  analyzed  to  be  certain  that  the  boundary 
condition  that  the  tangential  electric  field  vanishes  on  the  conductive  surfaces  away  from  the 
source  and  loads  is  satisfied.  As  a  rule  of  thumb,  a  density  of  at  least  ten  segments  per  wavelength 
will  help  ensure  that  sampling  error  because  of  point  matching  is  not  significant  for  most  cases. 

In  many  cases,  the  condition  that  the  tangential  electric  field  vanishes  on  the  surface  of 
perfectly  conducting  wires  is  not  satisfied  in  the  vicinity  of  free-wire  ends,  even  when  Equation  7 
is  applied.  This  is  illustrated  for  a  loaded  antenna  in  Figure  6.  The  reasons  for  this  are: 

1 )  Fields  due  to  currents  on  the  wire  ends  were  not  taken  into  consideration. 

2)  Computational  errors  such  as  roundoff  can  also  lead  to  inaccurate  results. 

In  most  cases,  the  latter  is  of  little  consequence,  since  word  lengths  in  excess  of  32  bits  are 
commonplace  on  most  modern  computers.  For  thin  wires,  it  should  be  noted  that  the  axial  current 
should  be  very  small  in  the  vicinity  of  free-wire  ends:  hence,  the  fact  that  the  boundary  condition 
on  the  tangential  electric  field  is  not  exactly  satisfied  near  the  free  end  is  usually  of  little  or  no 
consequence  when  computing  either  the  input  impedance  or  the  radiation  field.  However,  one 
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AXIAL  E-FIELD  VS.  HEIGHT 


HEIGHT  IN  METERS 


Figure  6.  Scattered  axial  electric  near-field  over  the  entire  axis  for  the  antenna  of  Figure  4  with  a 
2500  resistive  point-load  2.5  meters  above  the  base.  H  =  10  meters.  V,„  =  1.0  volt. 
H/X  =  .l.andH/a  =  100.  Twenty  uniform  segments  with  thin-wire  kernel. 
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should  be  very  cautious  about  modeling  a  load  near  a  free-wire  end.  In  this  case,  there  may  be  a 
significant  current  in  the  load  which  may  lead  to  erroneous  results  in  the  presence  of  the  non-zero 
tangential  electric  field  near  the  end. 

3.  COMPUT.ATION  OF  SYSTEM  PAR.\METER  OF  WIRE  ANTE.WAS 

By  solving  Equation  (4)  one  obtains  an  approximation  for  the  current  along  the  axis  of  the 
antenna;  however,  knowledge  of  the  current  by  itself  is  of  little  practical  value  to  the  antenna 
engineer.  The  parameters  which  are  of  interest  are  those  which  are  used  to  evaluate  the  antenna  in 
an  operating  environment,  such  as  input  impedance  and  radiative  efficiency.  -All  of  these 
parameters  are  related  to  and  computed  from  the  current  distribution;  hence  one  must  exercise 
caution  to  prevent  errors  in  the  computed  current  distribution  from  manifesting  themselves  in  the 
aforementioned  parameters.  What  follows  is  a  brief  discussion  of  some  techniques  that  can  be  used 
to  compute  these  parameters  along  with  an  assessment  of  their  sensitivity  to  small  variations  or 
errors  in  the  computed  current.  Once  again,  results  obtained  using  NEC  will  serve  as  the  basis  for 
comparison. 

3.1  {COMPUTATION  OF  THE  INPUT  IMPEDANCE  OF  THE  THIN  WIRE  ANTENNAS 

One  of  the  most  widely  sought  parameters  for  antennas  with  a  two-terminal  feed  is  the  input 
impedance  looking  into  the  feed  port.  This,  after  all.  relates  the  antenna  as  a  radiating  element  to 
the  antenna  as  a  circuit  element.  The  input  impedance.  .  to  a  two-terminal  network  as  shown  in 
Figure  7  is  defined  as  the  ratio  of  the  input  voltage.  V,„  .  to  the  input  current.  .  in  the  limit  as  the 
spacing  between  the  two  terminals  approaches  zero  [2].  Previously  it  was  shown  that  for 
computational  purposes,  that  the  voltages  across  the  input  port  must  be  modeled  using  a  field- 
theory  generator  as  an  axial  electric  field  incident  over  a  very  small  portion  of  the  structure.  By 
invoking  the  principle  ol  reciprocity  (2].  [6],  one  can  show  that  the  relation  between  the  terminal 
voltage  and  current  and  the  field-theory  source  is 


l^-mes  939 


ntOCEEMHGS  OF  THE  1989  AHTEWM  APPLICATIONS  SVHPOSIUH 
VOLUHE  2(U>  ROHE  AIR  DEVELOPHENT  CENTER  GRIFFISS  AFB  NV 
DEC  89  RADC-TR-89-242-V0L-2 


2/4 


UNCLASSIFIED 


r/6  9/S 


V  s 


(5) 


w-;.  ---w 


\\„  I,„  =  f  E,iz)J{z)dz  (5) 

where  E,  (r  )  is  ihe  axial  electric  field  of  the  field-theory  source  in  ( 1 )  and  I(z)  is  the  true  current  in 
the  feed  region.  Integration  is  confined  to  the  portion  of  the  antenna  over  which  I  0.  By 

invoking  reciprocity  once  more  we  obtain 


/  E,  (z  )I  (z  )dz  =  -  f  E,  {z  )/  iz  )dz  +  f  r-iz  )Z  (z  )  dz  (6) 

A„  =-  f  EJz)  I(z)dz  +  J  lHz)/{z)dz  (7) 

where  L  is  the  axial  scattered  field  due  to  the  computed  current.  I  (z  ).  and  Z(z)  represents  any 
loading  impedance  on  the  antenna.  Integration  is  over  the  portions  of  the  antenna  over  which 
I/;  I  ^0.  The  latter  equation  is  of  utmost  importance  since  it  represents  the  problem  which  was 
actually  solved,  whereas  (5)  represents  the  problem  that  would  have  been  solved  if  (1)  could  have 
been  solved  exactly.  From  C7)  one  readily  obtains 


Z,„  =  /  E,  (z  )/  (z  )-  iKz  )7  (r  )  dz 


where  Z,„  is  the  input  impedance  looking  into  the  antenna  terminals  and  l,„  is  the  current  at  the 
feed.  Through  the  calculus  of  variations  it  can  be  shown  that  (8)  is  highly  insensitive  to  small 
errors  in  the  computed  current  and  thus  is  an  accurate  representation  of  the  input  impedance.  If 
the  feed  region  is  modeled  as  a  perfect  conductor  and  the  computed  scattered  axial  electric  field  is 
indeed  highly  localized  to  the  immediate  vicinity  of  the  source,  as  is  often  the  case  with  NEC.  (8) 
reduces  to  the  familiar  form 


Z..  =  -i  f  E,(.z)Hz)dz 
I.r. 


where  integration  is  confined  to  the  feed  region  over  which  \E,  I  ^  0.  Due  to  the  added 
complexity,  one  does  not  always  wish  to  compute  Z,„  using  (9)  and  approximations  to  this  integral 
are  frequently  sought.  One  approximation  that  is  frequently  used  is  to  take  the  ratio  of  the  applied 
source  voltage  of  the  field-theory  source  to  the  computed  current  at  the  center  of  the  feed  region. 
This  approximation  is  used  in  NEC  and  the  conditions  for  its  validity  will  be  derived. 
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It  was  staled  previously  that  Z,„  =  — — .  in  the  limit  as  the  spacing  between  the  input 
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terminals  of  the  antenna  approaches  zero.  If  such  is  the  case,  the  input  voltage  is  confined  to  an 
infinitesimally  small  region  over  which  the  current.  Kz).  is  constant  over  the  feed  and  equal  to  /,„ . 
In  the  case  of  the  applied  field  (Alii  )  .source,  the  computed  current.  /  (2  ).  is  nearly  constant  over 
the  region  that  I/-.',  I  O;  hence.  (9)  can  be  approximated  by 


I,„  f  (7 )  d: 


(10) 


The  computed  current  and  near  field  for  the  AEF  source  are  shown  in  Figure  8  for  vertical 
monopoles  with  H/a  =  KKf  and  KKK).  The  integral  in  (10)  can  also  be  interpreted  as  a  voltage.  V',  . 
where 

\\  =  f  !■:,(:)  d:  (11) 

If  the  numerical  solution  is  a  good  approximation  to  (1)  then  V,  =  —V',,  .  where  V,„  is  the  terminal 
voltage  specified  by  the  user  in  NL(i.  and  (10)  reduces  to 


where  /,„  is  specified  in  NEC  as  the  current  at  the  center  of  the  segment  over  which  the  AEF  source 
is  applied.  In  Table  1.  a  comparison  is  made  between  Z,„  computed  using  (12)  and  by  numerically 
integrating  (9)  for  unloaded  and  resisiively  loaded  vertical  monopoles  of  H/a  »  KKK).  100.  The 
ground  was  modeled  as  a  perfect,  infinite  ground  plane.  The  results  presented  here  indicate  a  very 
close  agreement  between  the  input  impedance  of  thin-wire  antennas  computed  using  the  stationary 
form  of  (9)  and  the  NEC  approximation  (12).  In  each  case,  the  current  was  nearly  constant  over 
the  source  region  and  the  axial  electric  field  integrated  to  -1.00  volts  over  this  region.  The  assumed 
input  voltage  was  1 .00  volts.  The  .source  region  for  this  problem  w'as  taken  as  the  lower  most  three 
segments,  or  1 .5  meters,  of  the  antenna.  Beyond  this  region.  E,  iz  )  has  negligible  magnitude  as 
indicated  by  the  plot  in  Figure  6. 
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computed  using  NEC  approximation  (12)  and  stationary 
ar  th#  first  thraa  segments  (1.5  m).  H  =  10  m,  H/X 


A  similar  comparison  is  shown  in  Table  2  for  the  same  antennas  with  the  CSD  feed.  The  axial 
near  field  and  the  corresponding  current  over  the  feed  region  are  shown  in  Figure  9  for  vertical 
monopoles  with  H/a  =  100  and  1000.  respectively.  When  H/a  =  100.  the  current  is  not  constant 
over  the  feed  region  and  the  scattered  electric  field  exhibits  irregular  behavior  there  as  well.  The 
integrated  voltage  over  the  "main  lobe"  of  the  electric  field  is  -l.lO  volt.  If  we  include  the  second 
and  subsequent  lobes  in  the  feed  region  the  integrated  voltage  is  -0.80  volts.  The  input  voltage  is 
still  assumed  to  be  1.00  volts.  Due  to  the  erratic  behavior  over  the  feed  region,  we  do  not  expect 
good  agreement  between  the  input  impedances  computed  using  the  NEC  approximation  and  the 
stationary  form.  This  is  evidenced  in  Table  2.  The  erratic  behavior  over  the  source  region  for  H/a 
=  100  can  be  explained  in  part  by  noting  that  the  discontinuity  in  the  charge  due  to  the  CSD  source 
model  is  accurate  only  for  filamentary  wires  with  infinitesimally  small  feed  regions  [l]  and  [3]. 
For  H  a  =  KKK).  the  feed  region  is  much  more  confined  and  the  current  is  nearly  constant  over  this 
region  as  shown  in  Figure  9.  Hence,  as  one  would  expect,  agreement  between  the  NEC 
approximation  to  Z„  (12)  and  the  stationary  form  (9)  is  better  for  H/a  =  1000. 


3  2.  CO.MPL'TATION  OF  THE  RADIATIVE  EFFICIENCY  OF  THIN- WIRE  ANTENNAS 

The  second  parameter  of  concern  in  this  paper  is  the  radiative  efficiency  of  wire  antennas.  The 
radiative  efficiency.  7).  is  defined  as  the  ratio  of  the  average  radiated  power  to  the  average  input 
power  of  the  antenna. 


In  \f:('  the  input  power.  P,„  .  is  computed  using  the  formula 

C  =  y  Ke[V,„/,„‘]  (14) 

w  hich  assumes  that  ihe  feed  region  is  infinitesimally  small.  The  radiated  power  is  computed  as 

=  P,n  -  Tu...  (15) 
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computed  using  NEC  approximation  (12)  and  stationary 
'er  the  first  three  segments  (1.5  m) .  H  «  10  m,  H/ A 
lire  kernel.  E  and  I  sampled  at  40  points  per  segmei 
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where  7’,,,^  is  the  total  power  dissipated  in  the  structure  [l]  The  pouer  dissipated  in  the  structure 
is  computed  directly  as  the  sum  of  the  average  power  dissipated  in  each  of  the  M  loads  in  the 
structure  as  stated  by  liquations  ( 16)  and  (17). 

=  Z  (16) 

■1  =1 

where 


P„,  =1  I/.,  R„ 


(17) 


and  I,„  is  the  current  through,  and  is  the  real  part,  or  resistance  of  the  m'*  load.  Equations 
(13)-(  17)  can  be  combined  to  yield 


Re  IV„  ‘  J 


(18) 


Computations  of  the  efficiency  using  (18)  are  inherently  sensitive  to  errors  in  the  assumed  input 
voltage  and  current  unless  the  current  over  the  feed  region  is  nearly  constant  and  the  integral  of 
the  scattered  electric  field  over  the  feed  region  equals  -V,„  .  Such  is  usually  the  case  in  NEC  when 
the  -AEF  source  is  used  and  when  the  CSD  source  is  used  for  wire  structures  with  H/a  ^  1000. 
Four  alternate  methods  of  computing  the  efficiency  are  presented  below.  All  of  these  expressions 
involve  integration  of  the  current  and  thus  are  inherently  less  sensitive  to  small  errors  in 
computing  the  current. 


A  statement  of  conservation  of  energy  is  obtained  by  considering  the  real  parts  of  the  complex 
power  of  both  sides  of  ( 1 ). 


1  Re  <  /(,  .  /*  >  =  -1  Re  <  £, .  /»  >  +  1  Re  <  /  ^^  .  /»  > 

2  2  2 


(19) 


Stated  more  simply,  the  input  power  equals  the  sum  of  the  radiated  power  and  the  power 
dissipated  in  the  siriisture  The  input  po\n  er  is  given  by 
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P,.  =  ^  Re  <  £, .  /*  > 


=  ±Re 


f  £,(:  )/*  (2) 


the  radiated  power  is  given  by 


P..,  =  -y  Re  <  £, .  /*  > 


=  -yRe 


H 

f  £,(.2)  P  (2  )dz 


(20) 
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and  the  power  dissipated  in  the  structure  is  given  by 


Po..  =  y  Re  <  /Z  .  /♦  > 


=  ±Re 


/  l(z)P(.2  )Z  iz)dz 


(22) 


When  point  loads  are  used.  (22)  reduces  to  the  discrete  summation  of  (17).  Since  (1)  could  not  be 
solved  exactly,  computation  of  the  input  power  using  (20)  is  ill-advised.  As  an  alternative,  the 
input  power  should  be  computed  for  the  problem  that  was  solved  approximately  by  inverting  (4). 
This  is  done  in  Equation  (23) 


/ 
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E,(.z)P(z)-Hz)P{z)Z(z)  dz 


(23) 


where  the  tilde  indicates  that  the  input  power  is  computed  using  the  current  obtained  by  solving 
the  discretized  Equation  (4).  If  there  are  no  loads  over  the  source  region  and  £,  is  localized  to  the 
source  region.  (23)  reduces  to 


f  £Az)r(z)dz 
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(24) 


II  the  current  is  constant  oser  the  source  region  and  the  integral  of  the  near  electric  field  over  this 
region  equals  -V’,  then  (24)  reduces  to  (14),  thereby  validating  the  assumption  of  a  point  feed 


'A  hen  compuiing  the  radiative  efficiencj.  .  I'he  conditions  for  the  validity  of  the  assumption  of  a 
point  feed  are  identical  whether  the  input  impedance  or  radiative  efficiency  is  to  be  computed. 

When  (24)  is  substituted  inti'  (IS)  for  P,.,  .  an  expre.ssion  for  the  efficiency  is  obtained  which 
assumes  pointv.  ise  loading  hut  makes  no  restrictions  on  the  feed  other  than  that  it  is  localized  to  an 
unloaded  portion  of  the  antenna.  Since  the  e.xpre.ssion  obtained  using  this  procedure  involves 
integration  ot  the  current  it  is  inherently  less  .sensitive  to  small  variations  of  the  computed  current 
about  its  true  value. 

The  second  method  utilizes  direct  computation  of  the  average  radiated  power  using 


where  t),.  12(>7r  c'hms  is  ihe  inlnnsic  impedance  of  free  space  and 
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is  the  electric  tar  field  from  a  z-direcled  current  at  a  distance  R  from  the  antenipa  [3].  Integration 
over  6  in  (25)  is  from  0  to  •jr/2  if  a  ground  plane  is  present  as  is  the  case  for  the  Vertical  monopoles 
described  herein.  The  efficiency  could  then  be  computed  by  substituting  (25)  and  (17)  into  (15) 
and  (13)  to  obtain 


^  -  1  T  (/v„.77^:,T 

A  third  method  of  computing  the  efficiency  v/hich  is  also  less  sensitive  to  computational  errors  in 
the  current  involves  computation  of  the  input  and  radiated  power  directly  using  (24)  and  (25)  and 
substituting  these  expressions  directly  into  ( 13). 

In  the  fourth  method  presented  herein  for  computing  the  radiative  efficiency,  the  radiated 
power  is  computed  directly  by  computing  the  inner  product  of  the  scattered  axial  electric  field  and 
the  conjugate  of  axial  current  I'ver  the  entire  siructure  as  shown  by  Fquation  (21).  The  radiative 
efficiency  is  i  hen  computed  b-,  subsi  n  ui  ing  (21  )  amf  ( 24)  into  (13)  to  yield 


'.'■.vJ-’-'j 


V 

s'  vV-' 

•x.v'c-c 


■  ■  ">'  ■  •  ■  •  '  •  ■  »*•  '.s','”  . 


(28) 


1)  = 


Re[  /  E,  (z  )/*  (z  )  dz  ] 
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'I'his  meihod  yields  an  accurate  compulation  of  the  radiative  eff\cienc>  for  the  modeled  antenna 
because  neither  point-sources  nor  point-loads  are  assumed  as  they  were  in  the  two  methods  staled 
previously  by  I-quaiions  (18)  and  (27).  In  addition,  both  the  input  and  radiated  powers  are 
computed  using  integrals  and  thus  are  inherently  less  sensitive  to  small  variations  of  the  computed 
current  about  its  true  value  The  only  constraints  are  that  the  source  region  be  localized  and  free 
of  loads.  If  more  than  one  source  exists  and  the  source  regions  are  not  free  of  loads,  a  more  general 
expression  for  the  efficiency  is 
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dz 


A  comparison  is  made  in  Table  3  of  the  efficiency  ol  resisiively  loaded  whip  antennas  of  H  a  = 
100  and  KiOO  using  liquation  (28)  and  the  efficiency  computed  using  the  .NliO  approximation  (18) 
for  both  the  AFT  and  CSl)  source  models.  These  results  show  that  the  .\FC  approximation  is  valid 
when  the  AEF  source  is  used  The  NFC  approximation  is  shown  to  yield  marginal  results  for  the 
CSD  source  when  IT'a  =  1000  and  even  le.ss  acceptable  results  when  H/a  =  100. 

1\  ('ONCl  I  SKINS 

It  IS  shown  herein  that  w'hen  wire  antennas  are  sufficiently  thin,  some  simplifying 
assumptions  cun  l)e  used  m  conjunction  with  the  meihod  ol  moments  to  obtain  a  set  of  equations 
for  the  unknown  structural  current  which  can  be  .solved  with  relali\e  ease  using  a  digital 
computer  When  the  input  impedance  and  efficiency  are  computed  using  simple  expre.ssions  that 
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lative  efficiency  of  lO-meter  reslstlvely-loaded  verticol  whip  above  perfect  ground  computed  using  AEF  and  CSO  sources 
h  the  NEC  approximation  (18)  and  the  less  sensitive  form  (28).  Numerical  Integration  performed  using  E  and  I  sampled 
points  per  segment.  250(1  point  load  located  2.75  meters  above  the  base.  Twenty  uniform  segments  with  thin-wire  kernel 


assume  a  poiniwise  feed,  as  in  NEC.  these  two  parameters  are  very  sensitive  to  errors  in  the 
tompuled  current  Alternate  formulas  for  both  the  input  impedance  and  the  efficiency  that  are  less 
sensitive  to  small  variations  in  the  current  about  its  true  value  are  presented.  It  is  shown  that 
hen  the  modeled  feed  region  is  small,  the  computed  current  is  nearly  constant  over  the  feed 
region,  and  the  integrated  near-electric  field  over  the  feed  region  equals  the  negative  of  the  input 
\oltage.  the  alternate  expressions  reduce  to  the  simple  forms  used  in  \EC.  Comparisons  between 
the  input  imp>edance  and  efficiency  of  thin-wire  whip  antennas  computed  using  \EC  alone  and  NEC 
in  conjunction  with  the  alternate  forms  are  presented  herein.  Eor  the  simple  geometry  discussed 
herein,  it  is  shown  that  NIiC.  alone,  yields  highly  accurate  results  when  the  applied  electric  field 
(.AEF)  source  is  used  in  conjunction  with  wires  with  H-a  ^  100.  The  current-slope-discontinuity 
(CSD)  source  is  shown  to  yield  marginally  acceptable  results  when  used  to  model  wire  antennas 
with  H/a  <  1000. 

The  agreement  or  disagreement  between  the  data  sets  is  shown  to  be  directly  related  to  how 
well  the  aforementioned  criteria  defining  the  validity  of  the  simple  expressions  used  to  compute  Z,„ 
and  7)  in  NEC  are  met.  It  is  highly  recommended,  therefore,  that  one  should  closely  examine  the 
computed  current  and  near  field  over  the  surface  of  the  modeled  structure  and  then  determine 
what  simplifying  forms,  if  any.  are  valid  when  computing  the  terminal  properties  of  antennas. 
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I  his  paper  is  concerned  with  the  analysis  ol  a  plane  wace  incident  upon  a  slab 
of  lossy,  dielectric  honeycomb.  Two  methods  are  used  in  the  analysis.  The  first 
replaces  the  structural  medium  with  a  simple  dielectric  using  effective  medium 
theories.  'Lhe  second  is  a  lull  wave  analysis  that  takes  the  structure  into  account. 
These  methods  are  compared  to  each  other  and  to  experimental  results. 

1 .  Introduction 

A  lrei(iiciu  \  select  i\e  surLice  used  in  .satellite  applications  is  structurally  .sup- 
porteii  b>  a  dielectric  slab,  as  shown  in  L'igure  1.  that  must  be  light-weight  yet 
strong.  Structured  dielectrics  are  the  preferred  support  material  since  they  have  a 
higher  strength  to  weight  ratio  than  do  simple  dielectrics.  One  such  structured 
dielectric  slab  is  the  dielectric  honeycomb.  .As  show  n  in  I  igure  2.  the  honeycomb 
backing  allect'-  the  characteristics  of  the  frequency  seleciice  surface  by  lowering  the 


k. 


trequency  ol  transparency  by  approximately  lO^'i.  It  is  therefore  important  to  be 
able  to  characteri/e  the  behavior  ol  an  electromagnetic  wave  in  the  presence  of  a 
slab  of  honevcomb  material. 
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Figure  1.  F.S.S.  (left)  supported  by  a  dielecirit  honeycomb  (right) 


In  this  paper,  methculs  are  vuilineil  in  analyze  an  eleclroniagnelic  plane 
wave  intideni  on  a  slab  of  lossy,  dielectric  honeccomb.  The  first  method  is  very 
simple.  It  uses  a  quasi  static  model  in  order  to  remove  the  structure  in  the  honey¬ 
comb  and  replace  it  with  a  homogenecuis  dielectric  slab  having  an  efleclive  permit- 
tivitv.  This  can  be  easily  cascaded  with  the  characteristics  ol  the  frequency  selec¬ 
tive  surface.  The  second  method  is  a  full  vcave  solutic'n  and  is,  therefore,  more 
complex.  It  is  based  on  an  application  ol  I'loquet's  theorem  and  the  generalized 
scattering  matrix.  These  two  methods  cc  ill  Iv  compared  to  each  other  and  to  experi¬ 
mental  results. 

2.  Problem  (ieometry 

The  geometry  of  the  problem  under  consideration  is  shove  n  in  figure  3.  A 
plane  wave  is  incident  on  a  .slab  of  dielectric  honeycomb  which  is  infinite  and 
periodic  in  both  the  x  and  the  y  directions  but  finite  in  the  /  direction.  The  axes  of 
the  honeycomb  are  aligned  along  the  z  axis.  The  unit  cell,  shown  in  f  igure  4.  is  the 
structure  that,  when  repeated  in  the  x  and  y  direction,  gives  the  honeycomb  its 
periodicity.  'The  clireclion  of  propagation  lor  the  incident  wave  makes  an  angle  6 
w  ith  the  positive  x  axis  and  (f)  with  the  positive  y  axis  Since  an  arbitru  y  field  can 
be  broken  into  a  superposition  of  wav'es  transverse  electric  (  T1  )  and  transverse 
magnetic  ('TMl  to  a  given  direction,  the  scattered  field  response  is  calculated  for  an 
incident  field  'TM  to  x  and  'Tf  to  x. 

.3.  1  Hectic e  .Meilium  .Approach 

This  method  was  useil  to  gain  an  initial  insight  into  how  the  honeycomb  slab 
behaves.  In  past  work  with  resistive  strips,  it  has  Iteen  established  that  at 


sullkienilv  lovv  frequency,  a  structured  medium,  such  as  a  honeycomb,  exhibits 


cfiaracleristics  similar  to  a  simple  dielectric  with  constitutive  parameters  equal  to  a 
weijL’hteil  a\era^e  I'l  the  parameters  of  the  tumeycomb  walls  and  the  parameters  of 
the  tree  space  that  tills  the  honeccomti  1  urther.  we  expect  a  slight  amount  of 
anisotropic  hehac  lor  fe  >€,  =€,  )  sirite  tlie  walls  <>1  the  honeyci'mb  are  parallel  to  z 
and  the  electric  field  in  the  /  direction  t/-.  )  is  more  strongly  influenced  by  the 
dielectric  than  A,  or  . 

\  number  ol  different  ways  of  averaging  the  constitutix  e  parameters  have 
been  proposed  in  literature  tor  application  to  random  inhomogeneous  media  such  as 
smoke,  hail  and  metallic  powders.  We  propose  to  apply  three  of  these  effective 
medium  theories  ( Maxwell-Clarnet.  Bruggeman.  and  Maxwell-darnet  for  spheroids) 
to  the  honeycomb  to  obtain  an  effective  permittivity.  In  practice,  one  often  applies 
a  number  of  these  theories  to  the  problem  at  hand  and  determines  the  validity  of  a 
particular  theory  by  comparing  the  theoretical  results  to  experimental  measure¬ 
ments. 

rtie  honeycomb  can  tie  viewed  as  being  compo.scd  of  strips  (labeled  the  inclu¬ 
sion)  with  a  permittivity  embedded  in  a  matrix  of  free  space  with  a  permit- 
tivitv  €  as  slu'wn  in  I  igure  5. 


^Inc 

f  igure  5.  Definition  of  the  inclusion  and  matrix 
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I  l^^ure  />  Isolaletl  inclusions  requiring  the  application 
ol  the  MawceirCiarnel  iheorv 


C\r\ 


I  li^uro  7  Reduction  of  cells  in  the  honeycomb  lo  equivalent  spheroids 


Simihirlc  the  .f.  eraije  polarization  is 

</’,„(r)>  =(!-/  )<r„,,{r)>  +  f  </5,,„(r)>  (2) 

Assuming  that  the  matrix  and  inclusions  are  lx)lh  linear  and  isotropic,  the  average 
/•-'  and  average  /’  are  related  by 

</’,„>  =e„X  <F:„„,>  (3) 

</V.  >=  €..\  ,  </?,„.  >  (4) 

<r  ,,  >=  €,\,  ,  >  (5) 

where  V. »  defined  as  the  axerage  susceptibility  ten.sor.  (Combining  equations  (1) 

l  hrough  ( .“i )  we  obtain 


•  ."xi 


s-iv.vi 


V  si 


(1-/  )(?,  -  e  T)  (r)>  +  /  '€.<  -  €  ,  T)  </7„,  (r  )>  =d  (6) 

Iki-ause  Ihe  o.er.ill  meilium  is  assume!  lu  l)e  sUiiistically  liomogeneous.  we  must 
obtain  a  linear  relationship  (jet'^een  ,  >  am!  <A,,,  >.  II  we  look  at  an  iso- 

lateil  splteroiii  in  a  unil  orm  liekl  ean  >;et  tlie  relationship  [3] 


i 

■m 

--i 


</  .  >  =  K  <  }  ■:  .  > 


w  here 


and  I.  are  the  ilepolarizalion  factors  [-1],  Maxwell-Ciarnel  for  a  s'pheroid  becomes 

(1-/  -  €„„„  T)  T  +  /  ^  T)  X  =  0  (9) 

II  the  spheroid  becomes  a  sphere  then  /.^  =  1/3  and  the  above  expression 
reduces  to 

€,  »  -  €„, ..  -  e,. ,, 

which  is  the  original  Maxwell-darnet  expression  for  the  efTective  permittivity  in 
terms  of  the  fill  factor  and  the  permittivity  ol  the  inclusion  and  the  matrix. 
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3  2  llruggeman  I'heorv 

Hruggeman  theory .  developed  m  1935.  applies  to  the  aggregate  structure 
shown  in  l  igiire  X. 


1  igure  X  .Aggregate  structure  requiring  the 
application  of  Hruggeman  theory 


is 


rile  mt I Li'i  li  e  ni.t'nx  niu^i  ''<■  tii-.ileil  \v,ilii  '•uui.Uiral  equivalence  so  the 

unit  cell  1"  r'  lliicisi  i-  i  '•p!  ere  .1  railin''  .1  iinmerseii  m  a  ilielecinc  w  ith  a  permit 
uv  it  .  ot  £  ..  Tile  ''phere  iut'-  a  prchahilii ot  t  i|';il  laetor)  that  Us  permiUiv  ity  is 
e  aiiil  1  1  lh.it  its  perm itt iv  It V  ic  £  l  lie  dipole  moment  of  each  of  these 

spheres  is 


P. 


=  Atra  *€.  f,  - 
P,,  „  =  47rfl  '€,  ^ 


£  +  2€. 

-e, - 


I- 


(11) 


(12) 


+  le,.f 

L'sinij  tlie  delinition  cl  an  efleclive  medium,  the  dev  iation  produced  by  the  polariza 
tion  must  vanisti  in  'he  average. 


//V  +(!-/  )a. 


(13) 


or 


/■ 


(€„,  +  2€,,  )  +  2€,./1 

w'hich  is  the  Mruggeman  expression  for  an  efleclive  permilliviiy. 


(14) 


3.3  Small  ('ontrast  Limit 

11  the  contrast  between  the  inclusion  and  the  matrix  becomes  small  i.e. 
A  =  €,.  —  £  <  1  then  all  ol  the  aliove  theories  reduce  to  the  relationship 

€.f  =  /  +  (1-/  )e  ,,  (15) 

which  is  just  the  average  til  the  two  components. 


3.4  I  llective  Medium  Melhcuf:  Conclusion 


riie  advantage  to  this  ap,  'oach  is  that  it  permits  the  use  of  relatively  simple 
theoretical  lormulas  as  compared  to  those  derivetl  Irom  a  complex  boundary  value 


problem.  Ihe  seconti  aii^;inuii,y  is  thai  Uie  results  yieUi  a  simple  dieleetru  slab 
\v)iit)i  IS  eas\  to  cascade  w  ith  1  ret[uenc\  selei.ti\e  surlaces  I'he  disatl vantage  is  that 
this  method  does  not  sltov.  hi'w  the  striuiure  iisell  alleets  tlie  \^aves.  lor  example 
to  what  extent  cross  i.oupltng  exists  betxxeen  11  anil  1  \1  inculent  xcaxes.  l  or  this 
information  \xe  must  turn  to  the  more  complex  full  \xaxe  analxsis 

4.  I'ull  \\'ave  Anaix  sis 

In  order  to  obtain  a  full  wave  solution  lor  the  honeycomb  slab,  an  integral 
equation  will  be  formulated  where  the  unknown  quantity  is  the  total  e’ectric  field 
across  the  unit  cell  showm  in  l-'igure  4  This  equation  w  ill  then  be  solved  approxi¬ 
mately  by  the  method  of  moments. 

4.1  Integral  liquation  formulation 
We  know  that  in  general 

/.■"  (r  )-/?'  (T)  =  h  '  (F)  (16) 

Where  /•.  is  the  unknown  total  field.  /•  ’  is  the  unknown  fielil  scattered  by  the 

dielectric  and  /.  is  the  kni  vxn  plane  wa\e  iiiiideni  upon  the  honeycomb.  We 
must  now  find  a  relationship  between  /:  and  /: '  1  irst  the  dielectric  in  the 

problem  is  replaced  b\  an  etjuivalent  polan/ation  current.  J^\  .y  .z  ).  radiating  in 
free  space.  '  (F  )  max  (hen  be  found  bv  loncnlx  ing  J  (x  .z  )  over  Ihe  unit  cell 
with  the  periodic  (ireen's  lunation  (i  (.\  .y  .r  '.v’.y'.r’)  which  is  defined  as  the  elec¬ 
tric  field  prodticed  by  an  inlinite  array  ol  periodic  point  sources  radiating  in  the  x-y 


I  o  ilen.t*  I  Ik-  perioclii.  (ireen  lunttion  ^v•e  note  lluu  the  geometry  is  perioiiK 
m  X  and  \  (period  ol  c  and  d  respet  lively )  and  that  the  incident  fiehi  is  a  plane 
wave  with  propagation  constants  ol  .k,  .k  .  I'loquet  c  theorem  slates  that  the 
fields  in  the  \  -y  plane  are  periodic  in  x  and  y  except  l('r  the  phase  shill  introduced 
by  k,  and  k  of  the  incident  wave  We  can.  therefore,  dehne  a  I  oiirier  series  for 
the  function  Kx.y)  bv 


I  'I  =  L  L  /■' 


(IS) 


where 


‘ . 


dx  d\ 


and 


_  ,  ^  Inm  ,  ,,  , 

a.„  =  *r ,  +  — —  m  =  — oi . —1,0.1 . oo 

c 


i3„  =k,  + 


2vn 


n  =  —CO . —1 ,0.1 . c!o 


(19) 

(20) 
(21) 


The  array  of  point  sources  located  at  .v  '.y '  in  each  unit  cell  and  represented  by 


(-V  .y  )  = 


^  21  '—me  )6(y  — y  '—nd  )S(r  — r  ')  (22) 


can  be  equivalenilv  represented  as  an  infinite  set  of  current  sheets 


7.,  (.V  .V  )  = 


cd 


I  I 


(23) 


through  the  use  ol  the  l  ourier  series 


The  electric  field  response  ;ii  /  due  to  the  single  current  sheet 


(a  .y  )  =  — - -  k  y  x 

V/i,.i  ±oc,„y„.,  ±0  y,,^  o,* +0,;  - 


where 


+  'Jk  — a,"  ~0,r  it  ^  .T  ^  (of.n 't'0  r  ^ 

"  -y  V(ai+0,f)--t7'  .<*.:< (a, ^+0,-) 

A-,,  is  ihe  free  space  propagation  constant.  Rewriting  the  alx)'ve  equation  in  shorter 
form  vce  obtain 


llAx  .y  ,r  )  =  d{a„,  .0,  )  v  c  ^  '  e  -  '  '  (27) 

'I  he  response  doe  to  the  array  of  point  sources  (the  periodic  (ireen's  function)  is  a 
sum  of  tlie  respv'oses  to  the  current  slieels  according  to  equation  (27) 

(7(A,y,r  i  .x',  V.-)=  -L  f  f  a  (o...0.,)e^ ^Ji) 

,tl  =~-fjOll  =— Ok. 

The  double  summation  can  be  viewed  as  the  inverse  I’ourier  transform  to  return  to 
the  space  domain. 

Substituting  the  periodic  (Jreen's  function  shown  above  into  etjuation  (17)  and 
perli'i  niing  the  vcmcoluiion  we  obtain 


.  2  *  “  --OU//  =-«x. 


J  (.V  '.V  ' ‘)dx  './v 'd:  ' 


=  Z  z  (7  (q,„.0„  )i- Q,„  A  . "  'rfr'  (30) 

>*.  here  J  io  0  .r  ')  is  ihe  I  ourier  transform  of  the  current  in  the  x  y  direction.  We 
use  the  boundary  cc'ndilion  (7] 


J(x  .y  =  j  we.,{€,  (.V  ,v  )-l  )£""  ix  ,v  .r  ) 


J(a,  .0,,  ..  )  =  I\r.  |(€,  U  .y  f-l)/?""  (x  .y  --  )  |  (32) 

to  obtain  the  needed  relation  lietvxeen  J  and  .  I'.T.  is  the  two-dimensional  I'ast 
I  ourier  iranslorm  with  respect  to  a,„  and  0,,  .  .41so  note  that  from  0  to  Ar  .  £,  is 
invariant  in  y 

Substituting  ecjuation  (32)  into  (30)  and  this  into  the  original  equation  (16) 
we  obtain  the  tlesired  integral  equation. 


/?"■'  (x  ,v  .r  )  -  /(«€„  Z  Z  ) / I'T  (e,  (.V  '.y ')-!  )/"">'  (.t  '.y  '.z  ')  x  (33) 


f  '  ~  'd:  V '  ‘  '  ’  =  /?'■'•  (.V  .V  ) 


4.2  Solution  h\  Method  of  .Moments  [S] 


We  w  ill  expand  the  total  F.  fiehi  across  the  cell  as  a  sum  of  kl  basis  functions. 


F'"  {x  .y  .z)=  /7j 


liquation  (33)  Itecomes 


X 


(35) 


;  <*>«..  it  L  f  ^ '  •>■ )- J 


,  -.'  >/>,■  '  -  '  y  .  ‘  '>V,  •  +/<„  '  i 


Next  we  form  a  matrix  equation 


=  /?''■  (  a  .v.r  ) 


Za  =b  (36) 

by  performing  an  inner  product  on  both  sides  of  equation  (35)  with  kl  different 

testing  functions  T i  (x  .y  .z  ).  a  is  the  column  vector  of  unknown  coefficients 


a  =  (a  1  fly  iV  , 

6  is  a  column  vector  representing  the  forcing  function 

6,  =  <7“,  (a  .y  )^'"-  (a  ,v.r  )> 
and  /  is  the  impedance  matrix  where 


I  =  \  k] 


(37) 


(3S) 


Zn  =  <Ti(.\  .y  .z  )Si  {.X  .y  )>-J(i)€„<fi{.\  .y  ^  0  («,„  ,6e  )x  (39) 


f/r 


(e,  (v  ,V  )-l)«y  (a  ,v.z') 


A.  .  iz 


~  f  f  f  * '  '  '  ^  ^  ^  ~  Jo)€,fid  f  it 


1(  .  : 


/ft. 


(e,  (  a  .y  )— 1  )^y  (  v  .y  .z  ') 


j  r  (a,,.  A  .z  )dz  Wz 

and  w  here  T  (q„  .0,  .z  )  in  the  second  term  is  the  conjugate  of  the  two-dimensional 
1  ourier  iranslorm  of  the  testing  lunction.  The  tlouble  summation  can  now  tie 
vieweii  as  a  discrete  inner  protiuct  of  two  discrete  1  unctions. 

W  c  will  now  clioose  our  basis  function  to  be  a  1  loquet  harmonic  in  x  and  y 
and  a  pulse  in  /  of  length  Az  . 
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/};  ,  i contains  mlorniaiion  about  the  components  of  the  I.  held  normalized  to 
units  lor  the  particular  I  loquet  harmonic  chosen.  These  components  depend  on 
whether  the  harmonic  is  I  I  to  x  or  I'M  to  x.  Moth  types  must  Ire  included  in  the 
basis  functions.  I'lie  testing  function  is  chosen  idenlicalls  to  the  basis  function 
except  that  it  exists  only  in  the  dielectric.  This  truncation  serves  to  spread  out  T 
over  all  values  of  a,,,  and  /3„  . 

f  (.V  ,v  .r  )  =  f, ,  ,2(6  (a  ,v  )-l  )e' /’(Ar  )  (42) 

where 


€  .V  .V  )  = 


2  0  il  -V  .y  IS  li'cated  in  the  dielectric 

1.0  if  .X  .y  is  net  located  in  the  dielectric 


(43) 


With  the  ,ib(.'ve  choice  of  basis  and  testing  function,  term  1  in  the  impedance 
matrix  expression  (equation  40)  becomes 


i,<  2 


Ar  (44) 


Term  1  =  cdfii  1 ,  2  1 1  1 1  2  (•''  •>' 

where  I  I'.  is  the  I  ast  I  ourier  transform  ol  the  testing  function  (11,12)  evaluated  at 

the  mode  numbers  of  the  basis  function  (kl,k2).  Term  2  is  evaluated  by  perform¬ 
ing  the  integrations  in  z  analytically,  ('are  must  be  exercised  to  use  the  proper 
(7(a.,  .fi„  )  depending  on  the  position  of  the  observation  point  (z)  with  respect  to  the 
source  point  (c'l.  I  et  (/  =  for  c  >r'  and  G  —(.’ir  (o*"  ^  <-  ’•  Then  Term  2 
becomes 


Term  2  =  j  (u€„r./  £  I  1 1  2  ^.  ,  ^  r/'  |(e,  (.v  ,y  )-l  )f  ^ ''  ’ 
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riie  right  hand  side  of  the  expression  (36)  Iwcomes 


I  T  (€  (a  ,v  )-iV' 


r'-"y  2 

wliere  I'.T.  is  the  I  asi  1  ourier  iranslorm  ol  ihe  lestin;;  1  unclion  evaluated  al  ihe  0,0 

mode  niiml^er. 

4.3  Generalized  Seaiiering  Matrix  [9] 

L'pon  solution  of  the  matrix  equation  (.36)  for  the  coefTicienls  a  .  we  have  an 
approximate  expression  for  the  total  field  acro.ss  the  unit  cell.  Since  we  approxi¬ 
mated  the  field  in  the  z  direction  by  a  single  pulse,  this  solution  is  only  valid  for 
thin  slabs.  In  order  to  find  the  characleri.stics  of  a  thicker  slab,  we  must  cascade  a 
series  of  these  identical  thin  slabs  using  the  generalized  scattering  matrix. 

The  first  step  in  calculating  the  scattering  matrix  is  to  determine  the  scattered 
fields  on  either  side  of  the  slab  (al  0  and  al  Ar  )  by  substituting  (.v  .y  ,r  )  into 
the  expression  for  the  scattered  field  and  find  the  scattered  fiekls  lor  each  of  the  kl 
I'loquel  harmonics  specified  as  basis  functions. 

C"  (a  .y  .r  )  =  j  a>€.,/ {f(a,,,  .0,„  )c  ''  '  (47) 

A 

I  T.  (t,  (a  .y  )-l  )/7'-'  (a  .y  V ' '  ' 
w'here  m.n  are  the  mode  numbers  of  the  kl  basis  functions.  The  scattered  fields  on 
the  left  side  of  the  slab  are  the  reflected  fields  w  hich  lead  to  the  .S'n  submatrix.  The 
scattered  fields  on  the  right  side  are  added  to  the  incident  field  to  give  the  transmit¬ 
ted  field  which  lead  to  the  S^i  submairix.  I'he.se  fields  are  then  broken  up  into 
waves  Tl-  to  x  and  TM  to  x  for  each  mode. 


W'e  now  define  the  generalized  scattering  matrix  as 


ss 


<•  _ 

-''ll  ^12 

P 

.S  — 

.S '■>  1  S  t-> 

K. 

^  . 
'»•*, 
V. 


(48) 


\\  here  .S  n  is  the  power  re!letiecl  in  the  /.  direction  at  z=()  normalized  by  the  incident 
power  in  the  v  direction  at  z=0.  .S^i  ihe  transmitted  power  in  the  z  direction  at 
r  =Ar  normalized  to  the  incident  power  in  the  z  direction  at  z=0.  Because  the 
geometry  ol  the  slab  is  symmetrical  with  respect  to  z.  A'i2  =  S2i  and  S2s  =  Si,. 
I'ach  ol  the  elements  in  the  above  matrix  is  itself  a  sub-matrix.  For  example 


.Si,=  = 


■''ll.;/;; 

1 1  ;/ /;  1.'  II ;  \; /;  \/ 


(49) 


where  'I'l  TM  represents  the  'I'l  response  given  a  I'M  incident  wave.  Each  element 
in  the  sul)  m.itrix  is  a  A  IXA'  1  sub-^sub-matrix  where  Kl  is  the  number  of  modes 
taken  as  basis  functions.  The  N  |./V 2  element  of  this  sub-sub-matrix  is  given  by 

•12  J  l2 


^  /?-.//: 


0„  I  -  1  /  J  I'-n.  I  •,  1 


zds 


(50) 


/  / 


(tic  Zdi 

-T  /?-:/  /2  ‘  I 

where  l  repre.sents  the  transmitted  wave  u.sed  for  .Ssj  calculations  and  r  the 


reflected  wave  for  the  .S'n  calculations.  The  numerator  is  the  .square  root  of  the 
scattered  p<iwer  passing  through  the  unit  cell  in  the  z  direction  for  the 
.V|  =  m  l+ii  1  mode.  a,„  1 ,  1  contains  the  phase  information  for  the  scattered  field. 
The  denominator  is  the  square  root  of  the  incident  power  passing  through  the  unit 
cell  in  the  z  direction  for  the  N  ^  =  m  2+n  2  mode. 


.5.  Experimental  Methods  ( 10] 

.A  waveguide  technique  was  used  lo  obtain  the  permittivity  of  a  honeycomb 
slab.  As  show  n  in  Figure  9,  the  honeycomb  .sample  in  a  waveguide  is  equivalent  to 
a  honeycomb  panel  illuminated  by  a  plane  wave  with  an  incident  angle  of 
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(a)  Honeycomb  sample  in  \>.aveguide 


Honeycomb  Slab 


hquivaleni  honeycomb  slab  illuminaled  by  a  plane  vv'ave 
at  incident  angle  9  =  sin“‘(X/X.  ) 


I  igure  y.  Waveguide  measurement  technique 


6  =  sin 

K 


(.■>1) 


where  X  is  ihe  wavelength  of  the  incident  wave  and  X,.  is  the  cutoff  wavelength  of 
the  waveguide.  Therefore,  by  measuring  the  insertion  phase  delay  and  the  insertion 
loss  of  the  honeycomb  sample  in  the  waveguide,  we  can  determine  the  effective  per¬ 
mittivity  of  the  corresponding  honeycomb  panels. 

The  measurement  technique  is  as  follows:  .X  \VR175  waveguide  without  the 
honeycomb  sample  is  connected  to  a  network  analyzer  (HP8409)  and  the  system  is 
calibrated  over  the  frequency  band  from  10  to  15  (illz.  The  honeycomb  is  inserted 
into  the  waveguide  and  the  amplilucie  and  pha.se  of  the  parameter  is  measured. 
From  this  measurement,  an  equivalent  simple  dielectric  is  calculated  which  exhibits 
the  same  insertion  loss  and  phase  delay  (.V,,)  for  the  given  frequencies  and  incident 
angles.  Results  of  these  measurements  will  be  given  in  the  results  section  (Table  3). 

6.  Results 

\ke  will  examine  the  S  parameters  for  the  case  of  normal  incidence  on  two 
different  structured  slabs  using  the  methods  developed  above.  Case  one  is  meant  to 
demonstrate  the  additional  information  provided  by  the  full  wave  analysis  over  the 
effective  medium  theories  (in  this  case  the  small  contrast  theory).  Case  two  is 
meant  to  compare  the  two  analytic  methods  with  experimental  results  for  an  actual 
honeycomb.  The  slabs  in  both  cases  are  0.(K)lm  thick  and  the  frequency  is  10  GHz. 
The  unit  cell  lor  slab  one  is  a  square  dielectric  filling  one-quarter  of  a  square  cell. 
The  cell  side  is  0.(K)f>93m  long.  The  mi  trix  is  free  space  and  the  inclusion  is  a  loss¬ 
less  dielectric  with  e,  =  2.0.  The  unit  cell  slab  for  two  is  the  honeycomb  shown  in 
Figure  4  with  Sl=0.()()3(i7m  and  dd  =  l.()2xl0“^m.  The  matri.x  is  free  space,  but 
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ihe  inclusion  is  slightly  lossy  (€,  =  3.29—  y.0757)  Iteing  a  composite  reinforced 
with  Kevlar  49  and  bonded  with  I'pon  S2S  epoxy. 

For  the  effective  medium  approach,  the  fill  factor  is  0.25  for  slab  one  and 
0.0255  for  slab  two  The  spheroid  is  defined  by  taking  the  length  of  the  major  axis 
to  be  equal  to  the  length  of  the  unit  cell  (O.(K)lm).  The  area  of  the  inclusion  in  the 
x-y  plane  is  collapsed  into  a  circle,  and  the  diameter  of  this  circle  is  taken  to  be  the 
length  of  the  minor  axis  of  the  spheroid  Once  the  effective  epsilon  is  calculated,  the 
scattering  parameters  are  calculated  for  the  simple  dielectric  slab  0.(K)lm  thick. 

In  the  full  wave  approach,  we  must  decide  how  many  points  should  be  taken 
in  the  two-dimensional  F.F.T.  to  describe  the  unit  cell's  structure,  and  how  thick 
the  slab  can  be  so  that  a  single  pulse  describes  the  field  in  the  z  direction.  Both  of 
these  questions  are  answered  numerically.  The  number  of  points  are  determined  by 
deciding  on  how  accurate  the  scattered  fields  should  be  and  then  increasing  the 
number  of  sampling  points  until  this  accuracy  is  achieved  Since  the  approximate 
scattered  fields  will  oscillate  about  the  true  scattered  fields  with  increasing  numbers 
of  sampling  points.  \cc  must  increase  the  number  of  points  beyond  where  the  accu¬ 
racy  criterion  first  appears  to  be  satisfied  in  order  to  insure  true  accuracy.  In  unit 
cell  one.  25x25  points  were  taken  leading  to  an  accuracy  of  10"^  while  in  unit  cell 
two.  45  points  were  taken  in  the  x  direction  and  77  in  the  y  direction  to  achieve 
10“-  accuracy  (thinner  structure).  The  slab  thickness  (O.OOlm)  was  arrived  at  by 
comparing  the  program  solution  of  a  .simple  dielectric  slab  to  the  analytical  solution 
in  order  to  understand  the  relation  of  slab  thickne.ss  and  relative  dielectric  constant 
to  accuracy.  The  thickness  that  gave  10”'  accuracy  for  a  slab  composed  entirely  of 
inclusion  was  used  as  a  worst  case  estimate  lor  the  honeycomb  slab.  The  results  are 
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TE/TE 
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(-.183i:-02.+  .237P-02) 

(+.0O01-+(K),+  .00OE+00) 
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(-.294E  02.+.382P-02) 

(  +  .0O0P+00.+.0O0E+00) 
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+0.-1 
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(+.0(X)E+(X).+.(XX)E+00) 
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(-.516E-02.+.670i;-02) 

( +  .000P+00,+ .0(X)1:+(K)) 

Table  2.  .Sjj  parameter  calculated  by  the  lull  wave  analysis  compared  to  the  small 
contrast  theory  for  a  square  block  of  dielectric  filling  1/4  of  a  square  unit 
cell.  The  side  of  the  unit  cell  is  ().()0693  m.  e,  =  2.0  Only  TE/TE  and 
TM/TE  results  are  shown  since  TE/TM  and  TM/T.M  are  similar. 


Solution  Method 
Small  Contrast 
Maxwell-Garnet 
Bruggeman 
^^M^G_S£heroid_ 
^__^uB_Wave__ 
Experimental 


_ SIJ _ 

(-.148E-02.-.585E-02) 

(-.787i:-03.-.334E-02) 

(-.809E-03.-.344E-02) 

(-.660E-03.-.283E-02) 

(-■329E-02.-.125E-01) 

(-.258E-02.-.687E-02) 


_ S2J _ 

(+.977E+00.-.214n+00) 
(+  977E+()0.-.2nr.+00) 
(+.977E+00.-.212E+00) 
(-t.977E-KX).-.211E-t-00) 
(-s.975E400.-.221E-f00) 
(+.976E+00-.215E+00) 


Table  3.  Sn  and  S21  calculated  by  all  the  effective  medium  theories  and  the  full 
wave  analysis  using  only  the  0,0  mode.  These  are  compared  to  measured 
experimental  results  for  a  honeycomb  unit  cell.  S1»K).00367  m. 
dd  =  1.62xl0-^m  .  and  e,  =  3.29-/0.0757. 


shown  in  Table  1  through  ^ 


1.  Conclusion 

We  have  seen  that  of  the  effective  medium  theories,  the  small  contrast  approx¬ 
imation  provides  the  best  agreement  to  both  experiment  and  to  the  0.0  Floquet  har¬ 
monic  of  the  full  wave  .solution.  Further,  the.se  results  are  obtained  .several  orders 
of  magnitude  quicker  than  a  full  wave  solution  and  are  easily  cascaded  with  the 
F’.S.S.  parameters.  The  full  wave  solution,  on  the  other  hand,  provides  information 
on  how  higher  order  modes  behave  and  how  the  structure  changes  the  polarization 
of  the  incoming  wave.  We  found  that  .some  of  the  higher  order  modes  are  as  large 
as  the  0.0  mode  implying  that  structure  does  play  an  important  role.  Al.so  note  that 
in  Table  1.  the  scattered  fields  are  showing  two-fold  symmetry  around  the  0,0  mode 
due  to  the  two-fold  symmetry  of  the  unit  cell. 

One  of  the  major  problems  with  the  full  wave  analysis  is  the  need  to  take  so 
many  sampling  points  to  describe  the  thin-shell  structure,  although  we  do  not  need 
to  zero  pad  as  we  do  with  non-periodic  structures.  Future  work  will  concentrate  on 
increasing  the  efficiency  of  the  full  wave  analysis,  possibly  through  the  use  of 
acceleration  techniques  for  the  slowly  convergent  .series. 
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ABSTRACT 

Space  antenna  designs  have  progressed  significantly  during  the  past  decade. 
Recent  designs  combine  requirements  for  high  spatial  (sidelobe)  isolation  and  high 
polarization  isolation  for  dual  polarized  use. 

Conventional  far  field  ranges  offer  various  disadvantages.  Ground  reflections 
may  contaminate  low  sidelobe  antenna  measurements.  Cross  polarization  measurements 
reuqire  a  source  correction  when  the  cross  polarization  level  of  the  source  antenna  is 
comparable  to  or  worse  than  that  of  the  antenna  under  test.  Intermediate  ranges  (less 
than  2D^/x)  may  require  correction  for  near  field  effects.  Furthermore,  the  high  price 
of  land  and  the  considerable  expense  of  maintaining  equipment  out-of-doors  make  the 
construction  and  operaton  of  such  ranges  a  very  ambitious  endeavor. 

In  addition  to  offering  an  attractive  solution  to  these  problems,  near  field 
antenna  measurement  offers  other  advantages.  Since  a  near  field  facility  can  be  indoors, 
all-weather,  RF  shielded  operation  is  possible.  The  facility's  considerable  automation 
capability  allows  unattended  operation.  Also,  the  antenna  under  test  at  a  planar  near 
field  range  need  not  move  during  the  measurement. 

This  paper  reviews  development  by  the  Hughes  Aircraft  Company's  Space  and 
Communications  Group  (SCG)  of  a  planar  near  field  range  for  testing  satellite  antennas 
at  SCG's  El  Segundo,  California,  facility.  Features  of  the  range  are  discussed,  and 
certain  problems  encountered  in  its  operation  are  described. 


INTRODUCTION 

In  early  1982  the  Hughes  Aircraft  Company  Space  and  Communications  Group 
(SCG)  was  awarded  a  contract  for  the  Intelsat  VI  communications  satellite.  Several 
problems  were  identified  with  using  the  existing  far  field  range  facility  for  systems 
testing  on  the  Instelsat  VI  antennas  (Figure  I).  The  most  significant  problems  were: 
1)  range  reflections,  2)  near  field  effects  due  to  short  range  length,  3)  inadequacy  for 


making  cross  polarization  measurements,  4)  RF  interference,  and  5)  mechanical  support 
of  the  antenna  to  simulate  a  zero  gravity  environment.  These  problems  and  the  high 
estimate  of  cost  of  developing  an  adequate  far  field  range  at  the  SCG  El  Segundo, 
California  site  led  to  the  decision  to  construct  a  planar  near  field  range  for  satellite 
antenna  measurements. 


FACILITY 

The  near  field  test  facility  adjoins  the  satellite  integration  (high  bay)  area;  large 
double  doors  between  the  two  areas  facilitate  movement.  The  facility  is  large  enough 
(50  feet  long  by  40  feet  wide  with  a  40  foot  ceiling  height)  to  allow  testing  of  the 
Intelsat  VI  satellite  with  all  antennas  deployed. 

Although  the  facility  was  not  designed  as  an  RFl  free  enclosure,  its  walls  are 
double  lined  with  metal  sheets,  and  the  inside  is  covered  with  RF  absorbant  material. 
Problems  previously  encountered  with  the  outdoor  range  concerned  externally  generated 
signals  detected  in  the  passband  of  the  satellites  being  tested.  Tests  show  that  this 
facility  provides  an  effective  isolation  of  approximately  60  dB  in  the  frequency  bands  of 
interest. 


MECHANICAL  DESIGN 

As  mentioned  previously,  the  facility  is  designed  for  planar  near  field  scanning. 
Figure  2  shows  the  scanner  structure  in  the  assembled  state.  The  scanner  can  measure 
an  area  21  feet  on  a  side,  (horizontally  and  vertically).  Horizontal  scanning  is  provided 
by  moving  the  large  truss  structure  (carriage)  on  two  cylindrical  rails  located  at  the  top 
and  bottom  of  the  scanner.  Vertical  scanning  is  accomplished  by  moving  the  probe 
platform  (to  which  the  near  field  probe  antenna  is  attached)  along  two  smaller  rails 
mounted  to  the  carriage.  Both  sets  of  rails  have  supports  that  are  adjustable  in  three 
dimensions  for  precise  alignment.  Two  dc  servo  motors  serve  as  drives.  One  motor 
moves  the  carriage  horizontally;  the  o^’her  moves  the  probe  platform  vertically.  For 
movement  along  each  axis,  motive  force  is  translated  from  each  motor  by  a  gear  reducer 
and  a  chain  drive;  a  controller  for  each  motor  ensures  uniform  motion  and  precise 
positioning;  and  a  laser  interferometer  system  for  each  axis  provides  precise  position 
information  (Figure  3). 

MEASUREMENT  SYSTEM 

Theory  requires  that  the  near  field  amplitude  and  phase  of  the  antenna  under  test 
be  measured  over  a  planar  surface.  The  system  is  designed  to  time  multiplex 
measurements  while  switching  frequency,  polarization,  and  antenna  beam.  Up  to  six 
beams  may  be  sampled  in  any  sequence.  Many  of  the  satellite  antennas  built  by  Hughes 


SCG  feature  multiple  beam  designs.  One  of  these,  the  Intelsat  VI,  hemi/zone  antenna, 
achieves  sixfold  frequency  reuse  through  a  combination  of  spatial  and  polarization 
isolation.  Figure  4  is  a  block  diagram  of  the  measurement  system.  The  amplitude  and 
phase  measurements  are  made  by  a  network  analyzer,  and  the  analog  outputs  are 
digitized  for  computer  processing. 

The  transform  used  requires  that  the  near  field  measurement  points  be  evenly 
spaced  in  two  dimensions.  Also,  for  efficiency,  near  field  measurements  are  made  while 
the  measuring  probe  is  in  motion.  To  satisfy  these  requirements,  the  position  of  the 
probe  is  constantly  monitored  by  the  computer.  The  outputs  of  the  two  interferometers 
are  read,  and  analog-to-digital  conversion  of  the  network  analyzer  output  is  initiated,  at 
the  appropriate  time.  The  system  is  monitored  for  late  measurements  to  allow  for 
optimum  adjustment  of  the  scan  rate. 

Many  of  the  antennas  that  must  be  measured  have  requirements  for  low  side- 
lobes  and  high  polarization  purity.  Therefore,  the  measurement  system  is  required  to 
have  a  large  dynamic  range  (>50  dB).  To  achieve  this,  low  noise  amplifiers  are  placed  as 
close  to  the  output  of  the  antenna  under  test  as  is  practical. 

It  is  important  for  the  signal  level  to  remain  above  the  system  threshold. 
Figure  5  shows  the  effect  of  measuring  with  a  signal  level  that  drops  below  the  threshold 
during  the  scan.  In  this  case,  the  phase  output  of  the  network  analyzer  becomes  nearly 
constant  at  the  low  levels,  resulting  in  distortion  at  boresight  when  the  Fourier  transform 
is  applied.  To  eliminate  the  distortion,  the  scan  area  was  reduced  to  exclude  the  low 
levels.  Figure  6  shows  the  resulting  pattern. 

Another  problem  resulting  from  the  large  dynamic  range  is  RF  leakage  that 
bypasses  the  near  field  probe  and  (or)  the  antenna  under  test.  Because  of  this  problem, 
care  must  be  taken  with  those  shielding  points  at  which  radiation  may  occur  (for 
example,  at  coax  or  waveguide  interfaces). 

NEAR  FIELD  PROBES 

A  key  element  in  planar  near  field  antenna  measurements  is  the  near  field 
probe.  This  antenna,  which  is  usually  designed  for  low  or  medium  gain,  is  scanned  over 
the  measurement  plane  to  characterize  the  antenna  under  test.  The  algorithm  used 
requires  measurements  of  two  orthogonal  polarizations.  Figure  7  shows  dual  polarized 
probes  for  performing  near  field  measurements.  With  this  probe,  the  two  orthogonal 
measurements  are  made  simultaneously,  and  a  second  measurement  scan  is  not  required. 

To  correct  the  patterns  of  the  antenna  under  test  for  the  probe  pattern,  the 
probe  must  be  accurately  calibrated.  Figure  8  shows  the  results  of  correction  for  the 
probe  pattern.  Clearly,  the  effect  of  correction  on  the  cross  polarized  pattern  of  the 
antenna  under  test  is  significant.  In  this  case,  it  results  from  the  cross  polarization  ratio 
of  the  probe  being  of  the  same  order  as  that  of  the  antenna  under  test.  Near  field  cross 


polarization  measurements  in  many  cases  are  likely  to  give  more  accurate  results  than 
far  field  results  because  they  include  source  (probe)  correction,  which  typically  is  not 
done  on  a  far  field  range. 

TEST  CONTROL  AND  ANALYSIS 

An  HP  1000  minicomputer  is  used  for  automated  test  control  and  test  data 
analysis.  Test  data  is  collected  on  the  system  disk  and  is  processed  at  the  conculsion  of 
the  measurement. 

The  devices  that  are  interfaced  to  the  computer  for  automated  measurement  and 
test  control  include; 

•  a  two-axis  motor  controller 

•  a  two-axis  laser  interferometer 

•  a  frequency  synthesizer 

•  analog-to-digital  converters 

•  RF  switches 

•  programmable  attenuators 

•  power  meters 

•  a  time  code  generator 

The  pattern  measurement  software  is  designed  to  run  unattended,  and  test 
progress  can  be  monitored  remotely  via  modem.  Automatic  test  restart  of  the 
measurement  is  initiated  in  the  event  of  a  laser  fault. 

Data  analysis  consists  of  the  fast  Fourier  transform  (FFT),  inverse  FFT,  the 
probe  correction,  and  the  antenna  pattern  graphics. 

Near  field  to  far  field  data  transformation  requires  several  steps.  Before 
transformation,  the  near  field  data  is  zero  padded  (that  is,  the  near  field  amplitude  is  set 
to  zero)  to  the  required  array  size.  The  FFT  is  applied,  and  the  data  is  transformed  into 
far  field  vector  (K)  space.  The  transformed  data  is  windowed,  and  an  inverse  FFT  is 
applied.  The  data  is  again  zero  padded  to  the  original  array  size,  and  the  FFT  is  applied 
a  second  time.  This  technique  results  in  higher  resolution  in  the  region  around  which  the 
data  is  windowed  without  additional  sampling  in  the  near  field.  Figure  9  shows  a 
transform  pair  (near  field  pattern  versus  far  field  pattern). 


Application  of  the  FFT  to  the  near  field  data  results  in  a  far  field  pattern  as  it 
would  appear  (theoretically)  at  an  infinite  range  length.  This  transformation  offers  some 
advantages  over  far  field  measurements  which  often  occur  in  practice.  For  high  gain 
antennas  (such  as  the  Intelsat  VI  hemi/zone  antenna),  testing  at  a  range  length  less  than 
2D^/X,  can  result  in  noticable  pattern  distortion  due  to  phase  error.  Also,  boresight 
errors  can  occur  when  parallax  is  improperly  accounted  for  on  ranges  of  a  finite  length. 
Parallax  correction  becomes  particularly  complicated  when  several  antennas  are 
mounted  on  the  same  structure. 

In  addition  to  antenna  patterns,  the  near  field  system  developed  has  the 
capability  to  derive  satellite  system  parameters,  such  as  for  EIRP  and  flux  density.  The 
techniques  for  these  measurements  use  the  near  field  scan  data,  the  transformed  data, 
and  the  probe  input  (flux  density)  or  output  (EIRP)  power.  These  techniques  are  very 
useful  because  they  do  not  require  moving  the  satellite  to  a  far  field  antenna  range  to 
derive  the  desired  parameters. 


TEST  RESULTS 

Several  satellite  antennas  have  been  measured  with  the  near  field  test  facility. 
Figure  10  shows  the  Morelos  T-1  satellite  during  near  field  tests.  This  satellite,  which  is 
designed  to  provide  domestic  communications  to  Mexico,  features  two  reflector  antennas 
for  4  GHz,  6  GHz,  and  12  GHz  beam  coverage  and  a  planar  array  antenna  for  14  GHz 
coverage.  Note  that,  although  minimal  mechanical  support  is  required  for  near  field 
testing,  the  requirements  for  far  field  testing  can  be  significant.  (See  Figure  1.)  This  is 
because  the  antenna  under  test  remains  stationary  during  the  measurement.  Figure  11 
compares  near  field  pattern  measurements  with  far  field  measurements  for  a  typical 
Hughes  satellite  antenna. 


CONCLUSION 

Applying  the  concept  of  planar  near  field  techniques  and  subsequently  developing 
this  facility  for  testing  satellite  antennas  provides  the  Hughes  Space  and 
Communications  Group  with  an  effective  alternative  to  developing  a  larger  far  field 
range.  Several  difficult  problems  associated  with  testing  state  of  the  art  satellite 
antennas  have  been  eliminated.  The  facility  already  is  highly  automated  and  offers 
further  growth  potential  to  satisfy  future  requirements. 


w. 
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1.  ABSTRACT 

A  highly  specialized  VHF  pattern  measurement  range  has  been 
constructed  for  the  measurement  of  two  circularly  polarized 
antennas.  Unusual  measurement  problems  were  created  by  the 
satellite  size  and  the  low  VHF  frequencies  of  the  antenna  system. 
The  two  antenna  systems  operate  over  bandwidths  of  nearly  fifty 
percent  and  are  capable  of  receiving  both  right  and  left  hand 
circular  polarization.  Due  to  the  long  wavelengths  Involved, 
surface  currents  on  the  satellite  and  Its  solar  panels  are  a 
significant  factor  In  the  net  antenna  performance. 

The  maximum  dimensions  of  the  satellite  are  nearly  equal  to 
the  twenty  foot  wavelength  of  the  lowest  frequency  to  be  tested; 
therefore,  although  the  antenna  directivity  Is  low,  the  aperture 
to  be  tested  may  be  as  large  as  one  square  wavelength. 

Even  though  the  range  shown  In  Figure  1  appears  to  be  a 
conventional  horizontal  configuration,  there  are  two  specialized 
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features  that  have  significantly  reduced  measurement  errors  due  to 
ground  reflections  and  multiple  scattering  between  the  two  closely 
spaced  towers. 

2.  RANGE  CONFIGURATION 

Since  a  remote  controlled  battery  powered  source  transmitter 
Is  used  Inside  the  satellite  to  eliminate  the  possible  effects  of 
a  detector  cable  on  pattern  results,  that  which  Is  normally 
considered  the  source  antenna  Is  referred  to  as  the  receiving 
array  In  this  configuration. 

The  horizontal  range  distance  was  set  at  102  feet,  which  Is 

2 

20  /x  at  the  highest  test  frequency  for  the  maximum  aperture 
presented  by  the  satellite.  It  was  desirable  to  use  an  array  to 
place  a  pattern  null  In  the  specular  reflection  zone  of  the  range 
as  shown  In  Figure  2. 

This  array  was  Initially  a  four  (4)  element  array  consisting 
of  four  (4)  crossed  dipole  elements  fed  In  quadrature  and  backed 
by  a  reflector  with  parasitic  elements  between  the  dipoles  and  the 
reflector  for  Increased  bandwidth  (designed  by  Howard  King  and  Or. 
Jimmy  Wong  of  Aerospace  Corp.).  The  dimensions  of  this  array  are 
shown  In  Figure  3. 
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Q 

c 

□ 

E 
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REFLECTOR 

SIZE 

LOW  BAND 

8.46' 

D 

1.16' 

.45' 

3.52' 

Q| 

27.5'  X  27.5' 

The  array  spacing  was  set  to  a  dimension  that  maintained 
the  energy  in  this  zone  at  least  25  dB  below  the  pattern  maximum 

throughodt  the  50  percent  band  to  be  tested.  This  array  was 

tested  and  as  portrayed  in  Figure  4,  performed  as  expected. 

From  an  analysis  by  Dr.  C.  J.  Chen,  of  Rockwell  Inter¬ 
national,  multiple  scattering  (depicted  in  Figure  5)  between  this 
receiving  array  and  the  satellite  can  introduce  significant 
measurement  errors.  The  worst  case,  for  example,  is  when  the 
satellite  solar  panels  are  at  0  degrees,  which  presents  the 

largest  plane  reflecting  surface  area  normal  to  the  range 
boresight  axis. 

The  multiple  scattering  is  then  magnified  if  the 
receiving  array  antenna  is  backed  by  a  reflecting  screen. 
Possible  interference  to  signal  ratios  for  such  conditions  are 

shown  in  Figure  6.  Since  the  measurement  of  axial  ratio  is  the 
most  sensitive  to  error  caused  by  reflections,  it  is  desirable  to 
significantly  reduce  the  multiple  scattering  cross  section  of  the 
receiving  array.  Figure  7  shows  the  desired  interference  to 
signal  ratio  to  achieve  from  0.5  dB  to  1.5  dB  measurement 
accuracy.  To  reduce  the  receiving  array  scattering  aperture 
enough  to  achieve  these  levels  requires  both  eliminating  the 
reflector  screen  and  the  use  of  non-resonant  array  elements, 
approximately  0.2  x  long.  Figure  8  shows  the  scattering  aperture 
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Figure  7.  Axial  Ratio  vs.  Interference  to  Signal  Ratio 
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for  elements  from  0.1  to  0.5  a.  As  shown,  if  normalized  to  0.5  x 
resonant  element,  the  scattering  aperture  of  a  0.2  x  element  is 
24.1  dB  lower. 


Using  non-resonant  elements  does  lower  the  array  efficiency 
due  to  VSWR  mismatch,  but  transmitter  power  levels  from  -^20  dBm  to 
-^30  dBm  proved  to  be  adequate  for  the  102'  range  spacing. 

From  this  analysis,  and  the  weight  limitation  of  the  mobile 
tower,  a  non-backplaned  array  with  non-resonant  elements  was 
built,  tested  and  utilized.  This  array  is  also  made  of  four  (4) 
crossed  elements  fed  in  quadrature,  through  a  remote  controlled 
transfer  switch,  to  permit  changing  polarization  sense  during 
testing  as  shown  in  Figure  9.  One  spacing  of  this  array  was 
chosen  to  provide  the  lowest  energy  possible  in  the  specular 

reflection  region  over  the  entire  band.  Representative  patterns 
of  this  array  are  shown  in  Figures  10,  11  and  12,  and  the  energy 
in  the  specular  reflection  zone  is  graphed  in  Figure  13. 

3.  SATELLITE  TESTS 

Two  methods  are  used  to  derive  axial  ratio.  The  test 

equipment  set-up  is  shown  in  Figure  14.  In  the  first  method,  the 

satellite  is  spun  on  axis,  against  a  fixed  linear  receiving 

antenna.  This  is  equivalent  to  the  standard  spinning  linear  method 
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Figure  13.  Energy  in  specular  reflection  zone  for  non-backplaned 
shortened  element  array 


tests 


of  determining  axial  ratio.  This  method  Is  only  valid  when  the 
satellite  Is  oriented  normal  to  the  range  boresight  axis  where  the 
resulting  pattern  Is  the  equivalent  of  a  polarization  pattern. 

These  patterns  are  taken  with  the  receiving  array  phased  for 

horizontal  linear,  and  then  again  for  vertical  linear.  The  axial 
ratio  Is  determined  by  averaging  the  value  of  these  two  patterns. 

Off  axis,  however,  pattern  amplitude  variations  are 

superimposed  on  the  polarization  ellipse.  Therefore,  In  the 
second  method,  the  axial  ratios  for  the  off  axis  angles  are 

computed  from  the  Isolation  levels  between  the  right  and  left  hand 
polarization  senses.  Since  the  receiving  array  does  not  exhibit 
perfect  circular  polarization,  axial  ratio  cannot  be  directly 
computed  from  the  isolation  data;  therefore,  the  analysis  by  Dr. 
Allen  Love  of  Rockwell  International  (depicted  In  Figure  15),  Is 
used  to  bracket  the  actual  value  of  axial  ratio. 

The  set-up  for  measuring  gain  Is  shown  In  Figure  16.  For  gain 
testing,  the  array  Is  re-configured  to  provide  either  horizontal 
or  vertical  linear  response  to  eliminate  some  problems  In 
accurately  correcting  for  polarization  loss  between  two  (2) 
elHptlcally  polarized  antennas.  Gain  Is  measured  for  both 
horizontal  and  vertical  linear  states  of  the  receiving  arrey  to 
arrive  at  an  averaged  level.  Two  element  Yagi  antennas  were 
calibrated  for  use  as  a  standard  gain  reference. 

8SSE2334B 


493 


AXIAL  RATIO  IN  dB  OF  RECEIVING  ARRAY 


receiving  array  Axial  Ratio  as  parameter 


AOTOHATED  THUX-ANTENIIA  POLARIZATION  MEASDSEMENTS 
USING  DIGITAL  SIGNAL  PROCESSING 

John  R.  Jones  and  Doren  W.  Hess 
Sclentlf ic-Atlanta,  Inc. 

Abstract 

In  this  paper  we  present  a  three-antenna  measurement  procedure  which 
yields  the  polarization  of  an  unknown  antenna  to  an  accuracy  comparable  to 
that  of  the  improved  method  of  Newell.  The  complete  method  is  based  on  step- 
scan  motion  of  the  two  polarization  axes  on  which  the  antenna  pairs  are 
mounted.  As  a  special  case  this  step-scan  procedure  includes  the  usual  single 
axis  polarization  pattern  method  of  polarization  measurement. 

This  three-antenna  polarization  measurement  method  can  be  readily 
automated  and  is  carried  out  straightforwardly  with  the  assistance  of  a 
minicomputer  for  data  acquisition  and  data  reduction.  The  data  reduction 
method  is  based  on  conventional  digital  Fourier  transform  techniques  and  has 
the  advantage  of  Inherent  noise  rejection.  It  utilizes  a  large  number  of 
sample  points  which  greatly  overdetermine  the  parameters  to  be  measured. 

The  method  has  been  verified  experimentally  with  measurements  made  on 
multiple  overlapping  sets  of  three  antennas,  as  is  conventional  for  this  kind 
of  procedure.  The  data  are  presented  for  broad-beam  antennas  of  the  type  used 
as  near  field  probe  horns. 

1.  Introduction 

The  purpose  of  a  polarization  measurement  is  to  determine  the  tilt  angle, 
axial  ratio  and  sense  of  the  electric  field  polarization  ellipse  by  means  of  a 
procedure  carried  out  on  an  antenna  measurement  range.  In  general,  these 
parameters  are  functions  of  direction,  and  therefore  in  general  a  polarization 
measurement  of  the  field  of  an  antenna  must  be  made  at  each  pattern  direction 
of  interest.  Such  solid  angle  measurements  of  polarization  require  knowledge 


of  the  polarization  of  the  range  antenna  with  which  the  measurement  Is  made. 

When  the  range  antenna  Itself  has  unknown  polarization  characteristics,  a 
three-antenna  polarization  measurement  Is  performed  to  determine  the 
polarization  characteristics  of  the  range  antenna.  Because  of  the  complexity 
of  the  data  reduction  associated  with  a  three— antenna  polarization  measurement 
procedure,  such  a  procedure  Is  usually  limited  to  determining  the  polarization 
of  an  antenna  In  only  one  direction  of  the  radiation  pattern.  By  performing 
polarization  pattern  measurements  on  each  of  the  three  possible  pairs  of 
antennas  In  a  three  antenna  set  one  Is  able  to  ascertain  the  axial  ratio,  tilt 
angle,  and  sense  of  the  polarization  ellipse  for  each  of  the  three  antennas. 
Usually,  the  polarization  for  the  direction  of  the  peak  of  the  beam  for  each 
of  the  three  antennas  Is  measured.  The  only  a— priori  knowledge  required  Is 
the  approximate  tilt  angle  (within  45*)  of  one  of  the  antennas  In  the  three 
antenna  set  and  the  fact  that  one  of  the  antennas  Is  reciprocal. 

A  number  of  different  procedures  have  been  devised  In  the  past  for 
performing  three-antenna  measurements  of  the  complex  polarization  ratios  of 
antennas.  Newell  and  iCems^^^  and  Joy  and  Paris Independently  developed 
three-antenna  phase/amplltude  polarization  measurement  methods  In  the  early 
1970's.  Newell,  Baird,  and  Wacker^^^  described  the  NBS  extrapolation 
technique  for  the  measurement  of  polarization  and  gain  In  1973  and  gave  some 
experimental  results  obtained  using  the  method.  All  polarization  measurement 
methods  up  to  this  time  had  a  requirement  for  accurate  measurement  of  phase. 

In  fact.  In  certain  unusual  measurement  situations,  all  the  polarization 
Information  derives  from  the  phase  measurement. 

In  1975^^1  Newell  recognized  that  a  three-antenna  polarization 
measurement  method  could  be  devised  such  that  almost  all  of  the  Information 
from  the  measurement  derives  from  the  amplitude  measurement.  Newell's  method, 
as  It  has  come  to  be  called,  requires  only  that  the  sense  (either  positive  or 
negative)  of  the  180*  phase  change  associated  with  the  null  In  the 
polarization  pattern  be  known.  This  Insensitivity  to  the  measurement  of  phase 
Is  the  primary  advantage  of  Newell's  method,  and  is  responsible  for  the 
Improved  accuracy  of  the  method  compared  to  those  described  in  references 
The  method  of  polarization  measurement  described  In  this  paper,  like 


Newell's  method,  depends  on  the  fact  that  the  form  of  the  functional 
dependence  of  the  measured  quantities  on  the  desired  unknown  polarization 
parameters  Is  known.  Instead  of  assuming  phase  Information,  as  Is  done  with 
Newell's  method,  this  new  method  uses  the  measured  phase  Information  but 
selectively  removes  measurement  errors  by  digital  signal  processing,  thereby 
Improving  the  accuracy  of  the  measurement. 

The  selective  removal  of  errors  Is  accomplished  by  a  one-  or  two- 
dimensional  discrete  Fourier  transform  (DFT)  which  filters  the  measured  data 
as  a  function  of  angle  In  the  spatial  frequency  domain  and  leaves  only  the 
spatial  frequency  components  that  vary  with  angle  as  polarization  patterns 
Ideally  should.  A  method  of  polarization  measurement  using  a  one-dlmenslonal 
discrete  Fourier  transform  was  Implemented  In  1981  by  workers  at  the  Technical 
University  of  Denmark^^^,  but  the  method  described  herein  based  on  a  two- 
dimensional  discrete  Fourier  transform  has  not  been  described  before  In  the 
literature.  The  mathematical  foundation  for  the  method  Is  described  briefly 
In  the  next  section. 

II.  Msthautlcml  Basis 

The  three  antenna  method  Is  carried  out  by  measurement  of  the 
polarization  patterns  between  successive  pairs  of  antennas.  The  measurement 
configuration  has  one  antenna  transmitting  an  Illuminating  wavefront  that  Is 
received  by  the  second  antenna.  For  a  given  pair  the  relative  received  phasor 
or  normalized  voltage  response  Is  given  by  the  matrix  product 

V  -  (W^  Wj  -  (1) 

where  the  superscript  ^  denotes  the  complex  conjugate  of  the  transpose  of  a 
matrix  and  where  the  W,  and  W,.  are  the  column  matrices 


Figure  1.  Poincare  sphere  with  polar  angle  2y  and  longitude  6  shown. 

Each  point  on  the  sphere  is  associated  with  a  particular 
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(5) 


where  the  time  dependence  has  been  assumed.  The  angles  41  and  x  are  the 

angles  through  which  the  transmitting  antenna  and  receiving  antenna 
respectively  are  rotated  in  the  clockwise  direction  as  each  is  viewed  from  the 
other.  Figure  2  Illustrates  the  coordinate  systems  for  the  polarization 
measurements  and  defines  the  rotation  angles. 

Equation  (5)  will  be  recognized  as  a  simple  two-dimensional  Fourier 
series  In  the  angles  ^  and  x*  Only  two  spatial  frequency  components  In  each 
variable  are  present — the  positive  and  negative  single  cycle  per  cycle 
harmonics. 

* 

*  • 

The  coefficients  of  the  Fourier  series  can  be  computed  from  the  known 
parametric  dependence  of  V  on  ^  and  x  t}y  Inversion: 
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Taking  the  ratio  of  these  two  quantities  we  get  that 
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where  p  is  the  circular  polarization  ratio  for  the  transmitting  antenna 
»  ct 

and  is  the  receiving  circular  polarization  ratio  for  the  receiving 
antenna. 


Using  the  relation  for  a  reciprocal  antenna  between  the  receiving 
circular  polarization  ratio  and  the  transmitting  circular  polarization  ratio 
that 
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cr 


ct 
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we  see  that  this  result  (8)  becomes  simply 


c  c 


(10) 
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where  p  is  written  a  and  o  has  replaced  o 

ct  ^c  ^c  ^cr 


Now,  both  circular 

polarization  ratios  for  the  two  antennas  refer  to  transmission  and  the  prime 
differentiates  the  two  antennas. 

Once  we  have  the  products  of  the  complex  circular  polarization  ratios  for 
each  of  the  three  antenna  pairs  determined  from  successive  polarization 
measurements,  we  can  compute  each  of  the  circular  polarization  ratios  from 
three  equations  In  three  complex  unknowns  derived  as  follows.  For  simplicity 
of  notation  we  define  the  products  of  the  complex  circular  polarization  ratios 
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where  the  quantities  and  Mg^  are  complex  numbers  derived  from  the 

data  reduction. 

The  solutions  for  the  circular  polarization  ratios  are  given  by 
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^AC 

^AB 

(12) 


Resolution  of  the  sign  ambiguity  Inherent  In  taking  the  square  root  follows 
from  an  a-prlorl  knowledge  of  the  approximate  tilt  angle  of  one  of  the  three 
antennas. 

Once  the  complex  circular  polarization  ratio  for  an  antenna  Is  known,  the 
axial  ratio,  tilt  angle,  and  sense  of  polarization  may  be  computed  from  It. 

The  details  of  this  computation  are  given  In  Hollis,  Lyon,  and  Clayton. 


III.  Bzperlaeatal  Procedure 

The  experimental  procedure  for  any  three’-antenna  polarization  measurement 
method  Is  to  measure  the  relative  received  phasor  voltage  for  the  three 
antennas  taken  pairwise.  The  voltage  as  a  function  of  rotation  about  the  line 
of  sight  Is  recorded.  This  rotation  is  described  by  the  angles  ^  and  x  as 
shown  in  Figure  2.  The  measured  data  is  reduced  using  expressions  for  the 


iV 


polarization  parameters  in  terms  of  the  relative  received  phasor  voltage  such 
as  those  given  in  Section  II. 

Initial  alignment  of  the  antenna  coordinate  system  is  Important  with  any 
three-antenna  polarization  measurement  method  because  of  the  need  to  maintain 
a  constant  orientation  reference  for  measurement  of  the  tilt  angle.  The 
coordinate  systems  of  the  antennas  in  the  measurement  must  be  defined  by 
mechanical  reference  devices  such  as  bubble  levels  mounted  on  the  antennas. 

The  measurement  procedure  described  herein  is  used  to  measure  directly  the 
polarization  of  antennas  in  the  coordinate  systems  so  defined.  No  coordinate 
system  redefinition  is  required.  In  practice,  this  means  that  the  antennas 
are  mounted  with  their  mechanical  reference  devices  in  their  reference 
orientation  (e.g.  with  their  bubble  levels  level).  The  angular  offsets 
in  1^  and  x  corresponding  to  this  condition  are  noted  and  these  offsets  used  to 
define  the  zero  rotation  angles  for  the  coordinate  systems. 

The  three-antenna  polarization  measurement  method  described  herein  is 
based  on  step-scan  motion  of  the  two  polarization  axes  on  which  the  antenna 
pairs  are  mounted.  The  complex  voltage  (phase  and  amplitude)  at  the  port  of  a 
receiving  antenna  is  measured  by  a  phase/amplitude  receiver  and  digitally 
recorded  as  one  of  the  antennas  is  scanned.  This  scanning  consists  of 
rotating  the  source  antenna  through  an  angle  ^  *  360*  or  the  receiving  antenna 
through  an  angle  x  “  360®. 

For  single  scan  polarization  measurements,  which  are  reduced  to 
polarization  parameters  by  means  of  a  one-dimensional  DFT,  this  is  the  only 
measurement  that  is  required  for  each  antenna  pair.  The  sample  increment  has 
only  a  secondary  influence  on  the  accuracy  of  the  measurement.  Our 
experimental  results  indicate  that  no  further  improvements  in  measurement 
accuracy  result  from  use  of  sample  increments  less  than  I*. 

For  multiple  scan  polarization  measurements,  the  axes  are  rotated  in  a 
step-scan  fashion.  Following  acquisition  of  single  scan  data  as  described 
above,  the  polarization  axis  which  was  not  used  for  the  scan  (logically,  this 
axis  may  be  referred  to  as  the  step  axis)  is  rotated  through  a  discrete 
angular  increment  and  the  scan  is  repeated.  The  phase  and  amplitude  data  are 


recorded  as  functions  of  both  the  step  and  scan  axis  angles  (i^  and  x)  and  a 
two-dimensional  DFT  Is  used  for  reduction  of  the  measured  coupling  data  to 
polarization  parameters.  As  In  the  single  scan  measurements  neither  the  step 
nor  the  scan  angle  increment  has  a  profound  Influence  on  the  accuracy  of  the 
measurement.  Our  experimental  results  Indicate  that  no  further  Improvements 
In  measurement  accuracy  result  from  use  of  step  angle  Increments  less  than 
15". 

IV.  Bxperiaeotal  Sesnlts 

A  set  of  three  single-ported  antennas  and  one  dual-ported,  dual-polarized 
antenna  was  used  to  verify  the  measurement  method.  All  possible  combinations 
of  the  five  antennas  (the  dual-ported  antenna  Is  considered  to  be  two  antennas 
and  Its  two  ports  are  referred  to  In  the  tables  as  antennas  number  4  and 

number  5)  were  measured,  with  the  results  shown  In  Tables  I  and  II. 

For  the  single  scan  measurements  shown  In  Table  I  the  total  variation  In 
the  measured  axial  ratio  appears  to  be  very  large.  However,  because  the  axial 
ratio  approaches  Infinity  for  nearly  linear  antennas,  total  variation  in  the 
axial  ratio  is  a  deceptive  measure  of  accuracy.  Note  that  the  total  variation 
in  polar  angle  ly  for  all  the  antennas  is  0.29®  -  0.44*.  This  is  a  true 
indication  of  the  size  of  the  locus  of  uncertainty  on  the  PoincarS  sphere, 
while  axial  ratio  variation  is  not.  Note  in  particular  that  for  antenna  #2, 
with  an  axial  ratio  of  28.44  dB,  a  0.78  dB  total  variation  corresponds  to  a 
larger  locus  of  uncertainty  on  the  Poincar€  sphere  than  does  the  5.67  dB 
variation  for  antenna  number  3  with  an  axial  ratio  of  49.31  dB,  the  9.44  dB 

variation  for  antenna  number  4  with  an  axial  ratio  of  57.50  dB  or  the  2.74  dB 

variation  for  antenna  number  5,  with  an  axial  ratio  of  38.92  dB. 

The  single  scan  measurements  shown  in  Table  I  indicate  the  degree  of 
accuracy  achievable  using  a  one-dimensional  discrete  Fourier  transform  for 
reduction  of  the  polarization  pattern  data.  Such  a  DFT  tends  to  reject  random 
phase  and  amplitude  errors  in  the  measurement  such  as  stray  reflections  in  the 
test  volume.  However,  periodic  errors,  such  as  rotary  joint  phase  and  VSWR 
wow,  are  not  rejected  by  the  one-dimensional  DFT  since  they  have  the  same 


period  as  the  polarization  pattern  (2ir).  But  these  periodic  errors  (in 
particular,  rotary  joint  errors)  are  dependent  only  on  the  rotation  angle  of 
one  polarization  positioner  (the  one  in  which  they  are  mounted),  while  the 
polarization  pattern  is  dependent  on  the  rotation  angle  of  both  polarization 
positioners.  Therefore,  the  periodic  errors  in  the  measurement  do  not  vary  in 
the  same  way  as  the  polarization  pattern  does  as  a  function  of  both 
measurement  angles,  ^  and  x»  and  they  will  tend  to  be  rejected  by  the  two- 
dimensional  DFT  in  i)i  and  x  which  is  used  to  reduce  multiple  scan  polarization 
measurements. 

Table  II  shows  the  results  of  multiple  scan  measurements  on  the  same 
antennas.  For  multiple  scan  measurements,  periodic  as  well  as  random  errors 
tend  to  be  filtered  out,  resulting  in  a  2  or  3  to  1  improvement  in  measurement 
accuracy  over  single  scan  measurements.  Work  on  computer-simulated  data 
indicates  that  multiple-scan  data  reduction  (a  two-dimensional  instead  of  a 
one-dimensional  DFT)  reduces  periodic  errors  by  about  20  dB. 

V.  Cottcluslons 

We  have  demonstrated  an  automated,  highly  accurate  method  of  polarization 
measurement.  The  method  is  restricted  to  use  under  far-fleld  conditions.  The 
method  is  implemented  as  part  of  the  spherical  near-field  subsystem  of  the 
Sclentlf ic-Atlanta  2022B  antenna  analyzer  option  08A.  The  measurements 
required  are  simply  phase/amplitude  polarization  patterns  for  multiple  source 
orientations.  No  redefinition  of  coordinate  systems  (as  in  Newell's  method) 
is  required,  and  polarization  parameters  are  computed  automatically  and 
displayed  in  the  format  shown  in  Figure  3.  The  automation  makes  this  method 
easier  to  use  than  manual  methods.  This  method  has  the  further  advantage  that 
no  longitudinal  translation  is  needed  to  average  out  stray  signals.  Stray 
signals  are  rejected  by  the  step-scan  rotation. 


POLARIZATION  PARAMETERS  PROM  PILE  DP1MS3 


PREOUENCY  AXIAL  RATIO  SENSE  TILT  ANGLE 

(MHx)  (dB)  (DEGREES) 

\30gB.B  S3. 22  RHE  88.98 


ALTERNATIVE  POLARIZATION  REPRESENTATIONS 
BASIS  COMPLEX  POLARIZATION  RATIOS  POINCARE  SPHERE  COORDINATES 


AMPLITUDE 

PHASE 

(DEGREES) 

POLAR  ANGLE 
(DEGREES) 

LONGITUDE 

(DEGREES) 

CIRCULAR 

1 .0044 

177.9580 

90.25 

177.96 

LINEAR 

S5.6969 

-6.9888 

177.94 

-6.99 

DIAGONAL 

.9650 

-.2504 

CD 

-.25 

Figure  3.  Sample  polarization  parameters  output  list 
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ABSTRACT 


Simple  analytic  expressions  that  locate  the  phase  centers  of 
small  TE^q  mode  rectangular  and  conical  horns  are  presented. 

These  expressions  relate  the  phase  center  displacement  (from  the 
horn  aperture  center)  to  the  E  plane  and  H-plane  horn  dimensions 
and  to  the  horn  flare  lengths.  The  equations  compare  well  both 
with  more  elaborate  derivations  and  also  to  measurements  made  on 
several  different  sized  horns.  Because  these  expressions  are 
quite  accurate  and  not  mathematically  cumbersome,  they  are  useable 
In  many  practical  situations,  replacing  commonly  used  measurement 
techniques. 


Introduction 


Expressions  for  location  of  phase  centers  of  small  horn 
antennas  should  be  in  simple  form  to  be  useful  to  antenna  designers 
The  most  commonly  used  way  for  determining  the  phase  center 
location  is  by  measurement.  Determining  the  phase  center  location 
by  analysis  is  usually  very  complex  and  time  consuming. 

In  this  paper,  the  location  of  the  phase  center  in  the  E  and  H 
planes  is  obtained  in  a  simple  fairly  accurate  form  by  a  series 
of  approximations. 

The  approach  described  here  relates  the  aperture  phase  varia¬ 
tion  of  the  horn  to  a  spherical  wave  front  radiated  from  the  actual 
phase  center.  More  accurately  stated,  the  phase  center  is  deter¬ 
mined  by  matching  the  far  field  main  beam  phase  variation,  in  a 
best  fit  match,  to  a  spherical  wave  front.  The  center  of  that 
spherical  wave  is  the  phase  center. 

Rather  than  comparing  the  entire  phase  variation  over  the 
main  beam,  only  the  phase  at  the  angles  where  the  power  is  down 
by  3-dB  together  with  the  phase  at  boresight  are  used.  This  is 
a  three  point  approximation  to  the  spherical  wavefront.  The 
3-dB  angles  are  chosen  both  for  simplicity  and  because  they 
probably  match  the  angles  most  experimentalists  would  choose. 

Other  reasonable  angles  should  also  give  accurate  results. 

Phase  Center  Derivation 

Assume  that  the  field  pattern  phase  at  the  3-dB  angles  is 
This  phase  minus  the  boresight  phase,  gives  the  horn  phase 


deviation  with  respect  to  the  horn  aperture  center.  It  can  be 


expressed  as  a  physical  length,  e, 

2iie 

1  “  <{10  “  (1) 

X 

Figure  1  shows  the  geometry  of  the  problem, 

with  Rq  the  distance  from  the  horn  aperture  to  an  arbitrary  point  In 
the  far  field;  AL  the  displacement  of  the  phase  center  from  the 
horn  aperture;  a,  the  radius  of  a  spherical  wavefront  passing  through 
the  arbitrary  point  and  R  the  distance  from  the  horn  aperture  center 
to  the  spherical  wavefront  at  an  angle  8  from  the  boreslght  direction 
Let  83  refer  to  the  3  dB  angle  (283  is  the  horn  beamwldth. ) 

From  the  law  of  cosines, 

a2  -  (  AL)2  +  r2  -  2RAL  cos  (180-83).  (2) 

also 

a  -  Rq  +  AL  (3) 

and 

e  =  R-Rq.  (4) 

If  Rq  >  >  e,  then 


then 


e  =  AL  (I-CO893) 


(5) 


(Thp  condition  63  <  1  is  almost  met  by  any  horn  aperture  larger  than 
a  waveguide  opening).  The  phases  of  the  horn  at  boresight,  <))q, 
and  at  63,  41  p  are  easily  derived. 

In  Figure  2,  p  is  the  E-plane  flare  length.  The  two  flare  lengths 
(Eplane  and  H  plane)  will  be  slightly  different  in  general,  unless  the 
waveguide  heighc  and  width  are  proportional  to  a  and  b. 

Assume  that  the  far  field  pattern  can  be  expressed  independently 
in  the  E  and  H  planes.  Thus  ggCw),  the  far  field  pattern  for 
the  E  plane,  assuming  a  quadratic  phase  error  across  the  aperture, 
represented  by(6) 

ggCu)  “  a/2[ggQ(w)  +  j3gp^Q  (Vi)l  (7) 


where  6  is  the  phase  error  at  the  edge  of  the  horn  aperture  and  ggQ 
is  the  unerrored  pattern  of  a  uniform  aperture  of  widt''  =*  and  y  is 
is  TiaX  sin  6.  From  equation  (7)  the  phase  of  ggCw)  at  y  *  0 
(boresight)  is: 


The  phase  at  the  3  dB  point  is: 
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Since  Q  is  the  phase  difference  at  the  edge  of  the  horn  relative  to 
the  center, 


2it 


♦ 


(K 


where  =  Is  the  physical  length  from  the  flare  center  to  the  aper¬ 
ture  center,  Ag  is  the  length  from  the  flare  center  to  the  horn  edge 
In  the  E-plane  and  Ig  is  the  guide  wavelength  at  the  horn  aperture. 
Thus , 

“  /p^+(a/2)2  -  p  (i; 

a2 

8p  (1( 


Assuming  that  the  horn  phase  velocity  is  proportional  to  the  arc 
length  b  in  the  H  plane  at  the  aperture, 

A 


'g 


[l-(Ag/2b*)j 


1/2 


(!' 


2  sin  ^/2 


where 


b  -  b* 


(21 


and  is  the  total  horn  flare  angle.  Thus 

i|i2  ^2 

b  »  b'[l - -  ...  ]  =  b'tl- - J  .  (2 

24  24 


Expanding  the  expression  for  A^, 


Xg  2  X  [1  + 


(23) 


t 

t- 

1;^ 

I 


I 

k 

. 


h* 


X2  3  X2  x2  AiJ)2 

8b2  128  b2  ^b2  9 


If  i);  «  3/2  X/b  . 

Since  b  should  range  between  IX  to  6X  for  small  horns  and  i|> 
between  20*^  and  80® 

X2 

[1+— J 

8b2 

For  very  small  flare  angles  (4;  =  0),  a  better  approximation  is 

X« 


(24) 


so  that 


AL 

X 


®  /l-(X2b)2 

.135(a/X)^ 


(25) 


(26) 


(p/X  J(l+l/8(X/b)2 


For  very  small  flare  angles  and  b  1,  a  better  approximation  for  Xg/2 
gives 


X  2 

AL  ,135(a/X)^  1  -  2b 
X  p/X 

In  the  H'plane,  the  same  approach  is  used  only  the  basic  horn  illu¬ 
mination  is  now  a  cosine  taper.  This  gives  the  errorless  pattern 

cos  u 


(27) 


4p2 


(28) 
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slightly  different  than  the  p  for  the  E-plane,  and  for  small  flare 


angle  horns  with  small  be  dimension,  the  factor  1  +  —  should  be 

8b 

1/2 

replaced  by  1  +  — 

2b 

The  phase  centers  for  conical  horns  in  both  planes  can  also  be 
found  in  the  same  manner.  In  this  case  the  E-fleld,  is 


given  by:  (7) 


Ggo(h)  ct 


Ji(m) 


The  three  dB  angle  is 


63  -  .521  X/d 


where  d  is  the  conical  horn  diameter.  These  give 


AL  .0934(d/X)^ 


where  the  factor  1  +  —  (X/b)2  has  been  ignored.  It  could  be 

8 

replaced  by  1  +  (X/b)^(.29).  For  the  H  plane, 


given  by  (8) 


(u)  Ji'(  u) 

«H0  «  - 


wliere  p  is  the  first  root  of  Jj  (y).  For  this  pattern,  the  three 
dB  angle 


V-U- •  ^  **•  V>_ ■•‘••A  , 
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This  expression  again  assumes  X/g  ■  X. 

With  equations  (26)  and  (32),  the  E  and  H-plane  phase  centers  of 
a  pyramidal  horn,  and  with  equations  (35)  and  (38),  the  E  and  H  plane 
phase  centers  of  a  conical  horn  are  described. 

Comparisons  and  Experimental  Results 

Equations  26  and  32  are  used  in  Figure  4  to  compare  these  results 
with  square  horn  theoretical  results  of  Mueldorf(^).  The  results  of 

for  sectoral  horns  are  plotted  in  Figure  5  with  the  phase  centers 
determined  from  equations  27  and  32.  These  figures  show  quite  good 
agreement  between  our  approximate  results  and  a  more  rigorous 
derivation. 

Table  I  gives  the  comparison  of  the  experimental  measurements  of 
the  two  phase  centers  of  five  pyramidal  and  Table  2  a  conical  horn 
with  the  phase  centers  derived  from  the  appropriate  equation. 


SOLID  LINE  FROM  MUEHLDORF(4) 


X 

Figure  3.  Square  Horns 


SOLID  LINE  FROM  HuO) 


Figure  4.  Sectoral  Horns 


TABLE  I:  Experimental  Results 


Pyramidal  Horns 


a 

b 

pe 

ALg  meas 

ALg 

AL^  meas 

ALh 

1.84 

1 

1 

3.47 

4.84 

+.12 

.13 

.15 

+  .12 

.88 

.78 

4.68  j 

5.97 

7.56 

3.88 

3.93 

4.13 

4.37 

3.25 

8.38 

18.84 

.50 

.37 

8.44 

8.11 

2.13 

2.88 

5.39 

.31 

.23 

.19 

.30 

2.00 

3.44 

10.17 

.25 

.10 

.50 

_ i 

.36 

TABLE:  Experimental  Results 
Conical  Horns 

d  p  ^  measured  ^  ^  o^^sured  A 


8e00 

12.40 

14.81 

U.85 

7.13 

7.47* 

1 

4.00 

i  8.72 

1.88 

1  1.76 

.75 

.65 

*Corrected  for  inadequate  field  distance 

In  general  the  measured  results  compare  quite  well. 

Conclusions 

Four  simple  but  reasonably  accurate  expressions  have  been  derived 
for  the  phase  centers  of  pyramidal  and  conical  horns.  Control  of 
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antenna  astigmatism  in  reflectors  that  use  horns  as  a  primary  feed  is 
now  simple.  For  example,  any  rectangular  horn  with  the  H  field  dimen¬ 
sion,  b,  set  equal  to  1.22  times  the  E  field  dimension,  a,  will  give 
coincident  phase  centers.  As  a  second  example,  a  limit  of  1/16 
separation  between  phase  centers  and  an  equal  illumination  of  10  dB  in 
both  planes,  (this  sets  b  =  1.38a)  give  an  expression  for  the  smallest 
horn  length,  p,  (p/1  ■  .7(a/l)^)  with  the  dimension  "a"  set  by 

the  f/D  ratio. 

All  the  equations  have  been  derived  for  the  observation  point  at 
infinity,  the  use  of  a  horn  with  a  reflector  or  lens  that  is  at  a 
finite  and  possibly  close  distance  might  change  the  effective  phase 
center.  The  factor  to  account  for  this  effect  is  easily  derived.  The 
effective  phase  is 

AL  pig 

1  1  R1 

where  R  is  the  distance  from  the  horn  to  the  reflector  or  lens. 


References 


1.  Y.Y.  Hu,  "A  Method  of  Determining  Phase  Centers  and  Its  Applica¬ 
tions  to  Electromagnetic  Horns,"  Journal  of  the  Franklin  Institute 
vol.  271,  pp.  31-39,  Jan,  1961. 

2.  E.R.  Nagelberg,  "Fresnel  Region  Phase  Centers  of  Circular  Aperture 
Antennas,"  IEEE  Trans.  Antennas  Propag.,  vol.  AP-13,  No.  3,  pp. 
479-480,  May  1965. 

3.  1.  Ohtera  and  H.  Ujlie,  "Nomographs  for  Phase  Centers  of  Conical 
Corrugated  and  TE ^ »  Mode  Horns,"  IEEE  Trans.  Antennas  Propag., 
vol.  AP-23,  No.  6,  pp.  858-359,  Nov.  1975. 

4.  E.l.  Mueldorf,  "The  Phase  Center  of  Horn  Antennas,"  IEEE  Trans. 
Antennas  Propag.,  vol.  AP-13,  pp.  753-760,  Nov.  1970. 

5.  W.L.  Bowen,  L.J.  Chu,  "The  Theory  of  the  Electromagnetic  Horn," 
Proc.  IRE,  vol.  27,  pp.  51-64,  Jan.  1939. 

6.  S.  Silver,  Microwave  Antenna  Theory  and  Design,  Lexington,  MA, 
Boston  Technical  Publishers,  1964. 

7.  L.J.  Chu,  "Calculation  of  the  Radiation  Properties  of  Hollow 
Pipes  and  Horns,"  J.  Appl.  Phys, ,  vol.  11,  pp.  603-610,  Sept. 

1940. 

8.  C.W.  Horton,  "On  the  Theory  of  the  Radiation  Patterns  of  Electro¬ 
magnetic  Horns  of  Moderate  Flare  Angles,"  Proc.  IRE,  vol.  37, 

pp.  744-749,  July,  1949. 


SHAPED-BEAM  EARTH  COVERAGE  HORN  ANTENNAS*  t 


by 

Joseph  C.  Lee,  Walter  Rotman 

Massachusetts  Institute  of  Technology,  Lincoln  Laboratory 
Lexington,  Massachusetts 


Abstract 

Earth  coverage  antennas  are  used  on  communication  satellites 
to  cover  the  whole  visible  portion  of  the  earth.  Due  to  the  earth 
curvature  and  satellite  orbit  geometry,  there  are  path-length 
differences  within  the  covered  area.  To  compensate  for  this 
path  length  and  any  atmospheric  attenuation  differences,  more 
radiated  power  should  be  directed  toward  the  edge,  than  toward  the 
center,  of  Che  covered  area.  Conventional  pencil-beam  antennas 
with  unshaped  radiation  patterns  are  inefficient  for  this 
application  in  the  sense  that  more  power  is  radiated  toward  the 
nadir  where  the  path  length  is  least  and  that  less  power 
(typically  A  dB  down  from  peak)  is  radiated  toward  the  edge  of 
coverage  where  the  path  length  is  longest.  To  overcome  this 
difficulty,  antenna  beam-shaping  techniques  are  needed.  This 
paper  will  present  two  practical  shaped-beam  earth-coverage  horn 
antennas  for  a  synchronous  orbit  satellite.  One  is  a  double-flare 

*Thls  work  has  been  sponsored  by  the  Department  of  the  Air  Force. 

f’The  views  expressed  are  those  of  the  authors  and  do  not  reflect 
Che  official  policy  or  position  of  the  U.  S.  Government." 
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smooth  horn  utilizing  dual  modes  for  beam  symmetry  and  aperture 
phase  differences  for  beam  shaping.  The  other  is  an  angle-limited 
corrugated  horn  with  a  beam-shaping  zoned  lens.  Preliminary  data 
of  Q-band  breadboard  models  for  both  types  of  antennas  showed 
improvements  of  at  least  1.0  dB  In  minimum  edge-of-ear th  gain, 
compared  to  an  optimized  single  mode  corrugated  horn.  Design 
considerations  with  measured  patterns  and  gains  for  the  two 
antennas  are  included. 

I.  INTRODUCTION 

Earth  coverage  antennas  are  used  on  communication  satellites 
to  cover  the  whole  visible  portion  of  the  earth.  The  portion  of 
the  earth  visible  to  a  satellite  is  a  function  of  the  satellite 
altitude.  It  covers  less  area  at  low  than  at  high  altitude.  At 
geosynchronous  orbit  (Fig.  1)  of  22,300  miles,  the  visible  area  of 
the  earth  surface  covers  up  to  81.44®  in  latitude  for  the  north 
and  south  hemispheres,  an  area  very  close  to  the  limit  of  a  whole 
hemisphere.  When  the  antenna  is  boresighted  at  nadir  on  the  earth 
equator,  the  cone  angle  at  the  edge  of  the  earth  is  8.6®  from  the 
axis.  To  provide  some  margin  for  satellite  pointing  uncertainty, 
9®  will  be  used  as  the  cone  half  angle  for  earth  coverage  antenna 
evaluation.  Due  to  the  earth  curvature  and  satellite  orbit 
geometry,  there  are  path  length  differences  within  the  covered 
area.  To  compensate  for  this  path  length  and  any  atmospheric 
attenuation  differences,  more  radiated  power  should  be 
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direcr.ed  toward  the  edge  than  toward  the  center  of  the  covered 


An  ideal  earth  coverage  antenna  at  synchronous  orbit  should 


have:  1)  circular  symmetric  patterns,  2)  zero  spillover  beyond 


0=9°,  and  3)  an  appropriate  inverse  taper  between  the  central 


axis  and  edge  of  the  covered  earth  disc  to  compensate  for  the  path 


length  and  atmospheric  attenuation  differences.  The  path  length 


attenuation  difference  between  the  nadir  and  the  edge  of  the  earth 


is  1.32  dtt  and  the  atmospheric  loss  difference  is  about  l.A  dB  in 


Q-band  (•=  4U  GHz). 


With  a  finite-size  antenna,  it  is  impossible  to  attain  this 


ideal  beamshape.  But  it  can  be  approached  with  increasing 


aperture  size, 


II.  ANALYSIS 


Conventional  pencil-beam  antennas  with  unshaped  radiation 


pattern  are  inefficient  for  earth  coverage  application  in  the 


sense  that  more  power  is  radiated  toward  the  nadir  where  the  path 


length  is  least  and  that  less  power  (typically  4  dB  down  from 


peak)  is  radiated  toward  the  edge  of  coverage  where  the 


path-length  is  longest.  To  overcome  this  difficulty,  antenna 


beam-shaping  techniques  are  needed. 


Several  synthesis  techniques  for  generating  the  desired 


far-field  patterns  with  various  amplitude  and/or  phase 


distributions  over  a  finite  aperture  have  been  theoretically 
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evaluated  in  a  recent  paper Here  we  will  first  start  with  some 
basic  analyses  and  then  go  directly  to  present  two  practical 
shaped  beam  earth-coverage  horn  antenna  designs. 

The  gain  of  a  uniformly  Illuminated  sector  beam  radiation 
pattern  within  a  semi-flare  cone  angle  9g  and  with  no  side  lobes 
is  given  by: 

G  =  10  log 

For  0g  =  9®,  G  =  22.11  dBi.  This  value  can  be  used  as  an  upper 
bound  design  goal  for  geosynchronous  earth  coverage  antennas. 

A  single  corrugated  horn  hybrid-mode  horn  antenna  optimized 
for  earth  coverage  built  in  Lincoln  laboratory  has  a  minimum  edge 
gain  (0  ■  9°)  of  17.8  dBl.  This  horn  had  an  almost  ideal 
Gauss lan-shaped  beam,  and  represents  nearly  an  optimum  design  for 
a  conventional  beam  shape.  We  will  use  this  value  of  17.8  dBl  as 
our  reference  for  measuring  any  improvement  in  antenna 
performance  achieved  in  going  to  the  shaped-beam  approach. 

III.  DOUBLE  FLARE  HORN  ANTENNA 

Concave  antenna  beam  patterns  can  be  obtained  by  conical 

3 

horns  loaded  with  dielectric  sheets  or  dual-mode  corrugated 
horns  .  Here  we  will  discuss  the  double-flare  horn.  Timofeyeva 
first  introduced  this  type  of  horn  and  presented  its  near  field 
patterns  for  Cassegrain  feed  applications.  The  shadow  region 


formed  due  to  the  flare-angle  change  provides  the  aperture 


amplitude  taper  and  the  phase  difference  across  the  aperture 
indents  the  peak  of  the  radiation  pattern  which  otherwise  would  be 
a  pencil  beam.  The  horn  characteristics  can  be  analyzed  from  a 
dual-mode  point  of  view^.  The  TMii“/TEii®  phase  conversion  factor 
is  given  by  Tomiyasu^.  The  differential  phase  relation  between 
the  two  modes  is  found  in  Potter's  original  work  on  dual-mode 

O 

horns  .  This  mode  conversion  by  changes  in  the  flare  angle  allows 
one  to  use  a  much  shorter  differential  phase  correction  section 
(90°  versus  300°  or  660°  for  the  improved  and  the  original  Potter 
horn,  respectively).  Therefore,  the  bandwidth  is  broadened  to 
well  over  10%  . 

A  Q-band  double-flare  horn  was  designed  and  built.  Figure  2 
gives  the  pertinent  dimensions.  Radiation  patterns  and  gains  were 
iraasured  from  A2  to  A7  GHz  in  0.5  GHz  steps.  In  this  frequency 
range,  all  measured  radiation  patterns  have  a  1  to  2  dB  dip  in  the 
bores ight  direction  of  the  horn.  The  E-  and  H-plane  cuts  follow 
each  other  fairly  close  to  within  1  dB.  Figure  3  is  a  typical 
measured  pattern  at  the  frequency  of  44.5  GHz,  Preliminary 
measured  gain  on  the  boreslght  axis  and  the  E-  and  H-plane  average 
gain  at  a  9°  angle  from  axis  are  given  In  Table  1. 


A 

.'•'A' 
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The  measured  E-  and  H-  plane  differences  In  gain  at  9®  off 


axis  over  the  frequency  range  of  43.5  to  45.5  GHz  are  less  than 
0.5  dB.  The  average  gain  is  thus  a  good  approximation  to  the 
actual  gain.  Compared  to  the  9°  off-axis  gain  of  the  single-mode 
optimized  corrugated  horn  of  17.8  dBl,  a  minimum  Improvement  of 
1.1  dB  has  been  achieved  with  the  double-flare  horn.  Some  studies 
of  the  effects  of  parameter  changes  on  the  electrical 
characteristics  of  the  horn  are  being  continued. 

IV.  HORN  LENS  ANTENNA 

Another  shaped  beam  earth-coverage  antenna  is  an 

Q 

angle-limited  corrugated  horn  with  a  beam-shaping  zoned  lens  . 

A  Ji(x)/x  aperture  distribution  with  3  phase  reversal  zones  Is 
used  in  the  design  (where  x  »  6n,r,  ra  -  3,  Bn  -  10.17,  r  » 


normalized  aperture  radius).  Geometric  optics  and  ray  tracing 
techniques  are  used  In  accordance  with  scalar  fields  from  the 
feed.  Figure  4  Is  a  scaled  drawing  of  the  horn-lens.  The 
radiation  patterns  of  this  antenna  have  been  measured  from  39  to 


46  GHz.  At  the  high  frequency  end,  the  radiation  pattern  turned 
out  to  be  too  narrow.  The  best  frequency  range  over  a  2  GHz 
bandwidth  was  found  to  be  from  39.5  to  41.5  GHz.  Figure  5  shows 
the  measured  E-  and  H-  plane  patterns  at  40.5  GHz.  At  this 
freuency,  the  lens  diameter  to  wavelength  ratio  is  15.3.  The 
measured  on-axls  and  the  E-  and  H-  plane  average  gains  at  a  9° 
angle  from  boresight  axis  are  given  in  Table  II. 

TABLE  II.  Measured  Gain  of  Horn-lens  Antenna  at  9“ 


Frequency 

On-axls  Gain 

9“  Off-axis  Gain 

GHz 

dBi 

dBi 

39.5 

21.15 

19.65 

40.0 

21.00 

19.60 

40.5 

21.10 

19.60 

41.0 

20.90 

20. 1 

41.5 

21.75 

20.25 

Compared  to  the  9®  off-axis  gain  of  the  single-mode  optimized 
corrugated  horn  of  17.8  dBi,  a  minimum  improvement  of  1.8  dB  has 
been  achieved  with  the  horn-lens  antenna. 

V.  CONCLUSIONS 

Two  practical  broadband  earth-coverage  horn  antennas  were 
designed,  built  and  tested.  Definite  improvements  of  1.1  and  1.8 
dB  over  the  performance  of  the  optimized  simple-mode  horn  are 


achieved.  The  double-flare  horn  antenna  has  a  very  simple 
geometry  which  can  be  useful  at  submll lime ter  wavelengths.  The 
horn-lens  antenna  has  very  good  pattern  symmetry  and  broad 
bandwidth  characteristics.  Its  dielectric  weight  and  loss  can  be 
reduced  by  special  zoning  techniques,  which  would  however  reduce 
the  bandwidth. 
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IN  A  PARALLEL  PLATE  WAVEGUIDE  -  EXPERIMENT 
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ABSTRACT 

Circular  ring  arrays  of  coaxially-fed  monopole  elements  between 
two  circular  parallel  plates  with  E-plane  (radial)  flare  are  appealing 
for  applications  requiring  scanning  of  a  uniform  beam  360  degrees  in 
azimuth. 

The  performance  of  a  circular  array  of  coaxially-fed  monopole 
elements  radiating  into  an  infinite  parallel  plate  region  has  been 
analyzed  in  [I]  and  (2). 

This  paper  describes  the  experimental  effort  and  presents  measured 
data  for  active  impedance,  coupling  coefficients  and  element  patterns. 

A  6A-eleraent  circular  ring  array  of  coaxially-fed  monopoles  in  a 
parallel  plate  wa/eguide  was  constructed.  The  measured  values  show 
excellent  agreement  with  theoretical  results.  This  agreement  strongly 
suports  the  validity  of  the  analysis,  and  furnishes  a  firm  basis  for  a 
systematic  and  accurate  array  design. 

This  work  was  supported  by  the  Rome  Air  Development  Center  under 
project  No.  2305J303. 


1.  INTRODUCTION 


Circular  ring  arrays  of  coaxially-fed  monopole  elements  between  two 
circular  parallel  plates  with  E-plane  (radial)  flare  are  of  interest  for 
various  antenna  aplications.  In  particular  they  are  very  attractive 
for  radar  and  communication  applications  because  of  their  uniform 
circumferential  radiation  characteristic  and  360  degree  azimuth  coverage 
capabilities.  In  view  of  its  simplicity,  low  cost,  polarization 
purity,  reasonably  wide  bandwidth  and  power  handling  capability,  the 
coaxially-fed  linear  monopole  is  an  attractive  choice  for  an  array 
element  in  a  parallel  plate  waveguide. 

Performance  of  a  coaxially-fed  monopole  element  in  a  circular  ring 
array  radiating  into  an  infinite  parallel  plate  region  has  been  analyzed 
in  [1,2]. 

To  validate  the  theory,  a  64-element  circular  ring  array  of  coaxially 
-fed  monopoles  in  a  parallel  plate  waveguide,  and  a  one  and  two-half 
element  waveguide  simulator  were  constructed.  The  array  and  element 
construction  is  described.  An  illustrative  method  of  selection  of 
optimal  geometry  for  element  match  in  an  array  environment,  as  well  as 
the  dependence  of  Important  design  parameters  on  the  array  and  element 
geometry,  is  discussed.  The  measured  data  is  presented  for  active 
impedance,  coupling  coefficients,  and  element  patterns.  The  measurements 
show  an  excellent  agreement  with  theoretically  predicted  results. 

These  results  strongly  support  the  validity  of  the  previously  reported 
analysis,  and  furnish  a  firm  basis  for  accurate  design  of  a  circular 
ring  array  radiating  into  a  parallel  plate  waveguide  region. 


2.  ARRAY  DESCRIPTION 


The  experimental  array  is  shown  in  Figs.  1  and  2.  The  array  contains 
64  elements  between  two  circular  parallel  plates  with  an  E-plane  (radial) 
flare  as  schematically  shown  in  Fig.  3.  The  flare  is  to  simulate  an 
infinite  parallel  plate  region.  In  Fig.  3  the  top  plate  is  partially 
removed  to  display  the  section  of  the  array.  Fig.  4  shows  schematically 
the  top  and  side  view  of  the  array  with  relevant  dimensions. 

The  circular  ring  array  of  coaxially-fed  monopoles  of  length  I  *s 
in  a  parallel  plate  waveguide  of  height  h  (only  the  TEM  mode  propagates 
since  h  <  X/2).  The  ring,  of  radius  B,  contains  N  equispaced  identical 
monopoles,  each  located  coaxially  a  distance  s  above  a  perfectly 
conducting  circular  cylindrical  surface  of  radius  A.  The  probe  radius 
is  a  <<  \  while  the  inner  and  outer  radii  of  the  coaxial  lines  are  a 
and  b,  respectivelly.  The  coaxial  feed-lines  are  filled  with  teflon 
(er=2)  and  have  characteristic  impedance  Zc=50  ohms. 

The  array  element  spacing  is  chosen  to  be  d^O.A  (0.945")  and 
the  parallel  plate  separation  h=0.369  Xj.  (0.872")  (which  corresponds 
to  the  height  of  the  standard  C-band  rectangular  waveguide)  at  a  center 
frequency  f(,=5.0  GHz.  The  probe  length  and  probe  to  cylindrical 
ground  distance  were  selected  so  that  the  equivalent  infinite  linear 
array  (also  called  the  reference  array)  shown  in  Fig.  5  would  be  matched 
at  broadside  scan  l.e.:  t  =0.233  Xj.  (0.550"),  s  =  0.163  Xj.  (0.385"). 

The  array  elements  shown  in  Fig.  6  were  Omni  Spectra's  Flange  Mount 
Jack  Receptacles  with  a  =  0.01  X^,  (0.025")  and  b  =  0.034  X^.  (0.081"). 
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To  observe  Che  creeping  wave  grating  lobe  ripple  effect  in  the 
element  pattern,  an  Inter-element  spacing  d  =  0.6  Xj.  (0.945")  was 
also  chosen.  This  was  achieved  simply  by  Increasing  center  the  frequency 
from  1^=5. 0  GHz  to  f(,=7.5  GHz,  while  retaining  the  same  top  and  bottom 
parallel  plates  with  radial  flare.  Probes  and  the  cylindrical  ground 
(parallel  plate  spacer)  however,  were  replaced  in  order  to  match  the 
corresponding  reference  array  at  broadside  scan.  The  new  dimensions 
were:  I  =0.25  X^,  (0.393"),  h  =  0.369  X^  (0.580"). 

In  summary,  coupling  coefficients  and  element  patterns  of  a 
64-element  circular  ring  array  of  coaxial  monopole  elements  radiating 
into  a  parallel  plate  region  were  measured  for  two  array  element 
spacings  (a)  d  =  0.4  X^.  and  (b)  d  =  0.6  X^..  In  both  cases,  equivalent 
infinite  linear  arrays  were  matched  at  broadside  scan.  The  respective 
array  element  dimensions  are: 


r  V 
r.v?. 


(a)  fc  =  5.0  GHz  ( X^  =  2.360") 
d  =  0.4  Xc  (0.945") 

N  =  64  (kB  =  25.6) 

A  =  3.915  Xc  <9.240”) 

B  =  4.078  Xc  (9.625") 
h  »  0.369  Xc  (0.872") 

I  =  0.233  Xc  (0.550") 
s  =  0.163  Xc  (0.385") 
a  -  0.010  Xc  (0.025") 
b  =  0.034  Xc  (0.081") 

Zc  -  50  n  (Cr  -  2) 


(b)  fc  =  7.5  GHz  (Xc  =  1.574”) 
d  =  0.6  Xc  (0.945") 

N  =  64  (kB  =  38.4) 

A  =  5.870  Xc  (9.240") 

B  =  6.115  Xc  (9.625") 
h  -  0.369  Xc  (0.580") 

I  »  0.250  Xc  (0.393") 
s  =  0.245  Xc  (0.385") 
a  =  0.016  Xc  (0.025") 
b  -  0.051  Xc  (0.081") 

Zc  ■=  50  n  (Er  -  2) 
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3.  MONOPOLE  ELEMENT  DESIGN 


The  array  elements  sliown  In  Fig.  6  were  Omni  Spectra's  Flange 
Mount  Jack  Receptacles  (model  No.  ?04  CC,  part  No.  2052-1201-00)  with 
a  =  0.025"  and  b  =  0.081".  The  element  dimensions  ( i  and  s)  were 
determined  by  the  requirement  that  an  equivalent  infinite  linear  array 
would  be  matched  at  broadside  scan  (<i>Q=0°)  and  center  frequency  f^. 

For  this  purpose,  Figs.  7  and  8  display  a  contour  plot  of  the  magnitude 
of  the  active  reflection  coefficient  vs.  probe  length  (l/\)  and  vs. 
probe  to  ground  distance  s/X,  for  both  element  spacings  (a)  d/X^  =  0.4, 
and  (b)  d/X^  =  0.6,  respectively.  It  is  seen  that  in  case  (a) 
for  1/ X  =  0.233,  s/X  =  0.163,  and  in  case  (b)  for  i/\  =  0.250, 
s/X  =0.245,  Che  active  reflection  coefficient  =  0. 

Figs.  9  and  10  compare  active  impedance  dependence  on  scan  angles 
(jig  at  fg  of  a  circular  array  with  those  of  an  equivalent  linear  array. 

In  the  circular  array  case,  discrete  values  of  (t>o  are  given  by 
sin  *0  =  u/kB,  where  u  =  l,2,...,N-l.  The  Smith  chart  normalization 
is  50  ohms.  A  region  of  overlap  may  be  noticed  (0°  <  (>q  <  30°)  which 
justifies  the  above  element  matching  scheme.  For  >  30°,  curves 
gradually  depart  from  linear  ones,  due  to  the  curvature  effect. 

For  the  above  geometry.  Fig.  11  shows  the  active  impedance  of  the 
reference  linear  array  at  broadside  scan  vs.  frequency.  It  is  seen  that 
the  frequency  bandwidth  corresponding  to  a  VSWR  of  2:1  is  approximately 


4.  WAVEGUIDE  SIMULATOR 


To  validate  the  numerical  results  for  active  admittance  a  one  and 
two-half  element  waveguide  simulator  was  constructed  [3,4,51.  The  single 
mode  simulator  for  array  (a)  is  shown  in  Figs.  12  and  13.  The  simulator 
waveguide  dimensions  as  shown  in  Fig.  14  are  2d  x  h,  i.e.,  approximately 
1.872"  X  0.872"  which  corresponds  to  standard  C-band  rectangular  waveguide. 
The  waveguide  was  terminated  in  a  matched  load  with  a  VSWR  <  1.02  over 
the  frequency  band  4  to  6  GHz.  In  view  of  Fig.  15  and  the  relation 


sin  <t)o 


it  is  seen  that  the  device  simulates  scan  conditions  from  52°  off-broadside 
at  4.0  GHz  through  40°  off-broadside  at  *=  5.0  GHz  to  32°  off-broadside 
at  6.0  GHz. 

The  active  reflection  coefficient  measurement  was  performed  with  an 
HP-8410  Network  Analyzer.  The  measurement  test  setup  is  shown  in  Fig.  16. 

Fig.  17  shows  the  theoretical  and  the  measured  active  impedance 
vs.  frequency.  The  Smith  chart  normalization  is  50  ohms.  Excellent 
agreement  between  theoretical  and  measured  results  may  be  observed 
across  the  operating  band.  The  two  results  for  the  active  reflection 
coefficient  differ  less  than  IX  in  magnitude  and  less  than  3  degrees  in 
phase. 
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5.  COUPLING  COEFFICIENTS 


Coupli7ig  coefficients  for  the  circular  array  (b)  and  its  Infinite 
linear  reference  array  are  presented  in  Fig.  18.  As  expected  the  linear 
array,  coupling  coefficients  decay  monotonically  with  element  number  p. 
Also,  as  will  be  seen  in  the  next  figure,  for  elements  distant  from 
the  excited  monopole  (p=0)  the  coupling  is  primarily  due  to  the  parallel 
plate  guide  TEM  mode,  the  amplitude  of  which  decays  as  l//r.  As 
expected  in  the  circular  array,  the  coupling  coefficients  initially 
follow  that  of  the  linear  array.  For  elements  far  from  the  excited  one 
(p=0),  the  coupling  coefficients  decrease  exponentially,  which  indicates 
that  the  main  contribution  is  primarily  due  to  a  single  creeping  wave. 

Fig.  19  exhibits  the  phase  of  the  coupling  coefficients  for  the 
array  geometry  of  Fig.  18.  For  convenience  we  define  coupling 
coefficients  SP  as 

“j  itip 

SP  =  jsPj  e  (2a) 


where 


(Jip  =  kdp  +  Aiip 


(2b) 


The  second  term  in  (2b)  is  plotted  in  Fig.  19  while  the  first  term 
represents  the  phase  delay  (for  linear  array)  from  the  reference 
element  (p“0)  to  element  p  of  the  TEM  mode  in  the  parallel  plate 
waveguide.  Thus,  for  elements  close  to  the  excited  one  the  coupling 
is  due  to  the  TEM  mode,  plus  contributions  of  higher  non-propagating 


parallel  plate  waveguide  modes.  For  elements  further  removed,  the 
coupling  is  primarily  due  to  the  TEM  mode. 

Coupling  coefficients  were  measured  down  to  -55  dB  in  amplitude 
which  corresponds  approximately  to  the  seventh  element  from  the  excited 
reference  element  (p=0).  Beyond  this  low  amplitude  level,  reflection 
from  the  aperture  of  the  E-plane  flare  influenced  the  measured  results. 

Figs.  20  and  21  present  a  comparison  between  the  measured  (open 
circles)  and  theoretical  (solid  circies)  coupling  coefficient  values 
in  amplitude  and  phase,  respectively.  The  difference  between  the  two 
results  is  less  than  1.0  dB  in  amplitude  and  less  than  5  degrees  in 
phase. 

6.  ELEMENT  AMPLITUDE  PATTERN 

Fig.  22  shows  the  amplitude  of  element  (field)  pattern  of  circular 
array  (b)  and  its  equivalent  infinite  linear  array.  It  is  observed 
that  both  patterns  exhibit  a  substantial  drop-off  near 

X 

*EGL  “  si*'  ( - 1)  •  O) 

d 

In  the  linear  array  case  this  drop-off  is  caused  by  an  end-fire  grating 
lobe  condition  and  in  the  cylindrical  array  by  its  quasl-linear 
counterpart  as  discussed  in  [6]  for  a  cylindrical  array  of  axial 
strip-dipoles.  The  curves  of  Fig.  22  also  exhibit  a  ripple  in  the 
broadside  region,  whose  amplitude  diminishes  with  closer  element 
element  spacings  and  becomes  negligible  for  d/X  0.5.  The  ripple 
is  due  to  the  Interference  of  the  direct  single  element  radiation 


(with  linear  array  element  pattern)  with  grating  lobes  of  corresponding 
quasl-linear  subarrays  excited  by  the  guided  creeping  wave  as  illustrated 
in  Fig.  23. 

The  element  pattern  was  measured  on  the  far-field  range.  The  axis 
of  rotation  was  at  the  location  of  the  receiving  (reference)  element 
which  was  connected  to  a  20/20  Scientific  Atlanta  Antenna  Analyzer 
while  other  elements  were  terminated  in  50  ohm  coaxial  loads. 

Figs  24  and  25  show  a  comparison  between  measured  (solid  curve) 
and  the  calculated  (dashed  curve)  element  pattern  amplitude  for  both, 

(a)  and  (b)  arrays,  respectively.  ’='rom  these  figures  it  is  seen  that 
in  the  shadow  region  the  patterns  fall  off  exponentially  which  indicates 
that  it  is  primarily  due  to  a  single  creeping  wave  with  an  angular 
attenuation  constant  proportional  to  (kA)^/^. 

In  the  theoretical  result  the  ripple  In  the  =  180°  region, 
similarly  to  that  found  on  a  conducting  cylinder,  is  a  result  of  the 
Interference  of  two  creeping  waves  traveling  in  opposite  directions 
around  tie  cylinder.  Unfortunally  it  could  not  be  observed  experi- 
mentaly  because  of  reflections  from  the  aperture  of  the  E-plane  flare 
that  limited  reflection-free  amplitude  pattern  measurements  to  -30  dB 
level . 

For  better  exposition  of  creeping  wave  grating  lobe  ripples  in 
the  broadside  region  of  the  element  pattern  in  Fig.  25,  Fig.  26  shows 
the  same  curves  on  a  linear  scale.  Excellent  agreement  between  measured 
and  theoretical  results  may  be  observed. 


7.  ELKMENT  PHASE  PATTERN 


Design  of  a  cylindrical  array  requires  knowledge  of  the  element 
phase  center  location.  It  was  found  that  the  flattest  element  pattern 
phase  in  the  broadside  region  is  obtained  for  the  phase  reference 
location  near  the  monopole. 

Fig.  11  exhibits  the  calculated  element  pattern  phase  for  the 
geometry  of  Fig.  21,  One  observes  that  for  both  the  circular  and  its 
equivalent  linear  array,  the  phase  varies  only  by  a  few  degrees  up  to 
the  EGL  drop-off.  Thus  the  EGL  position  essentially  determines  the 
limit  of  usefulness  of  the  element,  both  in  amplitude  and  phase. 

Figs.  28  and  29  show  a  comparison  between  measured  (solid  curve) 
and  the  calculated  (dashed  curve)  element  pattern  phase  of  (a)  and 
(b)  circular  arrays,  respectively.  Excellent  agreement  between 
measured  and  the  theoretical  results  may  be  observed. 


8.  CONCLUSIONS 

The  close  agreement  of  the  experimental  and  the  theoretical 
results  for  active  impedance,  coupling  coefficients,  and  element 
patterns  strongly  supports  the  validity  of  the  analysis  previously 
reported,  and  furnishes  a  firm  basis  for  the  matched  element  design 
method  that  was  developed. 

The  knowledge  of  the  element  patterns  in  the  array  environment 
and  of  the  phase  center  location  which  has  been  experimentaly  verified 
to  be  near  the  element  serves  to  Improve  the  design  accuracy  of 
conformal  arrays. 
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Fig.  7  Linear  reference  array  (a)  active  reflection  coefficient 
-amplitude  vs.  probe  to  ground  distance  and  vs.  probe 
length  (d/Xj,  -  0.4,  (1*^  «  0°) 


Fig.  9  Active  impedance  of  circular  array  (a)  and  its  equivalent 


linear  array  vs.  scan 
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Fig.  15  Crossectlon  of  waveguide  simulator  pertaining  to 
relation  between  frequency  and  scan  angle 
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Fig.  18  Coupling  coefficients  amplitude  for  circular  array  (b) 


and  the  reference  linear  array 
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Coupling  coefficients  phase  for  circular  array  (b) 
and  the  reference  linear  array 
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Fig.  21  Theoretical  and  experimental  phase  of  coupling 
coefficients  for  circular  array  (b) 


30 


572 


element  gain  /y'27rd/x 


OBSERVATION  ANGLE,  <t>  (DEG) 


Fig.  22  Amplitude  of  element  pattern  of  circular  array  (b) 
and  its  equivalent  Infinite  linear  array 
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Fig.  24 
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Theoretical  and  experimental  element  pattern  amplitude 
for  circular  array  (a)  (log  scale) 
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Fig.  25  Theoretical  and  experimental  element  pattern  amplitude 
for  circular  array  (b)  (log  scale) 
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Fig.  28  Theoretical  and  experimental  element  pattern  phase 
for  circular  array  (a)  {d/\.  ■  0.4) 
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Fig.  29  Theoretical  and  experimental  elemnt  pattern  phase 
for  circular  array  (b)  (d/l^  «  0.6) 
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1 .  ABSTRACT 

The  VHF  antenna  system  described  here  will  be  on  the  Global 
Positioning  Satellite  System.  The  antennas  are  designed  to  re¬ 
ceive  circularly  polarired  signals  with  a  design  goal  of  a  6  dB 
axial  ratio.  Due  to  limited  space  on  the  satellite,  they  are 
electrically  small  in  size.  The  low-band  antenna  (43  percent 
bandwidth  at  VHF)  is  a  set  of  four  meandering  monopoles  (designed 
by  Dr.  Alan  W.  Love,  based  on  the  meandering  monopole  patent  by 
Mr.  Howard  E.  King  and  Dr.  Jimmy  L.  Wong)  fed  by  a  quadripole 
network.  They  are  arranged  on  the  rear  bulkhead  of  the  satellite 
to  provide  for  best  circular  polarization.  The  high-band  antenna 
has  a  bandwidth  of  32  percent  at  VHF  and  is  positioned  on  the 
forward  bulkhead  on  an  extended  ground  plane.  The  high-band 
antenna  is  a  crossed-dipole  turnstile  antenna  (designed  by  Dr. 
Alan  W.  Love)  fed  by  a  quadrature  hybrid.  Both  antennas,  using 
directional  couplers,  have  self-test  verification  provisions.  An 
interesting  feature  is  that  due  to  the  electrically  small  size  of 
the  low  band  antennas,  the  antennas  act  as  probes  which  sample 


currents  on  the  satellite  surface  and  solar  panels.  The  antenna 
system  receives  well  from  directions  toward  and  away  from  the 
earth.  The  VSWR,  polarization,  gain,  and  reception  patterns  must 
be  measured  with  the  antennas  mounted  on  a  full  scale  mockup  of 
the  satellite.  The  antennas  are  developed  to  have  a  large 
bandwidth,  good  axial  ratio,  and  an  acceptable  gain  and  VSWR. 

2.  HIGH  BAND  VHP  ANTENNA  SYSTEM 

2.1  Description  of  Antenna  and  VSWR  Performance 

The  high-band  antenna  has  an  operating  bandwidth  from 
109-151  MHz.  It  Is  a  crossed-dipole  turnstile  antenna  which  Is 
fed  by  a  quadrature  hybrid  as  In  the  circuit  schematic  Figure  1. 
It  Is  situated  on  the  forward  bulkhead  of  the  satellite  with  an 
extended  ground  plane  as  may  be  seen  In  Figure  2.  The  antenna 
dipole  element  Is  44"  In  length  which  Is  (.40)  wavelengths  at  109 
MHz.  The  parasitic  element  which  Is  provided  for  greater  band¬ 
width  Is  40"  In  length  which  Is  .37  wavelengths  at  109  MHz  .  It 
Is  seen  that  this  Is  electrically  short;  however,  this  Is  neces¬ 
sary  to  prevent  Interference  with  the  solar  panels.  The  high- 
band  antenna  Is  depicted  to  Figure  3.  The  Voltage  Standing  Have 
Ratio  was  within  the  specification  of  6:1.  This  was  accomplished 
by  using  an  adjustable  balun.  The  fully  extended  position  was 
chosen  for  best  performance  as  shown  In  Figure  4.  The  VSMR  of 
the  high-band  antenna  for  both  ports  is  within  the  specification 
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frequency  band.  Solar  panel  position  did  not  appear  to  have  much 

effect  on  the  axial  ratio.  The  axial  ratios  for  the  e  =  7*,  and 

e  s  14*  were  also  satisfactory  due  to  wide  beamwidth  and  are  not 
discussed  In  this  paper. 

2.4  High  Band  Antenna  Gain 

The  gain  measurements  are  also  within  the  specification  of 
-6  dB^.^.  The  gain  tends  to  be  highest  at  around  130  MHz  and 
tends  to  tail  off  as  the  outer  limits  of  the  frequency  band  are 
reached  as  shown  In  Figure  7. 

2.5  High  Band  Antenna  Summary 

The  high-band  crossed  dipole  turnstile  antenna  on  the 

extended  ground  plane  was  basically  a  well  behaved  antenna 

system.  The  parasitic  element  helped  create  a  larger  bandwidth, 
along  with  the  "tuning"  of  the  movable  balun.  The  solar  panels 
for  the  high  band  antenna  did  not  appear  to  have  any  major  effect. 

3.  LOW  BAND  ANTENNA  VHP  ANTENNA  SYSTEM 
3. 1  Description  of  Low  Band  Antenna 

The  low  band  antenna  has  an  operating  bandwidth  from  48-54 
MHz  and  73-75  megahertz.  The  low-band  antenna  Is  a  set  of  four 
meandering  monopoles  fed  by  a  quadripole  network  with  directional 
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Figure  7.  High-Band  Antenna  Gain  (dBic) 


couplers  to  provide  self  test  capability.  The  circuit  schematic 
is  shown  in  Figure  8.  The  antennas  are  configured  on  the  rear 
bulkhead  of  the  spacecraft  to  provide  for  good  axial  ratio.  How¬ 
ever,  the  antennas  are  set  apart  by  104*  and  76*  to  eliminate 
interference  with  solar  panels.  The  position  of  the  antennas  may 
be  noted  in  Figure  9.  The  meandering  monopole  is  utilized  to 
create  a  "longer"  electrical  length  (see  Figure  10)  for  the 
antenna  without  increasing  the  length  of  the  antenna.  The  actual 
length  of  the  vertical  elements  of  the  antenna  is  approximately 
24"  (.lx)  at  48  MHz.  The  total  length  of  the  meandering  monopole 
is  about  110".  The  extra  length  increases  the  efficiency  of  this 
extremely  short  antenna. 

3.2  Low  Band  Antenna  as  Current  Probe  on  Satellite 

The  meandering  monopole  must  be  physically  small  to  prevent 
interference  from  the  solar  panels.  The  concept  is  to  use  the 
monopole  as  a  current  probe  exciting  surface  currents  on  the 
satellite  and  solar  panels.  This  method  is  effective  because  the 
circumference  of  the  satellite  and  solar  panels  are  on  the  order 
of  a  wavelength  which  may  lead  to  resonance.  The  position  of  the 
solar  panels  is  also  extremely  critical  to  the  axial  ratio  and 
gain  measurements. 


3.3  Low  Band  Antenna  VSWR 

The  voltage  standing  wave  ratio  was  expected  to  be  6:1. 

However,  the  antenna  was  changed  from  Its  original  position  to 

allow  for  clearance  of  the  solar  panels.  The  height  of  the 
antenna  over  the  satellite  was  changed  from  five  Inches  to 

2-3/4".  The  antennas  were  also  shifted  from  a  separation  of 

ninety  degrees  apart  to  that  of  a  104*  and  76*  apart,  and  the 
monopole  was  shifted  In  over  the  Interconnect  assembly  Instead  of 
being  two  Inches  away.  These  conditions  resulted  In  a  voltage 
standing  wave  ratio  which  did  not  meet  expectations.  At  48  htiz, 
the  VSWR  was  12:1  and  from  73-75  MHz  the  VSWR  was  usually  8:1. 
The  VSWR  on  the  satellite  mockup  Is  provided  In  Figure  11. 

3.4  Low  Band  Antenna  Axial  Ratio 

The  meandering  monopoles  were  positioned  on  the  satellite  to 
provide  for  reception  of  circular  polarized  signals  with  a  design 
goal  of  6  d8  axial  ratio.  As  seen  In  Figure  12,  the  axial  ratio 
was  measured  at  the  various  solar  panel  positions  and  at  the 
on-axis  position.  Axial  ratio  was  also  measured  at  e  x  7*  and  e 
X  14*  and  was  slightly  degraded.  However,  they  are  not  discussed 
In  this  paper.  The  low-band  antennas  act  as  current  probes 
exciting  the  satellite  and  solar  panels.  With  right  hand  cir¬ 
cular  polarization,  the  worst  case  solar  panel  positions  appears 
to  be  when  the  solar  panels  are  perpendicular  and  parallel  to  the 
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satellite.  Best  case  conditions  appear  to  be  with  solar  panel 
angles  at  4b*  and  135*.  Perhaps  this  could  be  due  to  interaction 
of  the  fields  at  the  positions  and  some  of  the  vectors  cancelling 
out.  Left  hand  circular  polarization  seemed  to  follow  the 
general  trend  as  established  by  RHCP,  however,  the  axial  ratios 
are  slighter  higher.  Also  for  73-75  MHz,  the  axial  ratio  for 
solar  panel  position  at  0*  is  almost  linear. 

3. 5  Low  Band  Antenna  Gain 

The  gain  is  shown  in  Figure  13.  For  most  cases,  the  gain  is 
within  the  specification  of  -6  dBic.  48  MHz  has  lower  gain  due 
to  its  high  VSWR  and  73-75  ftiz  due  to  the  high  axial  ratio. 
Solar  panel  angle  seems  to  have  only  a  slight  effect  on  the  gain. 

3.6  Low  Band  Antenna  Self  Test  Verification 

The  self  test  verification  used  the  directional  couplers  in 
reversed  position.  The  result  should  be  as  indicated  in  Figure 
14  for  one  antenna  with  varying  VSWR  versus  three  identical  impe¬ 
dance  antennas.  With  all  antennas  matched,  the  self  test  verifi¬ 
cation  will  result  in  a  transmission  loss  of  approximately  -55  dB 
(see  Figure  15).  With  one  antenna  open,  the  transmission  loss 
is  -35  dB  (see  Figure  16).  Thus  a  failure  of  one  antenna  may  be 
detected.  However,  if  two  antennas  (or  all  four  antennas)  feed¬ 
ing  into  the  same  quadrature  hybrld(s)  fail  and  have  the  same 
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resulting  Impedance,  a  null  condition  will  also  result.  The  self 
test  verification  will  verify  If  any  antennas  fall  In  operation 
with  the  exception  mentioned  above. 


3.7  Summary  of  VHP  Satellite  Antenna  System 

The  high-band  and  low-band  antennas  perform  well  as  antenna 
systems  on  the  satellite.  The  low-band  antenna  does  not  meet  the 
VSWR  specification  due  to  Its  shift  In  position  In  order  to  pro¬ 
vide  clearance  for  the  solar  panels.  However,  the  high  and  low- 
band  antennas  for  most  cases  do  meet  the  system  requirements. 
Several  design  considerations  for  large  bandwidth  antennas  at  VHP 
are  the  shape  of  the  satellite,  the  position  of  the  solar  panels, 
to  provide  for  clearance,  and  the  position  of  the  antenna  above 
the  satellite  bulkhead. 
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A  SIMPLE  TECHNIQUE  TO  CORRECT  FOR  CURVATURE  EFFECTS  ON 
CONFORMAL  PHASED  ARRAYS 
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Peter  Franchl 
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Rome  Air  Development  Center 
Hanscom  AFB,  MA  01731 

ABSTRACT 

The  curvature  of  conformal  phased  array  antennas  distorts  the  radia¬ 
tion  pattern  when  the  beam  Is  scanned.  The  distortion  occurs  because  the 
amplitude  taper  appears  skewed  when  the  beam  is  scanned.  If  an  anti¬ 
symmetric  taper  is  added  In  the  proper  proportion  to  a  symmetric  taper, 
the  symmetry  of  the  original  taper  can  be  largely  recovered.  Typically, 
the  sum  and  difference  aperture  illuminations  of  an  antenna  can  provide 
the  needed  symmetric  and  anti-symmetric  functions.  A  relatively  simple 
and  practical  technique  of  using  the  sum  and  difference  networks  in  com¬ 
bination  to  correct  for  curvature  effects  is  demonstrated.  This  technique 
uses  a  variable  power  divider  to  distribute  Che  power  in  a  prescribed 
proportion  at  an  arbitrary  scan  angle  to  the  sum  and  difference  channels. 
The  re-distributed  power  at  the  elements  restores  to  a  large  degree  the 
original  aperture  illumination  for  the  scanned  direction.  The  optimum 
amount  of  signal  to  be  distributed  to  achieve  the  maximum  restoration  is 
found  as  a  function  of  scan  angle  and  curvature.  The  technique  can  be 
used  to  improve  any  sum  and  difference  illuminations.  In  this  paper  cos¬ 
ine  and  sine  illumination  functions  respectively,  are  used  for  the  sum 
and  difference  illuminations  to  demonstrate  the  technique. 
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1.  INTRODUCTION 

The  development  of  phased  array  antennas  has,  for  the  most  part, 
been  for  antennas  on  planar  surfaces.  The  large  phased  array  antennas, 
for  example.  Pave  Paws,  Aegis,  Cobra  Dane,  and  Cobra  Judy,  are  planar 
arrays.  Interest  has  developed  In  phased  arrays  on  curved  surfaces  be¬ 
cause  of  the  anticipated  need  for  antennas  conformal  to  aircraft.  When 
scanning,  curvature  has  a  detrimental  effect  on  the  antenna's  radiation 
pattern;  the  pattern  becomes  distorted  and  the  sidelobes  are  increased. 

The  effect  is  Influenced  by  factors  such  as  the  element  pattern,  the  ex¬ 
tent  of  curvature,  and  the  antenna  size,  and  becomes  especially  pronounced 
as  the  scan  angle  Increases. ^ 

The  purpose  of  this  paper  is  to  demonstrate  a  relatively  simple  and 
practical  technique  to  correct  for  curvature  effects.  A  straightforward 
approach  to  correct  the  antenna  pattern  would  be  to  adjust  the  amplitude 
and  phase  of  each  individual  radiating  element.  This  procedure  would  be 
unwieldy  and  highly  impractical.  Introducing  either  attenuators  or  power 
dividers  into  the  system  in  order  to  control  the  amplitude  of  each  element 
would  greatly  Increase  both  rf  network  and  controlling  network  complexity 
In  addition,  cost,  errors,  failure  rate,  weight,  and  space  limitation 
problems  would  considerably  Increase.  Attenuators  would  also  introduce 
power  loss. 


network.  Consideration  of  the  aperture  taper  in  the  projected  aperture 
plane  gives  insight  into  the  cause  of  the  pattern  distortion.  At  scan 
angles  other  than  zero  degrees,  the  projected  aperture  taper  appears 
skewed;  the  greater  the  scan  angle,  the  greater  the  skewness.  By  adding 
in  the  proper  proportion  an  anti-symmetric  taper  to  the  symmetric  taper, 
the  symmetry  of  the  original  broadside  taper  can  be  largely  restored. 

This  correction  can  be  easily  Implemented  with  the  sum  and  difference 
networks  of  the  conformal  phased  array  antenna  and  a  simple  variable 
power  divider. 

2.  EXAMINATION  OF  FIRST  ORDER  CURVATURE  EFFECTS 

The  approach  In  this  study  is  to  investigate  the  effect  of  curvature 
on  an  antenna  radiation  pattern  based  on  geometrical  considerations  only. 
Phenomena,  such  as  creeping  waves  and  mutual  coupling  will  be  neglected. 
The  physical  model  is  an  array  of  a  fixed  number  of  elements  (65  in  this 
case)  and  a  fixed  element  spacing  (X/2).  The  projection  of  the  curved 
array  on  to  a  straight  line  perpendicular  to  the  beam  direction  would  re¬ 
sult  in  a  projected  array  which  is  shorter  than  the  original  planar  array 
The  radiating  elements  on  the  curved  surface  are  positioned  differently 
in  three  dimensional  space  compared  to  the  elements  on  the  planar  surface 
this  difference  distorts  the  wavefront  from  the  equivalent  plane  surface 
thereby  distorting  also  the  radiation  pattern.  Primarily,  this  first 
order  effect  on  the  close-in  sidelobes  will  be  addressed  in  this  paper. 
Grating  lobe  effects  will  not  be  examined. 

Of  the  different  types  of  aperture  Illuminations,  only  one  will  be 
considered  in  order  to  demonstrate  the  effect  of  curvature  and  the 


correction  technique.  The  aperture  will  have  a  cosine  1 1 luniinat ion  with 
each  element  having  a  cos^*"^  (y)  pattern,  where  y  i^  a  general  angular 
variable.  This  element  pattern  approximates  measured  array  element 
patterns  and  is  accurate  for  patch  radiators,  narrow  slots,  or  dipoles  on 
a  cylinder  in  which  tne  E-field  is  along  the  cylindrical  axis.  Blockage 
effects  are  included.  The  element,  radiates  only  up  to  the  line  drawn  from 
it  through  the  next  succeeding  element  (very  nearly  the  tangent  line). 
Although  a  section  of  a  ring  array  is  being  investigated,  the  results 
apply  to  the  principle  plane  cut  of  a  section  of  a  cylindrical  array  form¬ 
ed  by  the  three  dimensional  extension  of  the  ring  array.  There  are  no 
cross  polarization  considerations  in  this  model. 

Figure  I  demonstrates  a  cosine  aperture  illumination  (showing  nine 
elements)  which  changes  at  a  scan  angle,  Bq,  In  the  broadside  direction, 
the  taper  is  a  slightly  distorted,  but  symmetrical  projection  of  a  cos 
(it  n/K)  aperture  amplitude  distribution  for  the  sum  pattern,  where  n  is 
the  element  number  in  which  the  center  element  is  zero,  and  K  is  the  total 
number  of  element  spacings.  At  the  scan  angle,  6q,  the  figure  shows  how 
the  taper  becomes  skewed  due  to  the  displaced  positions  of  the  elements 
along  the  line  of  projection.  For  example,  the  center  element's  position 
is  no  longer  at  the  center  of  the  projection  line.  The  center  element 
retains  the  peak  amplitude  value  of  the  cosine  taper,  and  therefore,  the 
maximum  amplitude  is  shifted  off  center.  The  other  elements  follow 
accordingly.  In  this  manner  the  original  taper  becomes  skewed.  Figure  2 
depicts  the  same  phenomenon  taking  place  for  a  sln(2iT  n/K)  amplitude 


Figure  1.  Cosine  aperture  illumination  on  a  conformal  phased  array  antenna 
At  the  scan  angle,  dg  ,  the  illumination  becomes  skewed. 


Figure  2.  Sine  aperture  illumination  on  a  conformal  phased  array  antenna. 
At  the  scan  angle,  ,  the  illumination  becomes  skewed. 


distribution  for  the  difference  pattern.  As  the  scan  angle  increases,  or 
as  S/R,  S  being  the  surface  length,  and  R  being  the  radius  of  curvature, 
increases,  the  center  element  is  shifted  further  off  center  and  the  asym¬ 
metry  Increases.  This  distance  off  center  serves  as  an  indicator  of  dis¬ 
tortion.  The  above  discussion  dealt  with  the  dependence  of  the  field  at 
the  aperture  due  only  to  excitation  coefficients  of  the  elements  and  not 
to  their  individual  patterns.  The  element  pattern  must,  however,  also  be 
included  in  order  to  obtain  the  total  field.  Therefore,  the  general  dis¬ 
cussion  of  the  field  at  the  aperture  would  include  a  multiplying  factor, 
Ag,  which  is  the  amplitude  of  the  element  pattern  and  is  a  function  of 
scan  angle. 

Figure  3  shows,  for  an  arbitrary  scan  angle,  the  position  of  the 
elements  on  the  projection  plane  and  the  distortion  of  the  amplitude 
taper  as  compared  to  a  cosine  illumination  function  in  that  plane.  The 
center  element  is  at  the  origin  of  the  x,  y  coordinate  system,  and  the 
other  elements  are  numbered  accordingly.  The  cosine  value  of  the  nth 
element  is  cos(tt  n/K)  where,  for  convenience,  the  total  number  of 
elements,  K  +  1,  is  chosen  to  be  odd.  The  distance  of  the  nth  element 
along  the  x  axis  is  found  to  be  given  by  the  following  equation: 

Xn  =  2R  sin(ns/2R)  sin(x/2  -  6  +  ns/2R)  (1) 

where  R  is  the  radius  of  curvature,  6  is  the  scan  angle,  and  s  is  the 
surface  length  between  two  consecutive  elements.  A  plot  of  versus  the 
cos(im/K)  illumination  (the  peak  amplitude  being  normalized  to  one) 
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reveals  the  kind  of  distortion  which  Is  depicted  In  Figure  3.  The  same 
type  of  distortion  argument  holds  for  the  difference  function,  which  Is 
given  by  sln(27m/K).  This  Is  shown  In  Figure  4.  The  cosine  curve 
In  the  projection  plane  Is  shown  In  the  X*,  Y'  coordinate  system.  The 


cosine  curve  here  Is  given  by: 

Y'  =  cos(TrX'/ll)  (2) 

Y'  =  cos[ir(Xfi  -  h)/ll]  (3) 

where 

h  =  R[1  -  cos(Ks/2R)]  sin  0  (A) 

1  =  2R  sln(Ks/2R)  cos  9  (5) 


h  being  the  distance  the  center  element  shifts  off  of  the  geometrical 
center,  I  being  the  length  of  the  array  In  the  projection  plane. 

There  Is  one  more  contributor  to  the  distortion  of  the  radiation 
pattern.  It  Is  the  change  of  the  density  of  the  elements  In  the  projec¬ 
tion  plane.  This  causes  a  change  In  the  field  density.  Increasing  the 
element  density  Is  equivalent  to  Increasing  the  field  that  would  other¬ 
wise  be  present  If  the  elements  were  uniformly  spaced.  Consider  the 
field  of  the  nth  element  to  be  E(n),  and  the  effective  field  E'(n).  F(n) 
will  be  referred  to  as  the  density  function  and  Is  given  by  E'(n)/E(n)  or 

F(n)  -  2t/K(Xn+i  -  Xn-l)  (6) 

with  Xn+i  and  Xn-l  defined  by  equation  1. 
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Figure  3.  Distortion  in  the  projection  plane  of  cosine  amplitude  taper 
on  curved  surface. 
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Figure  4.  Distortion  in  the  projection  plane  of  sin  2>' amplitude  taper 
on  curved  surface. 
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The  effective  amplitude  taper  In  the  projection  plane  Is 
F(n)  Ag  cos(im/K) 

3.  RESTORATION  OF  THE  BROADSIDE  APERTURE  TAPER 

The  cos  (y)  curve  and  the  sin  (2y)  curve  are  both  distorted  In  the 
projection  plane.  By  adding  these  two  functions,  the  original  cosine 
taper  can  be  approximated  In  the  projected  plane.  This  addition  Is 
accomplished  by  using  the  sum  port  and  difference  port  simultaneously. 

The  power  Is  distributed,  without  loss,  by  a  variable  power  divider.^  The 
Idea  behind  a  variable  power  divider  Is  to  split  the  signal  entering  In 
one  arm  and  then,  by  phase  adjustments,  control  the  amount  of  power  going 
Into  Its  two  output  ports.  Each  port  Is  part  of  a  network  which  generates 
across  the  array,  a  cos  (y)  and  a  sin  (2y)  distribution  respectively. 
Figure  5  shows  the  variable  power  divider  In  the  antenna  system.  The 
signal  from  the  sum  port,  I,  is  distributed  between  the  sum  and  differ¬ 
ence  ports,  T',  and  A',  respectively.  A  portion  of  the  signal  in  the 
difference  network  Is  added,  at  the  elements,  to  the  signal  from  the  sum 
network.  For  the  receive  mode,  the  principle  of  reciprocity  applies. 

The  optimum  fraction  of  power  to  be  distributed  to  each  port  can  be  found 
that  achieves  a  best  fit  to  the  cos  (y)  taper  In  the  projection  plane  at 
a  particular  scan  angle.  Generally,  If  A  and  B  are  designated  the  ampli¬ 
tude  of  the  energies  In  the  sum  and  difference  ports  respectively,  then 
A  and  B  will  be  chosen  such  that  A  cos  (y)  +  B  sin  (2y)  approximates 
cos  (y'),  where  y'  Is  the  variable  In  the  projection  plane.  The  best 
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Figure  5.  Variable  power  divider  Incorporated  Into  the  antenna  system. 

fit  would  be  to  obtain  the  minimum  error  between  the  shape  of  the  correct¬ 
ed  curve  and  the  undlstorted  cosine  curve*  Since  we  are  matching  curve 
shapes,  the  curves  are  to  be  normalized;  the  corrected  curve  Is  to  be 
normalized  to  one.  The  difference  between  the  two  curves  at  the  nth 
point  Is  expressed  by: 

dj,  -  G(B)F(n)Ae[A  co8(ira/K)  +  B  sln(2im/K)J  -  co8(irX*/t)  (7) 

G(B)  Is  the  normalizing  factor,  and  varies  with  B. 

G(B)  -  [l/[F(n)Ae(A  co8(im/K)  +  B  8ln(2im/K)) 1 (8) 

The  constraint  In  the  addition  of  the  two  networks  Is  that  the  total 
power  remains  constant.  If  the  Input  power  Is  normalised  to  one,  then 
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Ultimately,  the  equation  used  to  match  the  distorted  curve  to  a  cosine 
curve  with  uniformly  spaced  elements  Is  the  following: 


=  G(B)F(n)Ag[/i-B^  co8(irn/K)  +  B  8ln(2im/K)]  -  cos [ ir(X^-h)/ 1]  (10) 


The  total  deviation  between  the  two  curves  Is  the  sum  of  the  over 

all  n.  Since  the  absolute  value  of  Intererst ,  It  Is  mathemat- 

Ically  convenient  to  use  6^^  Instead  of  6^^.  The  value  to  be  calculated 
Is  the  sum  of  the  squares  of  the  deviations,  c,  which  Is  defined  as: 


K/2  ^ 

n— K/2 


(11) 


Taking  6e/fiB  »  0  yields  Bg,  the  value  of  B  for  which  e  Is  a  minimum. 

The  values  of  Bg  as  a  function  of  scan  angle,  0,  with  S/R  as  a  parameter, 
and  for  equal  to  cos’^^O  (a  cosine  to  the  1.5  power  radiation  pattern) 

are  shown  In  Table  1.  There  were  no  values  chosen  In  which  element  block¬ 
age  was  present.  Blockage  occurs  according  to  the  following  formula: 


>  w/2  -  S/2R  +  s/2R 


(12) 


Since,  when  S/R  -2,  (S  ■  32X,  R  ■  16X,  s  ■  X/2)  blockage  occurs  at 
33.6”,  the  case  of  S/R  «  2  and  6  «  40*  was  not  used.  Figure  6  shotrs  how 
the  error,  e,  varies  with  scan  angle  for  particular  values  of  S/R.  The 
graph  demonstrates  the  extent  of  which  the  error  Is  reduced  by  Introducing 
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Table  1.  Correction  factor,  B,  as  a  function  of  ccan  angle,  and  the 
curvature  factor,  S/R,  for  the  cosine  to  the  1.5  newer  element  radiation 
pattern  case. 


the  correction.  When  there  Is  no  correction,  or  Bq  equals  zero.  It  Is 
observed  that  the  curve  of  c  versus  6  follows  approximately  a  tan^B 
curve.  This  dependency  Is  to  a  large  degree  verified  by  the  graph  in 
Figure  6.  The  equation  Is  approximate,  but  reasonably  accurate.  It 
shows  the  variation  of  pattern  distortion  as  the  scan  angle  Is  changed. 
The  solid  curve  shows  the  extent  of  the  uncorrected  amplitude  taper 
error,  which.  In  turn  Indicates  the  degree  of  distortion  of  the  radiation 
pattern.  The  dashed  curve  Is  the  residual  error  after  correction.  Both 
the  uncorrected  and  the  corrected  error  at  broadside  are  finite  because 
the  broadside  curvature  effects  were  not  readjusted  initially. 
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4.  CURVATURE  EFFECTS  ON  RADIATION  PATTERN 

The  power  radiated  at  an  arbitrary  scan  angle,  9q,  can  be  obtained 
by  considering  the  scan  angle  direction  of  each  element,  and  taking  Into 
account  the  path  length  difference  of  each  element.  The  E-field  In  the 
scan  angle  direction,  0q,  Is: 


n=-K/2 


where  A^CBq)  is  the  nth  element's  amplitude  due  to  the  element  radiation 
pattern  in  the  direction  of  0q. 


Ae(0o)  “  cosP(9q  -  ns/R) 


p  being  a  constant  defining  the  element  pattern.  A^  and  are  the  magni 
tude  and  phase  respectively  of  the  nth  element's  excitation  coefficient. 

A  beam  is  steered  to  the  6q  direction  by  having  equal  to  RZ^COq). 

The  value  of  Z^CBq)  is: 


Zn(0o)  =  2R  sin(ns/2R)  co8(ff/2  -  0q  +  ns/2R) 


The  E-field  in  the  0  direction,  relative  to  a  0q  scan  angle  direction  is: 


E(0)  -  )  A^(0)A„eJ'^^n<V  e-J*^*n^®> 

n=-K/2 


The  value  of  the  power,  P,  is  obtained  by  using  equation  (16)  in  the 
following  equation: 


P  -  1/2  |E(0)|2 


(17) 


The  effect  of  curvature  on  the  radiation  pattern  is  predominantly  a  func¬ 
tion  of  Sq*  ns/R,  and  consequently,  Ks/R,  or  S/R,  the  array's  subtended 
angle.  To  demonstrate  the  effect  of  curvature,  the  radiation  pattern,  P, 
is  obtained  as  a  function  of  S/R  (with  6q  constant)  and  as  a  function  of 
Bq  (with  S/R  constant).  Both  the  number  of  elements,  K  -f  1,  and  the 
element  spacing,  s,  are  held  constant. 

K  reasonable  choice  for  the  element  pattern  is  p  >  1.5  since  a  patch 
radiator,  a  slot,  or  a  dipole  in  an  array  exhibits  very  nearly  this 
dependency.  Figures  7a,  7b,  and  7c  show  the  radiation  plots  using  a 
cos^*^(y)  element  pattern  and  cosine  aperture  taper  for  various  values 
of  scan  angle,  Og,  and  subtended  angle,  S/R.  The  number  of  elements 
used  is  65,  the  element  spacing  is  X/2,  and  X  is  .1  meters.  The  correc¬ 
tion  is  most  pronounced  as  the  scan  angle  is  increased  and/or  the  sub¬ 
tended  angle  is  increased.  The  figures  show  the  Increase  of  the  near-in 
sldelobe  levels  as  a  function  of  scan  angle,  and  subtended  angle;  the 
uncorrected  and  corrected  patterns  are  shown.  The  improvement  of  the 
distorted  pattern  is  substantial,  showing  the  correction  works  well.  The 


optimum  sum  pattern  results  In  an  improvement  of  the  difference  pattern 


The  crucial  criterion  of  a  difference  pattern  is  the  mainbeam  null  depth; 
the  deeper  the  mainbeam  null,  the  better  the  pattern.  One  example  of  the 
improvement  of  the  difference  pattern  is  shown  in  Figure  8. 

Table  (2)  depicts  the  Improvement  of  the  mainbeam  null  depth  which 
can  be  seen  to  be  significant.  The  improvement  decreases  slowly  as  the 
curvature  Increases.  There  is  not  a  significant  variation  with  a  change 
of  the  scan  angle.  A  better  difference  pattern  can  be  obtained  by  sacri¬ 
ficing  some  of  the  correction  of  the  sum  pattern. 


0° 

S 

MAINBEAM  NULL 
(-dB) 

IMPROVEMENT 

R 

UNCORRECTED 

PATTERN 

CORRECTED 

PATTERN 

IN  dB 

1.07 

23.20 

30.66 

7.46 

20° 

1.6 

19.70 

25.72 

6.02 

2 

17.78 

22.50 

4.72 

1.07 

19.71 

27.35 

7.64 

o 

0 

1.6 

16.20 

22.40 

6.20 

2 

14.24 

19.27 

5.03 

'  1 

1.07 

17.24 

24.92 

7.68 

0 

o 

1.06 

13.68 

6.38 

Table  2.  Values  of  mainbeam  null  of  uncorrected  and  corrected  difference 
patterns  as  a  function  of  scan  angle,  and  subtended  angle.  A  cosine  to 


the  1.5  power  element  pattern  was  used. 


5.  CONCLUSION 


The  results  demonstrate  the  value  of  the  correction  technique.  It 
Is  a  simple,  practical  technique  expending  no  additional  power  and  yet 
offers  quite  a  high  level  of  Improvement.  It  Is  worthwhile  to  compare  an 
amplitude  correction  curve  with  a  result  obtained  by  Ludwig.^  Ludwig  used 
an  Iterative  algorithm  In  which  the  complex  weights  of  each  element  was 
adjusted  In  order  to  synthesize  an  Ideal  pattern.  Although  he  used  a 
Taylor  aperture  Illumination,  the  two  results  are  still,  to  some  extent, 
comparable.  For  his  case,  using  a  40*  scan  angle  and  S/R  •  0.8,  the  peak 
of  the  corrected  amplitude  taper  shifts.  In  terms  of  his  array,  by  about 
three  and  one  half  elements.  Using  oui  technique,  for  Qq  -  40*  and  S/R  > 
1.07,  there  Is  approximately  the  same  shift.  The  realization  of  the 
Iterative  algorithm  In  an  actual  system  is  highly  improbable,  but  our 
corrector  can  be  easily  implemented.  We  demonstrated  our  technique  using 
a  cosine  aperture  taper.  If  a  Taylor  illumination  was  employed,  we  would 
use  a  Taylor,  Bayllss  addition.  Generally,  any  symmetric  and  antisymmetric 
functions  are  to  be  combined.  Another  study  may  be  to  use  the  technique 
in  order  to  Increase  the  power  In  the  malnbeam  for  large  scan  angles.  In 
that  way,  the  extent  of  large  angle  scanning  can  be  augmented  and  made 
possible  by  the  existence  of  curvature. 
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ABSTRACT 


A  method  for  maximizing  the  slgnal-to-interf erence  ratio  of 
an  antenna  array  while  maximizing  the  directivity  has  been  de¬ 
veloped  in  this  paper.  This  method  is  formulated  for  an  arbi¬ 
trary  array  in  an  arbitrary  Interference  environment.  The 
objective  is  achieved  by  defining  a  null  space  of  a  matrix  A 
associated  with  the  interference  power.  Numerical  examples  for 
planar  arrays  of  short  roonopoles  subjected  to  point  and  distri¬ 
buted  interference  sources  illustrate  the  method. 

I .  INTRODUCTION 

An  antenna  system  when  used  for  receiving  a  desired  signal 
is  facing  in  reality  a  hostile  environment,  which  is  usually 
contaminated  with  noises,  interferences,  and  even  deliberately 
made  jamming  signals.  In  a  communication  center,  such  as  a 
shipboard  system  in  particular,  many  antennas,  supporting  struc¬ 
tures,  and  other  obstacles  are  jampacked  in  a  small  space.  All 
of  them  can  cause  a  serious  deterioration  of  the  intended  signal 
quality.  In  fact,  the  signal  wave  reflected  from  the  turbulent 
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water  or  moving  objects  can  reach  the  antenna  and  have  an  unde¬ 
sirable  fading  effect  on  the  Incoming  signal  directly  from  the 
transmitter.  Even  for  transmitting,  it  may  also  be  desirable  to 
reduce  the  antenna  radiation  below  a  very  low  level  In  some 
directions.  Following  the  work  by  Lo,  Lee,  and  Lee^  on  the 
maximization  of  the  slgnal-to-nolse  ratio  (SNR),  this  paper  con¬ 
siders  the  optimization  problem  for  an  arbitrary  array  In  an 
arlbtrary  Interference  environment,  such  that  the  slgnal-to- 
Interference  ratio  (SIR)  Is  maximized. 

The  Interference  source  can  be  distributed  over  a  finite 
angular  sector,  such  as  the  reflected  signal  from  the  water  or  a 
large  obstacle,  or  concentrated  in  a  single  direction  as  If  It 
were  a  point  source.  The  latter  is  a  special  case  of  the 
former,  when  the  sector  Is  vanishingly  small,  and  the  Dirac 
delta-function  can  be  used  to  describe  its  distribution.  But 
this  difference  can  result  in  a  significant  difference  in  the 
property  of  a  matrix  (A)  associated  with  the  received  inter¬ 
ference  power.  A  can  be  positive  definite  or  positive  semidefi- 
nlte.  In  reality,  there  is  no  true  point  source  in  the  mathe¬ 
matical  sense,  and,  therefore,  the  division  of  these  two  types 
of  Interferences  is  somewhat  artificial  and  not  sharply  defined. 
Numerical  results  for  these  Interferences  can  be  meaningfully 
Interpreted,  only  if  the  computer  accuracy  and  system  noise  are 
knotm  beforehand.  This  problem  is  addressed  here. 


For  point  Interference  sources.  It  Is  well  known  that  the 


array  should  be  designed  to  have  nulls  in  all  the  Interference 
directions.  In  a  linear  array  with  equally  spaced  elements, 
this  is  easily  accomplished  by  imposing  certain  zeros  of  the 


polynomial  which  represents  the  array  pattern  function.  This 

2 


approach  is  used  by  Drane  and  Mcllvenna  ,  who,  in  addition, 
maximized  the  directivity  (D)  with  the  remaining  degrees  of 
freedom  of  the  problem.  However,  for  a  more  general  case,  such 
as  a  two-  or  three-dimensional  array  with  nonuniform  spaclngs 
and  with  the  Interferences  of  both  types  simultaneously  distri¬ 
buted  in  space,  this  problem  becomes  much  more  complicated. 

The  objective  of  this  paper  is  to  provide  solutions  to  these 
problems  with  a  few  typical  numerical  examples. 


2.  FORMULATION 

The  problem  discussed  here  is  that  of  finding  the  excita¬ 
tions  (J>)  for  an  antenna  array,  such  that  the  signal  to  Inter¬ 
ference  power  ratio  (SIR)  is  optimum.  The  equivalent  problem  of 


finding  the  optimum  signal-to-nolse  power  ratio  has  already  been 

I 


formulated  by  Lo  et  al.  Their  formulation  will  be  followed 
here.  In  addition,  there  are  cases  In  which  the  Interference 
can  be  reduced  to  zero  (or  SIR  =  «),  and  then  It  Is  obviously 
desirable  to  find  the  J>,  for  which  the  directivity  (D)  is  maxi¬ 
mum,  with  the  condition  that  SIR  -  «.  (It  will  be  seen  later 
that.  In  practice,  SIR  ■  <»  may  mean  that  It  Is  numerically 
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several  orders  of  magnitude  larger  than  all  other  quantities 


under  consideration.)  All  the  derivations  will  be  made  for  an 
arbitrary  array  and  for  any  arbitrarily  distributed  Interfer¬ 
ence,  which  can  be  coherent,  Incoherent,  or  a  combination 
thereof . 

Following  Lo  et  al.,  let  <J  =  (Jj ,  ...,  »  where  is 

the  excitation,  excluding  a  cophase  factor,  of  the  1th  antenna 
element  located  at  (x^,  y^ ,  z^).  The  array  factor, 

n 

P(0,(ji)  =  I  J.  exp  -  i|»°)  (1) 

1=1 

where 

=  x^  cos  ^  sin  0  +  sln(>  sln0  +  cos  0 
and  the  1th  element  cophasal  factor 

=  X  cos*  sln0  +  y  sin*  sln0  +  z.  cos  0 
1  1  o  o  1  o  o  1  o 

The  signal  direction  Is  specified  by  the  angles  0^  and 

Thus,  the  signal  power  la  proportional  to  <JCJ*>,  where  the 

product  defines  the  matrix  C  with  <Vj  =  (1,  1). 

it 

Indicates  the  complex  conjugate. 

SIR  -  <JCJ*>/<JAJ*> 

where  the  denominator  Is  proportional  to  the  Interference  power. 
For  Incoherent  Interference, 

<JAJ*>  «  |-  /  |P(0,t)|^|g(0,^)|^I(0,^)dn  (2) 

6 
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where  1(9,^)  is  the  magnitude  of  the  Interference  power  density. 
In  this  case, 

A  “  /  V*><VI(6,<>)|g(0,*)  l^dn  (3) 

4it 

where  <V  =  (expj  ( -  i|ij )  ,  expj(iji°  -  41^^)). 

3 

The  formulation  for  coherent  Interference  Is  similar. 

The  major  difference  between  coherent  and  Incoherent  Interfer¬ 
ences  Is  that  the  Integration  Is  performed  before  or  after 

A 

forming  the  product  V  ><V.  Ii  the  case  In  which  the  inter¬ 
ference  consists  of  both  types  or  several  coherent  regions  but 

Incoherent  among  the  regions,  Incorporating  them  all  Into  A  Is  a 

3 

straightforward  matter.  The  directivity  can  be  found  from  SIR 
by  setting  1(9,41)  “  I  for  all  9  and  For  this  case,  A  is 
denoted  by  B.  From  the  definition,  C,  A,  and  B  are  hermltlan. 

In  fact,  C  Is  positive  seraldeflnlte  and  B  Is  positive  definite. 

A  can  be  positive  definite  or  positive  semldef Inlte ,  depending 
on  the  Interference  distribution.  This  Is  explained  In  the  next 
chapter. 

3.  SOLUTION 

In  general,  the  Incoherent  Interference  distributions, 

which  will  be  considered  numerically  In  Chapter  4,  can  be 

classified  Into  the  following  four  different  types: 

1  —  Interference  over  a  non-vanlshlngly  small  solid  angle 
region; 


2  —  sliced  interference  source  from  to  ^2  with  a  fixed 

9,  or  from  9^  to  92  for  a  fixed 

3  —  point  interference  source;  and 

4  —  any  combination  of  1,  2,  and  3. 

Interferences  2  and  3  are  really  special  cases  of  1,  as  the  dis¬ 
tribution  in  one  of  the  angular  regions  becomes  a  delta  function. 
For  example,  interference  2  is  simply  type  1  with  a  delta  func¬ 
tion  distribution,  <S(9  -  9,^  ,)  or  S((p  -  .).  Similarly, 

fixed  fixed 

Interference  3  is  type  1  with  delta  functions  6( 9  -  9,,  ,) 

fixed 

Since  the  pattern  function  can  have  nulls,  for  interferences 
of  type  2  or  3,  it  is  possible  to  design  an  array  with  nulls  in 
the  directions  of  those  Interferences.  As  a  result,  for  these 
cases  SIR  can,  theoretically  speaking,  be  Infinity.  But,  In 
practice,  this  Implies  only  a  very  large  number.  Therefore,  an 
optimum  design  for  those  cases  should  be  one  with  directivity  in 
the  signal  direction  maximized,  while  the  nulls  in  the  prescribed 
interference  directions  are  maintained.  The  number  and  the 
directions  of  nulls  an  array  can  have  depend  on  the  number  of 
elements  as  well  as  on  the  array  geometry.  These  ideas  will  be 
set  forth  next  in  the  following  three  cases: 

Case  1 .  When  A  Is  positive  definite,  then  the  solution  of 
the  problem  of  optimizing  SIR  is  equivalent  to  that  of  optimizing 
the  slgnal-to-nolse  ratio.  This  problem  has  already  been  solved 


by  Lo  et  al.,  and  the  solution  will  be  used  here.  J>  for  optimum 
D  (DO)  is 

J>  =  Vj>  (3) 

and  J>  for  optimum  SIR  (SIRO)  Is 

J>  =  Vj>  (4) 

Case  2.  A  Is  positive  definite,  but  some  eigenvalues  ( X) 
of  A  are  so  small,  as  compared  with  other  eigenvalues,  that  one 
can  treat  A  as  if  It  were  positive  semideflnlte  within  the 
working  precision.  In  doing  so,  the  method  developed  below  for 
case  3  can  be  applied. 

Case  3.  When  interference  type  2  and/or  3  is  present,  then 
A  can  be  positive  semideflnlte  and  the  Interference  could  be 
nulled  out.  As  stated  above,  for  practical  purposes,  this  case 
and  case  2  may  be  found  equivalent  numerically.  This  is  dis¬ 
cussed  next. 

3.1  Matrix  A  Positive  Semideflnlte 

Let  A  be  positive  semideflnlte.  This  implies  that  <JAJ*> 
or  equivalently  <J*A*J>  can  equal  zero  for  some  nontrivial  J>. 
This  means  that  there  exists  a  subspace  of  excitation  vectors 
for  which  SIR  =  «•.  It  will  now  be  shown  that  the  eigenvectors 
of  A*  which  correspond  to  the  zero  eigenvalue  comprise  the  basis 
vectors  of  the  homogeneous  solution  space  of  <J*A*J>  ■  0, 


Since  A  is  hermitian  and  positive  semidef Inite ,  it  has 

*  * 

real  non-negative  eigenvalues,  2.  Thus,  <J^A 

it  it 

*  ^'^i^i'^i^  ”  ^i*''^i'^i^’  eigenvector  corresponding 

* 

to  Xj^.  For  a  nontrivial  X^<J  J^>  ■  0,  if  and  only  if  X^^  • 

*  it 

0.  The  dimension  of  the  null  space  of  <J  A  J>  equals  the  multi¬ 
plicity  of  the  zero  eigenvalue  X^.  (In  practice,  one  has  to 
decide  on  an  upper  bound  below  which  the  eigenvalues  could  be 
treated  as  zero.  This  upper  bound  could  be  10  with  q  set 
according  to  the  system  noise.)  But  all  the  eigenvectors  are 
Independent.  Consequently,  the  eigenvectors  with  corresponding 
eigenvalues  of  value  zero  form  a  basis  of  the  null  space.  Let 
the  columns  of  the  matrix  P  be  these  eigenvectors.  Now,  any 
linear  combination  of  these  vectors  is  also  a  solution  to 

it  it  it 

<J  A  J>  “  0  or  <JAJ  >  •  0.  Thus,  one  can  let 

J>  ■=  Py>  (5) 

where  y^ ,  the  1th  element  of  y>,  is  a  complex  constant  by  which 
the  eigenvector  in  the  1th  column  of  P  is  multiplied.  Thus,  it 
is  meaningful  to  define  the  optimum  J>  to  be  that  which  maxi¬ 
mizes  D  under  the  constraint  of  SIR  >  <»,  i.e.,  for  J>  in  the 
subspace  spanned  by  P.  In  other  words,  the  objective  is  to 
determine  y>  which  maximizes 


<yKy  > 

(6) 

where 

G  -  P^CP* 

(7) 

K  -  P^'BP* 

(8) 
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Both  matrices  K  and  G  are  herraltlan,  and  K  is  positive  definite 


i 


I 


f: 


r. 

I 


■ 


The  details  of  obtaining  the  solution  are  identical  to  those 


already  worked  out  by  Lo  et  al.,  for  the  similar  problem  of  opti¬ 


mizing  D  without  constraints.  The  results  are 


*-l  *L 

optimum  J>  =  PK  P 


(9) 


and 


D  =  <V^pV^S^> 


(10) 


where  <V  ^  =  (1 ,  . . .  ,  I ) . 


4. 


DISCUSSION  OF  RESULTS 

In  the  following  discussion,  A  with  at  least  one  eigenvalue 


of  the  order  of  magnitude  of  the  rounding  unit  (for  example, 
■14 


10  ,  for  the  Cyber),  will  be  called  positive  semidefinlte. 

Otherwise,  A  will  be  called  positive  definite.  The  computed  and 
graphed  absolute  values  of  the  array  pattens  for  various  inter¬ 
ferences  along  with  the  computed  values  of  the  SIRs  and  Ds  illus¬ 
trate  the  capabilities  and  limitations  of  the  theory  discussed 
in  Chapter  3.  The  arrays  considered  are  in  the  xy-plane.  The 
array  elements  are  short  monopoles  (sin9  pattern)  over  a  perfect 
ground  plane  and  parallel  to  the  z-axls.  In  addition,  they  are 
spaced  uniformly  with  spaclngs  X  along  the  x-axls  and  Y  along 
the  y-axls,  all  measured  in  the  wavelength  X.  The  spaclngs 
between  elements  have  been  kept  larger  than  0.5X  to  avoid  the 
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problems  of  supergain. ^  The  closed-form  solution  for  the  ele¬ 
ments  In  B  has  already  been  obtained  by  Lo  et  al. 


b  =  b  =  expj(-ij)°  ) 
nm  mn  2  nm 

for  n  <  m 


sin  2irp  cos  2np  sin  2iip 

_ nm  _ nm  _ nm 

2Trp  \2  f  1  .3 

nm  (2iip  )  (2  up  ) 

nm  nm 


(11) 


1 


(12) 


where 

o 

4’  " 

nm 

2itp 

nm 

sin  0  cos(  i  -  a 

o  0  nm 

(13) 

P 

nm 

[(x  - 

n 

.2  ,  n2,1/2 

+  (y„  -  yj  ) 

m  n  m 

(14) 

.  -1 
a  =  tan 
nm 

((yn- 

y)/(x  -  x)],0<a  <it 

m  n  m  —  nm 

(15) 

The  subscripts  n  and  m  are  referring  to  the  nth  and  mth  antenna 
elements,  respectively.  The  elements  of  A  are  found  by  numerical 
Integration.  Two  types  of  patterns  were  plotted  for  each  array. 
The  first  type  represents  patterns  for  an  array  excited  to  give 
the  optimum  SIR  (SIRO)  (Equation  (4),  for  A  being  positive  def¬ 
inite),  or  the  optimum  D  (DO)  (Equation  (5),  for  A  being  positive 
seraidef Inite )  with  the  constraint  that  the  received  Interference 
power  Is  at  a  numerically  very  small  value.  The  second  type 
represents  patterns  of  the  same  array  but  excited  for  maximum  D 
only  (Equation  (3)).  Although,  some  arrays  are  subjected  to 
different  combinations  of  interference  distributions,  only  one 
pattern  for  the  unconstrained  DO  needs  to  be  plotted.  To  pre- 


vent  confusion,  these  plots  do  not  have  interfrence  distribu¬ 
tions  drawn  on  them.  Unless  otherwise  noted,  the  plotted 
patterns  are  all  in  the  xy-plane,  i.e.,  9  =  90  and  all  angles 
are  in  degrees.  To  show  the  "nulls"  in  detail,  the  central  por¬ 
tion  is  expanded  in  all  plots  and  the  expansion  factor  can  be 
seen  from  the  scales  in  the  plots.  The  signal  (S)  direction  has 
always  been  set  along  the  positive  x-axls,  i.e.,  9^  =  90  and 
=  0.  The  signal  power  has  been  assumed  to  equal  unity,  always. 
For  comparison,  SIR  and  D  for  uniform  excitation  (U)  have  also 
been  calculated  for  each  array. 

The  information  for  all  the  examples  considered  is  sum¬ 
marized  in  Figures  I  -  3.  (All  figures  have  been  placed  after 
Chapter  5.)  When  following  the  discussion  for  Figures  4  -  15  to 
be  given  later,  the  use  of  these  figures  is  recommended.  In 
these  figures,  sketches  of  patterns  for  SIRO  or  for  DO  with 
constraint  (depending  on  which  is  applicable)  are  Included.  All 
the  data  under  the  heading  1  describe  the  type  of  Interference 
the  particular  array  has  been  subjected  to.  For  example,  in 
Figure  1,  in  the  first  column  under  I,  (2-1):  2  and  (2-2):  1 
stand  for  type  2  distributions  of  Intensities  2  and  1,  respec¬ 
tively,  (i.e.,  they  are  2  and  I  times  stronger  than  the  signal). 
The  numbers  1  and  2  after  the  dash  are  used  to  denote  the 
regions  of  interference  which  are  specified  by  the  angles 
following  each  Intensity  value.  There  are  M  rows  of  array 


elements  with  N  elements  each  parallel  to  the  y-axls.  The  rest 
of  the  notations  used  In  Figures  1-3  will  become  clear  in  the 
discussion  below. 


4.1  Matrix  A  Positive  Definite 

In  order  to  demonstrate  the  versatility  of  the  theory,  a 
2  X  10  array  subjected  to  all  three  types  of  interference  dis¬ 
tributions  (1,2,3)  is  considered.  The  distributions  with  the 
results  are  summarized  in  Figure  1(a).  It  turns  out  that  SIRO 
is  very  good,  but  D  for  SIRO  is  two  orders  of  magnitude  smaller 
than  DO.  This  situation  is  changed  drastically  if  a  type  2 
interference  of  small  Intensity  and  wide  distribution  is  added 
in  (Figure  1(b)).  Now  D  for  SIRO  is  80.4  and  SIRO  =  945.5, 
which  is  still  much  higher  than  the  SIR  of  49.1  when  D  is  maxi¬ 
mized. 

4.2  Matrix  A  Positive  Semidefinite 

As  discussed  before,  for  some  Interference  distributions, 

A  can  have  zero  as  an  eigenvalue.  Since  all  computers  have  a 
finite  accuracy,  numerically,  how  small  a  quantity  can  be  re¬ 
garded  as  zero  is  a  very  relevant  question.  To  this  end,  one 
needs  to  decide  upon  the  upper  bound  (10  **)  below  which  the 
eigenvalues  of  A  will  be  treated  as  zero.  At  the  same  time,  one 
must  consider  that  generally  the  larger  the  null  space  of  A,  the 
larger  the  corresponding  DO  with  constraint,  as  derived  from 
Equation  (6).  The  choices  of  q  •  7,  8,  10  have  provided  good 


results  for  type  2  (or  3)  interferences  as  discussed  below. 
Clearly,  the  SIR  so  obtained  depends  on  q,  but  for  sufficiently 
large  q  the  dependence  can  become  Insignificant.  Numerically, 
the  resulting  solution  can  thus  be  considered  as  optimum  In  the 
sense  that  the  signal  Is  maximum  under  the  constraint  that  SIR 
is  kept  sufficiently  large.  As  before,  DO  without  constraint 
will  be  referred  to  simply  as  DO. 

Figure  7  depicts  the  pattern  for  DO  with  constraint  for  a 
2  X  10  array  subjected  to  two  type  3  Interferences.  On  the  one 
hand,  SIRO  =  10^^,  which  Is  much  higher  than  SIR  =  53.7  for  DO 
and  than  SIR  =  66.3  for  U  as  shown  In  Figure  2(a).  On  the  other 
hand,  DO  with  constraint  =  102.3  compares  favorably  with  DO  » 
138.3,  and  D  ®  112.3  for  U.  For  this  computation,  q  =  13 
resulted  in  18  columns  for  P  in  Equation  (5),  In  spite  of  the 
fact  that  DO  with  constraint  Is  optimum,  the  asymmetry  of  the 
Interference  distributions  was  apparently  sufficient  to  shift 
the  main  beam  a  little  off  the  direction  of  the  signal. 

In  the  following  example,  the  element  spacing  of  a  3  x  20 
array  Is  assumed  to  be  0.6x.  A  type  2  Interference  of  intensity 
1000  Is  distributed  from  ^  “  175  to  185  (centered  on  the  nega¬ 
tive  x-axls)  In  the  0-90  plane  as  shown  In  Figures  2(b)  and  8. 
Thus,  the  signal  and  the  type  2  distribution  are  in  the  same 
plane.  Because  of  the  symmetry  in  this  example,  the  pattern  is 


symmetric  as  seen  In  Figure  8.  The  Inference  Is  virtually 

nulled  out  completely  (SIR  =10  ).  A  clearer  view  of  this  Is 

shown  In  the  magnified  center  portion  of  the  pattern  In  Figure 

8.  D  with  constraint  Is  214.6,  a  value  that  approaches  DO  = 

222.1  (Figure  2(b)),  leaving  D  =  144.6  for  U  far  behind.  Again, 

the  SIRs  for  DO  and  U  (8.3  and  0.6)  are  miniscule  by  comparison 
12 

with  SIRO  of  10  .  For  this  example,  q  =  10  with  a  resulting  51 

columns  for  P. 


One  may  be  Inclined  to  think  that  the  high  SIRO  in  the  pre¬ 
vious  example  was  only  possible  because  the  type  2  interference 
was  highly  confined.  This  is  not  the  case.  Figure  10  shows  the 
pattern  for  DO  with  constraint  for  a  3  x  22  array  that  is  sub¬ 
jected  to  a  type  2  Interference  distributed  over  the  entire 
second  and  third  quadrants  In  the  0  =  90  plane  as  indicated  in 
Figure  3(a).  It  turns  out  that  SIRO  Is  more  than  ten  million 
times  larger  than  the  SIRs  for  DO  and  U,  while  the  Ds  for  all 
three  excitations  are  on  the  same  order  of  magnitude  [Figure 
3(a)).  For  this  case,  q  =  8,  resulting  in  19  columns  for  P. 


As  a  final  example,  we  consider,  for  academic  interest,  a 


3  X  22  array  in  an  environment  with  type  2  Interference  of  an 
.3 


Intensity  10  times  stronger  than  that  of  the  signal,  and  dis¬ 
tributed  in  the  region  0  ■  -5  to  +5  at  0  *  80.  Recall  that  the 
signal  is  always  assumed  to  be  in  the  ^  =  0  and  0  =  90  direction 
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which  close  to  the  Interference.  The  results  for  optimum 
'ilR  as  well  as  those  for  optimum  D  are  summarized  In  Figure  3(b) 


with  detailed  patterns  shown  In  Figures  12  and  13.  As  It  should 
.9 


be,  SIRO  =  10  is  ten  orders  larger  than  the  SIR  of  0.26  for  the 
DO  case  and  0.25  for  the  U  case.  However,  the  D  associated  with 
SIRO  Is  only  2.9,  substantially  smaller  than  DO  =  481.6  and  D  = 
346.5  for  the  uniform  excitation.  It  is  of  particular  Interest 
to  examine  the  patterns  in  the  0  =  90  plane  for  two  cases:  SIRO 
and  DO  as  shovm  in  Figures  12  and  13.  For  the  former,  a  large 
lobe  (actually  the  main  beam)  is  directed  in  the  interference- 
free  region  in  the  left  half  space,  and  in  contrast  a  much 
smaller  lobe  is  in  the  signal  direction.  This  paradoxical 
result  can  be  explained  by  the  fact  that  the  integrals  in  SIR, 
being  weighted  by  the  zero  interference  Intensity  In  the  left 
half  space,  can  still  result  In  a  very  large  value  of  SIR.  For 
Figure  13  where  D  is  optimum,  the  main  beam  is  pointed  to  the 
signal  and  also  to  the  nearby  interference,  thus  resulting  in, 
not  surprisingly,  a  high  value  of  D  but  a  very  small  value  of 
SIR.  This  Is  more  clearly  seen  in  Figures  14  and  15  which  are 
the  patterns  in  the  conical  cut  with  9  ••  80.  The  main  beam  is 
practically  lying  completely  in  the  interference  region.  One 
should  not  be  surprised  by  these  results  since  optimizing  SIR  is 
a  problem  totally  different  from  that  for  optimizing  D,  and  one 
does  not  necessarily  imply  the  other,  particularly  when  the 
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signal  and  interference  are  close  to  each  other.  However,  to 
differentiate  them  more  effectively,  a  larger  array  should  be 
used.  On  the  other  hand,  to  avoid  a  large  back  lobe  shown  in 
Figure  12,  a  background  interference  (or  noise)  in  the  left  half 
space  should  be  assumed. 


In  the  above  study,  q  was  chosen  to  be  7  which  results  in 
57  columns  for  P.  One  may  wonder  what  would  happen  to  the 
results  if  q  is  decreased.  It  turns  out  that  a  choice  of  q  =  2 
increases  the  number  of  columns  in  P  to  60.  This  increase  pro¬ 
duces  a  marginal  change  in  DO  with  constraint,  from  2.9  to  2.95, 

9 

but  reduces  SIRO  dramatically  from  10  to  3277,  Thus,  for  a 
good  approximate  solution  to  SIRO,  as  stated  earlier,  a  suf¬ 
ficiently  large  q  should  be  used  in  order  to  obtain  a  better 
"null"  space  so  that  the  Interference  contribution  can  be  kept 
small . 


One  may  also  wonder,  whether  or  not  it  is  feasible  to  com¬ 
pute  accurately  the  eigenvectors  that  comprise  the  columns  of  P. 
It  may  not  be,  considering  that  the  corresponding  eigenvalues 
are  very  small.  Fortunately,  one  can  compute  accurately  the 
eigenvectors  associated  with  large  eigenvalues,  say  larger  than 
10  Then  the  null  space  can  be  accurately  determined  from  the 
vectors  orthogonal  to  these  eigenvectors. 


4.3  Sensitivity  Test  for  Optimum  Excitations 

The  sensitivity  of  the  optimum  excitation  giving  SIRO 

was  tested  for  the  three  examples  shown  in  Figures  2(a),  2(b), 
and  3(a)  by  injecting  a  random  error  uniformly  distributed  be¬ 
tween  -3%  to  +3%  Into  each  component  of  real  and  Imaginary 

parts,  independently.  The  perturbed  SIRO  Is  referred  to  as 
PSIRO.  For  completeness,  J>  for  DO  was  perturbed  similarly, 
resulting  In  a  perturbed  DO  (PDO).  The  results  are  shown  in 
Table  1.  Clearly,  the  PSIROs  are  very  much  smaller  than  the 
SIROs,  Indicating  that  the  astronomically  large  value  of  SIRO  In 
all  the  examples  should  not  be  taken  seriously.  However,  the 
PSIROs  are  still  much  larger  than  the  SlRs  for  PDO  or  DO.  Thus, 
the  ideal  values  of  SIRO  may  be  unobtainable  in  reality,  but  the 
optimum  excltatlni  Jq>  can  still  give  a  much  superior  SIR  as 
compared  with  otter  excitations.  As  expected,  the  random  errors 
do  not  affect  ot  ler  quantities  significantly  because  of  the 
large  element  spaclngs  used. 

5.  CONCLUSION 

An  optimization  technique  for  arbitrary  arrays  In  arbitrary 
interference  environments  has  been  developed. 


The  interference  environments  can  be  generally  separated 
into  two  kinds:  those  consisting  of  interference  distributions 
that  cannot  be  nulled  out  completely,  and  those  that  can.  For 
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TABLE  1.  Results  of  Perturbing  the  Real  and  Imaginary 

Components  of  J>  for  SIRO  and  DO  by  a  Uniformly 
Distributed  Error  of  a  Maximum  of  ±  3%. 


SIR  (Figure 

2(a)) 

D 

SIRO 

,o*>  3 

102.3 

PSIRO 

1.7  X  10c 

101.9 

3.2  X  10  c 

102.6 

15.1  X  10 

102.1 

PDO 

53.6 

138.2 

54.1 

138.2 

54.0 

138.2 

SIR  (Figure 

2(b)) 

n 

SIRO 

10^2 

214.6 

PSIRO 

7930.0 

214.3 

2510.0 

214.5 

4540.0 

214.5 

PDO 

7.8 

221.9 

8.3 

221.9 

8.5 

222.0 

SIR  (Figure 

3(a)) 

D 

SIRO 

10^ 

364.4 

PSIRO 

3399.0 

363.6 

5607.0 

364.4 

3848.0 

364.3 

PDO 

35.1 

482.8 

34.6 

482.9 

35.6 

482.9 

the  first  kind,  the  optimum  excitation  for  maximum  signal-to- 
interference  ratio  (SIR)  can  be  determined  from  the  conjugate 
Inverse  of  a  positive  definite  matrix  A  associated  with  the 
Interference  distribution.  For  the  second  kind,  since  the 
Interference  can  be  nulled  out,  the  optimum  excitation  Is 
defined  to  be  one  which  maximizes  D  (l.e.,  the  desired  signal) 
under  the  constraint  that  the  Interference  remains  nulled  out. 
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The  method  used  consists  of  first  finding  the  null  space  eigen¬ 
vectors  of  a  positive  semldeflnlte  matrix  A,  and  then  determin¬ 
ing  the  excitation  vector  In  this  space  which  maximizes  D.  In 
practice,  the  type  of  Interference  that  can  be  "nulled  out" 
(l.e.,  brought  below  a  prescribed  Insignificantly  low  level) 
need  not  be  limited  to  point  sources,  but  can  be  distributed 
over  a  sufficiently  small  solid  angle.  Thus,  at  the  heart  of 
the  method  described  here  lies  a  decision  about  which  eigenvalue 
of  A  can  be  treated  as  If  it  were  zero.  This  decision  should  be 
made  based  on  the  overall  system  noise.  Examples  for  inter¬ 
ference  distributions  of  various  types  have  been  examined 
numerically.  The  sensitivity  of  SIR  due  to  random  errors  in  the 
excitations  has  been  tested,  as  well.  The  results  show  clearly 
that  the  optimization  of  SIR  and  the  optimization  of  D  are  in 
general  incompatible.  In  fact,  the  maximized  SIR  can  be  many 
orders  of  magnitude  larger  than  the  SIR  obtained  if  only  D  is 
maximized.  Therefore,  the  criterion  for  designing  a  receiving 
array  is  different  from  that  for  designing  a  transmitting  array. 
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Figure  3  (a),  (b).  Summary  sheet  pertaining  to  radiation 
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f^adlacton  p«eccrn  of  a  2  «  tO  array  excltad  for  opclaun  SIR. 

For  daeatla  of  the  Incerferanca  dlscrtbucion,  I,  rafer  to  Figure  1(a) 
The  center  portion  le  expanded  at  the  upper  right. 
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Abstract 


n 

Antennas  to  provide  monopulse  direction-of-arrival  in  a  dense 
multi-polarized  signal  environment,  over  a  wide  instantaneous  RF 

m 

bandwidth,  are  a  challenge  to  the  designer.  Solutions  to  some  of  the 
associated  problems  are  presented,  along  with  a  description  and  per¬ 
formance  data  for  a  wideband  monopulse  antenna  system. 

The  monopulse  reflector  antenna  described  uses  a  unique  array  feed 
configuration  to  reduce  the  apparent  phase  center  separation,  thereby 
increasing  bandwidth,  along  with  absorber  techniques  to  control  the 
illumination  function  and  elevation  beamwidth.  Primary  and  secondary 
patterns  of  the  DF  antenna,  along  with  sidelobe  blanking  antenna  data, 
are  presented  for  an  operational  system.  Photographs  and  diagrams  of 
the  antennas  are  included.  This  antenna  configuration  is  currently 
deployed  as  part  of  an  ESM  feasibility  demonstration  system. 
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Techniques  for  monopulse  direction  finding  (DF)  have  been  used  for 
many  years  in  radar  systems.  Monopulse  allows  almost  instantaneous 
determination  of  direction  of  arrival  based  on  a  single  pulse  rather 
than  a  train  of  pulses.  A  narrow  RF  bandwidth  of  the  antenna  is  ade¬ 
quate  for  most  radar  applications  because  only  a  narrow  frequency 
range  is  needed  for  tracking.  Electronic  surveillance,  on  the  other 
hand,  requires  antennas  that  will  "listen"  over  a  wide  RF  bandwidth. 
Wideband  monopulse  direction  finding  antennas  for  surveillance  applies 
tions  present  some  challenging  problems  for  the  antenna  designer. 
These  problems  and  some  solutions  are  presented  here. 

2.0  DESIRABLE  FEATURES  OF  A  WIDEBAND  MONOPULSE 
SYSTEM. 


Desirable  features  of  a  wideband  monopulse  surveillance  antenna  sys¬ 
tem  are  as  follows.  The  antenna  beamwidth  must  be  large  enough  to 
cover  an  adequate  instantaneous  field  of  view  (FOV),  while  the  gain 
must  be  high  enough  to  provide  sufficient  system  sensitivity.  High  DF 
accuracy,  usually  the  most  important  characteristic  of  a  monopulse 
system,  is  achieved  through  careful  selection  of  pattern  shapes  and 
processing  techniques.  Since  the  incoming  signal  polarization  is  usually 
unknown,  circular  polarization  with  a  low  axial  ratio  is  normally 
required  so  that  the  DF  accuracy  is  not  a  function  of  the  incoming 


signal  polarization.  Also,  ambiguities  due  to  reception  of  unknown 
signals  through  the  sidelobes  must  be  eliminated,  usually  through  the 
use  of  auxili8u*y  sidelobe  blanking  antennas.  All  of  these  features 
should  ideally  be  obtained  over  a  sufficiently  wide  RF  bandwidth  to 
meet  the  surveillance  requirements. 


3.0  MONOPULSE  FEEDS  EXCITATION  TECHNIOUE 


Often  a  reflector  antenna  is  required  to  provide  high  gain  for  system 
sensitivity  and  a  limited  field  of  view  for  spatial  filtering  in  a  dense 
signal  environment,  rather  than  a  "wide  open**  DF  antenna  with  a  wide 
beamwidth  and  low  gain.  Assuming  a  parabolic  reflector  properly  sized 
for  the  desired  beamwidth,  the  critical  items  become  types  of  feeds 
and  methods  of  excitation  for  efficient  reflector  illumination  and 
optimum  DF  accuracy. 


Figure  1  shows  two  schemes  for  determining  direction  of  arrival  (DOA) 
in  one  plane  by  two-channel  comparison.  For  left-right  beams,  the 
amplitude  ratio  of  the  two  beam  outputs  is  a  function  of  the  coordi¬ 
nate  angle,  and  the  change  in  the  amplitude  ratio  versus  angle  is  rela¬ 
tively  constant  across  the  FOV.  For  sum-difference  beams,  the 
amplitude  ratio  of  the  two  outputs  is  a  function  of  the  coordinate 
angle,  and  their  phase  difference  indicates  left  or  right  of  the 
boresight  null.  Because  of  the  steep  difference  null,  however,  the 
change  in  the  amplitude  ratio  versus  coordinate  angle  increases  sharply 


near  boresight,  thus  providing  higher  DF  accuracy  than  the  left-right 
scheme. 

Since  the  feeds  are  used  individually  rather  than  in  an  array  in  the 
left-right  scheme,  use  of  frequency  independent  feeds  can  provide  fre¬ 
quency  independent  reflector  illumination.  In  the  sum-difference  array, 
however,  the  sum  pattern  illumination  is  narrower  than  the  difference 
pattern,  and  the  array  factor  further  narrows  the  illumination  patterns 
as  the  frequency  is  increased.  Thus,  although  the  sum-difference 
scheme  provides  higher  DF  accuracy  near  boresight,  which  is  advanta¬ 
geous  for  tracking  applications,  achieving  efficient  reflector  illumi¬ 
nation  over  a  broad  frequency  range  is  more  difficult  than  with  the 
left-right  scheme. 


angle  can  be  placed  such  that  the  base  separation  in  wavelengths  at 
the  lowest  operating  frequency  is  the  same  as  the  tip  separation  in 
wavelengths  at  the  highest  operating  frequency,  thus  providing  constant 
reflector  illumination  over  the  entire  spiral  frequency  range.  However, 
since  the  phase  center  of  the  array  moves  along  the  length  of  the 
conical  spiral  as  frequency  changes,  axial  defocusing  of  the  reflector 
results,  causing  a  loss  of  efficiency.  The  major  drawback  with  the  con¬ 
ical  spiral  approach  is  the  interaction  or  coupling  of  the  spirals  when 
placed  close  together,  which  results  in  deterioration  of  the  array  pat¬ 
tern  at  various  frequencies  across  the  band.  Coupling  between  adjacent 
flat  spirals  is  much  less  and  has  no  significant  effect. 

Thus,  while  conical  spirals  appear  to  be  well  suited  for  a  broadband 
monopulse  feed  array,  considerable  experimentation  is  required  to  solve 
the  coupling  problem,  if  it  can  be  solved!  On  the  other  hand,  a  tech¬ 
nique  that  makes  the  phase  centers  "appear"  closer  together  could 
cause  a  flat  spiral  array,  or  an  array  of  horns  or  some  other  elements, 
to  be  useful  as  a  feed  over  a  reasonably  broad  frequency  band.  Such  a 
technique  is  presented  in  paragraph  6.0. 

5.0  TWO-ELEMENT  FEED  ARRAY. 


A  simple  two-element  feed  array,  as  shown  in  Figure  3,  has  limited 
use  as  a  broadband  feed  as  explained  in  paragraph  4.0.  As  shown  in 
Figures  4  through  6,  the  sum  pattern  narrows  rapidly  as  frequency 


increases.  The  secondary,  or  reflector,  patterns  in  Figures  7  and  8 
show  a  flattening  and  then  a  dip  in  the  sum  pattern  caused  by  the 
first  nulls  and  first  sidelobes  of  the  feed  pattern  moving  onto  the 
reflector.  Measured  data  indicates  that  a  pair  of  contiguous  flat  spiral 
or  circular  horn  feeds  in  the  sum-difference  mode  will  illuminate  a 
reflector  over  a  useful  bandwidth  ratio  of  about  1.6:1. 


6.0  THREE-ELEMENT  FEED  ARRAY. 


A  technique  that  extends  the  bandwidth  of  a  reflector  with 
sum-difference  feeds  of  fixed  spacing  is  shown  in  Figure  9.  The  pur¬ 
pose  of  the  third  element  is  to  cause  the  element  spacing  S  to  appear 
smaller  for  the  sum  mode,  thus  broadening  the  sum  illumination  pat¬ 
tern,  while  leaving  the  spacing  unchanged  for  the  difference  mode, 
which  is  already  sufficiently  broad. 


The  network  in  Figure  9  is  easiest  to  analyze  in  the  transmit  mode, 
and  since  the  circuit  elements  and  antennas  are  reciprocal,  the  analy¬ 
sis  applies  also  to  the  receive  mode.  A  signal  injected  into  the  sum 
port  is  split  equally  between  the  two  circuit  halves.  In  the  left  half,  it 
is  split  again  with  half  going  to  antenna  B  and  half  going  to  the 
hybrid  that  feeds  antenna  A.  The  signal  in  the  right  half  is  similarly 
split  between  antenna  C  and  the  hybrid  that  feeds  A.  Since  the  signals 
into  the  hybrid  that  feeds  A  are  of  equal  amplitude  and  phase,  all  of 
the  energy  appears  in  antenna  A  (neglecting  hybrid  losses)  and  no 


energy  appears  in  the  difference  port  load.  Thus  for  the  sum  mode, 
the  circuit  elements  A,  B,  and  C  are  in  phase  and  have  a  1:2:1  power 
ratio,  respectively. 

A  signal  injected  into  the  difference  port  is  divided  in  exactly  the 
same  way,  except  that  all  of  the  energy  into  the  hybr  d  feeding  A 
appears  in  the  difference  port  load  with  no  energy  going  into  antenna 
A,  due  to  the  inputs  to  that  hybrid  being  of  equal  amplitudes  and 
opposite  phase.  Thus  the  difference  mode  excites  elements  B  and  C 
only  with  equal  amplitudes  and  180  degrees  phase  difference. 

Referring  to  the  triangular  array  in  Figure  9,  note  that  for  the  sum 
mode,  the  energy  in  B  combines  with  half  the  energy  in  A  resulting  in 
a  center  of  phase  half  way  between  elements  B  and  A.  Likewise,  the 
energy  in  C  combines  with  the  other  half  of  the  energy  in  A  resulting 
in  a  center  of  phase  half  way  between  elements  C  and  A.  These  two 
subarrays,  with  phase  centers  spaced  S/2  apart,  are  now  arrayed  to 
produce  a  sum  mode  pattern  equivalent  to  that  of  a  two-element  array 
of  spacing  S/2.  The  difference  mode  pattern  is  that  of  a  two  element 
array  of  spacing  S,  since  the  power  in  element  A  is  zero  for  differ¬ 
ence  mode  excitation. 

Patterns  of  an  array  of  three  flat  cavity-backed  spirals  fed  as  in  Fig¬ 
ure  9  are  shown  in  Figures  10  through  12.  Note  that  the  sum  and  dif¬ 
ference  patterns  have  approximately  equal  widths  across  a  large 


frequency  band,  compared  with  the  two-element  array  patterns  of  Fig¬ 
ures  4  through  6,  where  the  sum  patterns  are  much  narrower  than  the 
difference  patterns.  Secondary  patterns  for  the  three-spiral  array  in 
Figures  13  through  16  show  no  flattening  of  the  sum  pattern  over  a 
full  2.3:1  bandwidth  ratio.  Clearly,  the  three-element  array  extends  the 
useful  bandwidth  of  the  monopulse  reflector  by  making  the  "apparent" 
sum  mode  element  spacing  less  than  the  physical  spacing. 

7.0  FOUR-ELEMENT  FEED  ARRAY. 

In  a  monopulse  system,  the  sum  and  difference  beams  should  ideally 
point  along  the  same  boresight.  Analysis  of  the  triangular  array  in 
Figure  9  reveals  that  for  the  sum  mode  the  overall  phase  center  is 
located  half  way  up  the  vertical  centerline  of  the  triangle,  while  the 
difference  mode  overall  phase  center  is  at  the  midpoint  of  the  triangle 
base.  This  vertical  phase  center  offset  causes  the  reflector  sum  and 
difference  pattern  boresights  to  be  slightly  offset  in  the  vertical  plane. 
A  solution  to  this  problem  is  the  four-element  feed  array  shown  in 
Figure  17.  Note  that  the  only  difference  from  the  three-element  array 
is  that  the  energy  that  went  into  element  A  is  now  split  equally 
between  elements  A  and  D.  While  the  phase  center  separations  in  the 
horizontal  direction  are  still  S/2  for  the  sum  mode  and  S  for  the  dif¬ 
ference  mode,  the  overall  array  phase  centers  are  coincident  with  the 
array  physical  center  for  both  sum  and  difference  beams,  thus  produc¬ 
ing  colinear  boresights. 


Patterns  of  a  four-element  array  are  shown  in  Figures  18  through  20. 
The  array  elements  are  circular  waveguide  horns  with  dielectric  phase 
plate  polarizers.  Note  the  similarity  to  the  three-element  array  pat¬ 
terns  of  Figures  10  through  12.  The  reflector  patterns  in  Figures  21 
through  23  show  that  the  four-element  array  produces  uniform  second¬ 
ary  patterns  with  no  flattening  of  the  sum  pattern  over  a  bandwidth 
ratio  of  at  letist  2.7:1. 

Pictures  of  four-element  feed  array  networks  in  both  waveguide  and 
coaxial  configurations  are  shown  in  Figure  24.  The  waveguide  network 
uses  the  "magic  tee"  as  the  circuit  equivalent  to  the  180  degree 
hybrid. 

8.0  ELEVATION  PLANE  CONSIDERATIONS. 

In  the  foregoing  examples,  direction  finding  is  performed  in  the 
azimuth  plane  only,  and  the  elevation  plane  beamwidths  are  sufficient¬ 
ly  broad  to  provide  an  adequate  vertical  field  of  view.  Since  signals 
may  be  incident  from  any  angle  in  the  elevation  beam,  the  sum  and 
difference  patterns  should  ideally  be  independent  of  elevation  angle. 
This  implies  equal  sum  and  difference  elevation  beamwidths,  so  that 
both  beams  roll  off  at  the  same  rate  in  elevation,  thus  providing  a 
sum-difference  amplitude  ratio  versus  azimuth  function  that  does  not 
change  with  elevation.  Equal  sum  and  difference  elevation  beamwidths 


are  obtained  by  equal  illumination  tapers  across  the  vertical  dimension 
of  the  reflector. 

Figure  25  illustrates  illumination  contours  in  decibels  on  a  rectangular 
reflector,  produced  by  the  four-element  feed  array.  Note  that  in  Fig¬ 
ure  17,  elements  A  and  D  are  effectively  arrayed  with  combined  ele¬ 
ments  B  and  C  in  the  vertical  plane  for  the  sum  mode  (all  four 
elements  are  excited  equally),  while  for  the  difference  mode  the  verti¬ 
cal  plane  sees  only  the  equivalent  of  the  single  element  pattern  of  B 
or  C,  with  no  arraying.  As  a  result,  difference  pattern  illumination  is 
broader  than  sum  pattern  illumination  in  the  elevation  plane,  thus 
producing  narrower  difference  than  sum  pattern  secondary  beams. 

A  solution  to  the  elevation  illumination  problem  is  shown  in  Figure  26. 
The  proper  absorber  thickness  is  determined  experimentally  such  that 
the  difference  pattern  illumination  taper  is  effectively  increased  to 
equal  that  of  the  sum.  Figures  27  through  30  show  patterns  of  this 
configuration  at  a  fixed  frequency  for  various  angles  in  the  elevation 
beam.  The  sum  and  difference  pattern  shapes  are  seen  to  be  fairly 
constant  over  a  15  degree  elevation  sector. 


9.0  AXIAL  RATIO  CONSIDERATIONS. 

Careful  design  and  fabrication  of  monopulse  antennas  to  provide  low 
axial  ratios  over  the  entire  frequency  band  are  critical  to  accurate  DF 


performance.  The  effect  of  axial  ratio  on  DF  accuracy  is  illustrated  in 
Figure  31.  Unlike  other  errors,  the  angular  uncertainty  induced  by 
axial  ratio  cannot  be  calibrated  out  of  the  system. 

10.0  SIDELOBE  BLANKING. 

In  a  typical  monopulse  system,  signals  strong  enough  to  enter  the 
sidelobes  can  trigger  a  sum-difference  ratio  response  indistinguishable 
from  that  of  the  main  beam.  A  common  technique  to  remove  such 
ambiguity  is  to  use  one  or  more  omnidirectional  antennas  in  a  sidelobe 
blanking  scheme  as  shown  in  Figure  32.  The  sidelobe  blanking  antenna 
or  antennas  must  provide  a  level  greater  than  the  sum  pattern 
sidelobes  throughout  all  angular  regions  where  unwanted  signals  may 
appear. 

An  omnidirectional  antenna  that  provides  good  hemispherical  coverage, 
optimized  toward  the  horizon,  is  the  four-arm  conical  spiral  operated 
in  the  difference  mode.  Typical  patterns  are  shown  in  Figures  33  and 
34.  Except  for  a  narrow  null  directed  toward  the  zenith,  this  antenna 
produces  uniform  azimuth  patterns  for  elevation  angles  out  to  the 
horizon.  Variations  with  frequency  are  minimal  over  a  band  at  least  as 
wide  as  that  of  the  monopulse  reflector  described  in  paragraph  7.0. 


11.0  TYPICAL  WIDEBAND  MQNQPULSE  ANTENNA  SYSTEM. 


Figures  35  through  37  show  an  integrated  monopulse  system  that  uses 
the  four-element  feed  array  described  in  paragraph  7.0,  the  absorber 
technique  for  elevation  pattern  control  described  in  paragraph  8.0,  and 
omnis  for  sidelobe  blanking  of  the  type  described  in  paragraph  10.0. 
The  large  reflector  has  cavity-backed  flat  spiral  feeds  with  a  coaxial 
network,  while  the  small  reflector  uses  circular  horn  feeds  with 
dielectric  phase  plate  polarizers  and  a  waveguide  network.  A  four-arm 
conical  spiral  omni  for  each  frequency  band  can  be  seen  in  the  photo¬ 
graphs.  The  absorber  "shield"  around  the  periphery  of  the  small  reflec¬ 
tor  is  to  block  spillover  that  may  exceed  the  sidelobe  blanking  levels. 
This  antenna  system  provides  performance  similar  to  that  described  in 
paragraphs  7.0,  8.0  and  10.0  over  two  2.7:1  frequency  bands. 
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12.0  CONCLUSIONS. 

While  the  theory  of  monopulse  antennas  is  well  established,  the  practi¬ 
cal  problems  that  the  designer  of  a  wideband  monopulse  antenna  sys¬ 
tem  must  face  include  the  wideband  realization  of  efficient  reflector 
illumination,  high  DF  accuracy,  adequate  beamwidth  control,  low  axial 
ratios,  and  sidelobe  blanking.  Solutions  to  these  problems  and  resultant 
data  have  been  presented  which  involve  array  feed  configurations, 
absorber  techniques,  careful  attention  to  antenna  types  and  con¬ 
struction  details  for  low  axial  ratios,  and  selection  of  a  sidelobe  blank- 
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ing  system  matched  to  the  sidelobe  characteristics  of  the  DF  antenna. 
Although  each  wideband  monopulse  system  may  have  a  unique  set  of 
requirements,  the  design  considerations,  trade-offs,  and  development 
processes  bear  a  strong  thread  of  commonality. 
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Figure  1.  Monopulse  Feeds  Excitation  Techniques  Comparison 
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•  SFIRAL  FHASE  CENTER  SEFARATION  IS 
FIXED 

•  BANPWlOTN  IS  LIMITED  DUE  TO  NARROWHUO 
OF  ARRAY  FACTOR  BEAMWIOTH  AT  HIGH 
ENO  OF  BAND.  CAUSING  REFLECTOR 
UNOCRILLUMINATION 


Figure  2.  Z  -  h  Monopulse  Feed  Array  Types  Comparison 


Figure  8.  Two-Spiral  Array  Secondary  Patterns,  16  GHz 
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Figure  9.  Three-Eleraent  Feed  Array  Concept 


Figure  10.  Three -Spiral  Array  Primary  Patterns 
Azimuth  Cut.  8  GHz 
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Figure  31.  Effect  of  Axial  Ratio  on  DF  Accuracy 
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Figure  32.  Sidelobe  Blanking  Technique 


Figure  33.  Four -Arm  Conical  Spiral  Omni  Patterns 
Elevation  Gut 


Figure  34.  Four- Arm  Conical  Spiral  Omni  Patterns 
Azimuth  Cut,  60**  Elevation 
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MISSION 

of 

Rome  Air  Development  Center 

RAVC  p£a/ii  and  execate^  ^e^zoAck,  envelopment,  tei-t  and 
4cCe,ctcd  acijuc6AtAon  pAogAam6  f.n  ^uppoAt  Command,  CorvtAot 
CormurUcatioru  and  Iyit&lUgenc&  (C^I)  acZivOUiu.  TtchnLcal 
ayid  zngf.mzAing  iuppoAt  uiUJUn  ew-eew  o^  tecAnccat  competence 
x-4  pAQoidtd  to  ESO  VfiogKam  OUlau  {?06)  and  otkvi  ESV 
etemen^.  Tkc  p^xnetpot  technical  mciilon  aA.eai  cuie 
communAcaXioni ,  electAomagnetic  guidance  and  contact,  iuA- 
veAllance  ground  (uid  aeAoipace  objecti,  inteiUgence  data. 
co^eetton  and  handling,  In^ofmation  6yitem  technology, 
iolid  itate  4ciencei,  electAomagnetlcA  and  electronic 
AeltabiZity,  maintainability  and  compatibility. 
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